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THE  MISSION  OF  AGARD 


The  mission  of  AGARD  is  to  bring  togethei  the  leading  personalities  of  the  NATO  nations  in  the  fields  of 
science  and  technology  relating  to  aerospace  for  the  following  purposes: 

— Exchanging  of  scientific  and  technical  information; 

— Continuously  stimulating  advances  in  the  aerospace  sciences  relevant  to  strengthening  the  common  defence 
posture; 
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- Improving  the  coKiperation  among  member  nations  in  aerospace  research  and  development; 

— Providing  scientific  and  technical  advice  and  assistance  to  the  North  Atlantic  Military  Committee  in  the 
field  of  aerospace  research  and  development; 
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— Rendering  scientific  and  technical  assistance,  as  requested,  to  other  NATO  bodies  and  to  member  nations 
in  connection  with  research  and  development  problems  in  the  aerospace  field; 


— Providing  assistance  to  member  nations  for  the  purpose  of  increasing  their  scientific  and  technical  potential; 

— Recommending  effective  ways  for  the  member  nations  to  use  their  research  and  development  capabilities 
for  the  common  benefit  of  the  NATO  community. 

The  highest  authority  within  AGARD  is  the  National  Delegates  Board  consisting  of  officially  appointed  senior 
representatives  from  each  member  nation.  The  mission  of  AGARD  is  carried  out  through  the  Panels  which  are 
composed  of  experts  appointed  by  the  National  Delegates,  t.ic  Consultant  and  Exchange  Program  and  the  Aerospace 
Applications  Studies  Program.  The  results  of  AGARD  work  are  reported  to  the  member  nations  and  the  NATO 
Authorities  through  the  AGARD  series  of  publications  of  which  this  is  one. 

Participation  in  AGARD  activities  is  by  invitation  only  and  is  normally  limited  to  citizens  of  the  NATO  nations. 


The  content  of  this  publication  has  been  reproduced 
directly  from  material  supplied  by  AGARO  or  the  authors. 


K Published  April  1975 


FOREWORD 


There  has  been  a considerable  increase  in  interest  for  predicting  and  preventing  flutter  on 
figliter  aircraft  carrying  stores.  Because  of  the  importance  of  the  application  to  the  NATO 
countries’  air  forces,  the  Structures  and  Materials  Pan  .1  of  the  NATO  Advisory  Group  for 
Aerospace  Research  and  Development  (AGARD)  held  a one-half  day  Specialists’  Meeting 
on  “Wing-with-Stores  Flutter”  at  the  Structures  and  Materials  Panel’s  39th  Meeting  in  Munich 
on  9 October  1974. 

The  addition  of  the  stores  significantly  changes  the  dynamic  characteristics  and  often 
adversely  affects  the  flutter  properties  through  drastically  reduced  flutter  speeds.  Flutter 
speed  placards  are  thus  frequently  required  thereby  severely  restricting  potential  speeds  and 
mission  performances.  Moreover,  many  different  stores  can  be  carried  at  each  carriage 
station  requiring  evaluation  of  a very  high  number  of  potentially  flutter  critical  store 
combinations.  The  purpose  of  this  meeting  was  to  discuss  the  latest  state-of-the-art  in  the 
various  countries,  and  op  imum  methods  for  avoiding  restrictive  placards  and  for  rapidly 
and  economically  evaluating  the  many  store  combinations.  These  proceedings  contain  the 
papers  presenteo. 

The  AGARD  Structures  and  Materials  Panel  expresses  its  appreciation  to  the  authors 
and  to  the  Specialists’  Meeting  Chairman,  Mr  O.Sensburg  of  M.B.B.,  Germany,  for  their 
significant  contributions  to  this  important  figliter  and  ground  support  aircraft  topic. 


Walter  J.Mykytow 
Chairman,  Working  Group  on 
Aeroelasticity  and  Unsteady  Aero-Dynamics 
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CALCULATION  METHODS  FOR  THE  FLUTTER  OF  AIRCRAFT  WINGS  AND  EXTERNAL  STORES 


LI.  T.  Niblett  and  J.  C.  A.  Baldock 
Scructures  DeparCcaent 
Royal  Aircraft  Establishment 
Fsmborough,  Hampshire 
England 


SUMMARY 

This  paper  reports  theoretical  work  at  RAE  on  the  wing-uith-stores  problem  which  has  been  in  the 
fields  of  structural  representation,  the  solution  of  the  flutter  equations  and  the  prediction  of  flutter 
characteristics  from  structural  properties. 

The  subjects  covered  in  some  detail  are: 

(a)  a comparison  of  the  normal  modes  calculated  for  a wing-with-stores  from  some  of  the  normal  modes 
of  the  oare  wing  and  discrete-load  modes  with  those  calculated  ftor  the  full  flexibility  matrix; 

(b)  the  basis  of  a cocq>uter  program  which  traces  the  loci  of  constant  flutter  speeds  when  two 
structural  parameters  vary;  and 

(c)  the  interpretation  of  the  loci  of  constant  flutter  speed  in  terms  of  modal  shapes  and  frequencies 
with  the  object  of  assessing  the  ;iiost  critical  store  combinations. 


1 INTRODUCTION 

This  paper  describes  work  that  has  been  done  on  a few  of  the  techn.ques  of  the  theoretical  treat- 
ment of  the  flutter  of  wings  with  scores  favoured  at  RAE. 

Among  the  data  needed  for  wing-with-stores  flutter  calculations  is  a descriptior  of  the  structural 
properties  of  the  wing  which  is  adequate  enough  to  give  its  distortions  under  a wide  variety  of  score 
inertia  forces  yet  economic  in  computing  time  and  space.  That  such  a description  is  given  by  a combination 
of  normal  and  discrete-load  modes  has  been  shown  by  some  calculations  which  will  be  described. 

Since  this  combination  will  give  the  distortions  of  the  wing  under  practically  any  score  inertia 
forces  from  only  comparatively  few  coordinates  it  is  very  suitable  for  use  with  a computer  program 
which  automatically  follows  contours  of  constant  flutter  speed  against  two  of  the  mass  and/or  stiffness 
properties  of  the  store.  Such  contours  are  useful  in  that  they  give  the  most  critical  score  ermbinations 
and  the  basis  of  an  automatic  contour-following  program  is  given. 

The  contours  will  indicate  critical  combinations  of  store  parameters,  ^ method  is  given  for  a 
rapid  assessment  of  the  physical  basis  for  this  combination. 

2 DISCRETE-LOAD  MODES 

Comparable  calculations  have  been  made  of  the  normal  modes  of  an  encastrd  wing  with  stores.  The 
structural  system  consisted  of  the  bare  wing,  two  flaps  with  beams  connecting  them  to  the  bare  wing 
structure  and  two  stores  with  pylons  connecting  them  to  the  bare  wing  structure  (see  Fig.l).  The  object 
of  the  calculations  was  to  discover  how  close  an  approximation  to  the  'true'  modes  of  the  wing  with 
stores,  defined  as  the  results  of  calculations  using  all  the  information  available,  was  obtained  if  only 
a limited  number  of  the  normal  modes  of  the  wing-without-s tores  together  with  appropriate  'discrete-load' 
modes  of  the  wing  were  used.  The  'discrete-load'  modes  (sometimes  known  as  'junction'  modes)  are  the 
deflection  shapes  taken  up  by  the  wing  when  it  is  loaded  statically  at  the  pylon  o .achment  points  by  the 
types  of  load  the  stores  generate  dynamically  which  give  significant  distortions  of  the  wing,  i.e.  heaving 
force  and  rolling- and  pitching-moments . These  inodes  are  proportional  to  the  columns  of  the  flexibility 
matrix  for  the  wing  which  pertain  to  those  particular  types  of  loads  acting  at  the  pylon  attachment 
points . 

Thus  an  allowed  displacement  of  the  wing  can  be  written 


t‘^ab>Y^a^ac3K'<>c! 


where  is  a vector  of  wing  displacements, 

is  the  rectangular  matrix  of  wing  displacements  in  the  normal  modes  whose  generalised 
coordinates  are  qjj  , 

F^^  is  the  wing  flexibility  matrix, 

is  a non-square  permutation  matrix  which  selects  frcji  F^  the  columns  appropriate  to  the 
discrete-load  modal  coordinates  q^.  , 

ana  factoring  by  y simply  replaces  the  dimensions  of  flexibility  by  the  dimensions  of 
displacement . 
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The  inertia  matrix,  if  the  mass  matrix  associated  with  the  wing  deflection  points  is  A , is 
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ba  aa  aa  ac 
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(2) 


where  a^  is  a diagonal  matrix  of  the  normal  mode  inertias,  and  the  stiffness  matrix  is  given  by 
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(3) 


If  the  initial  data  is  in  the  form  of  a stiffness  matrix,  Bjj  , rather  than  a flexibility  matrix, 
the  discrete  load  modes,  , are  given  by  the  solution  of  an  equation  sUch  as 


E. jQ,  •*  yQj  P 

dd  dc  da  ac 


(A) 


the  rhs  giving  the  generalised  forces  in  the  stiffness  modes  due  to  the  discrete  loads.  The  total 
stiffness  matrix  is  then  given  by 


*^bd 

‘^bd 
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This  is  not  possible  if  is  singular.  One  of  the  usual  methods  of  overcoming  this  type  of 

problem  is  to  split  the  coordinates  into  coordinates,  q^  say,  which  involve  distortion  of  the 

structure  and  coordinates,  q^  say,  which  describe  rigid-body  displacements.  The  inertia  matrix 


and  the  stiffness  matrix 
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The  q^  coordinates  are  replaced  by  qj  coordinates  using  the  fans  formation 


(say) 


(7) 


where  the  qj  coordinates  are  fewer  than  the  qj  coordinates  by  the  number  of  rigid  body  freedoms 
present.  The  constraints  on  the  qj  modes  are  the  equivalent  of  the  addition  of  rigid-body  deflections 
to  the  qf  modes  such  that  the  new  inodes  are  orthogonal  to  further  rigid-body  deflections.  Eq.(4) 
transformed  and  premultiplied  by  ( is 
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which  multiplied  out  and  transformed  back  gives 
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dc 


<^rf 


e"'(Q,  + Q;  Q )P 

ff  ^fa  fr^ra  a 


(9) 


An  approximate  description  of  the  system  could  be  obtained  by  adding  some  of  cne  natural  modes,  including 
all  the  rigid-body  modes  to  these.  This  has  not  been  tested  numerically. 

Some  simple  exploratory  calculations  have  been  made  in  the  case  of  'discrete-constraint'  modes,  which 
is  another  approach  useful  when  the  stiffness  matrix  is  singular,  r.ie  term  is  used  to  denote  the  fundamen- 
tal modes  of  the  system  described  by  the  equation 


1^ 
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which  each  of  the  constraint  equations,  collectively  described  by 


“ra^d'^d 


- 0 


(10) 


(11) 


used  separately  as  a subsidiary  equation.  The  constraints  applied  are  the  equivalents  for  the  discrete 
loads  on  the  actual  system,  i.e.  a discrete-load  mode  under  a force  is  replaced  by  a discrete-constraint 
mode  with  position  constraint  at  the  point  of  action  of  the  force.  The  constraints  can  be  applied  before 
or  after  the  system  is  transformed  to  one  orthogonal  to  the  body  freedoms. 

The  'true'  modes  of  the  wing  with  flaps  and  stores  were  obtained  by  setting  up  the  dynamical  equations 
for  the  wing  primary  structure,  the  flap  beams  and  flaps  and  the  pylons  and  stores  and  solving  them 
simultaneously.  The  wing  primary  structure,  which  originally  had  98  degrees  of  freedom,  was  approximated  to 
by  its  lowest  70  normal  modes  and  6 discrete-load  modes.  Each  flap  beam  had  2 distortion  modes  and  the 
flaps  had  30  degrees  of  freedom  in  toto.  The  pylons  each  had  3 distortion  degrees  of  freedom  covering 
motions  in  pitch,  roll  and  torsion.  The  total  nundier  of  degrees  of  freedom  subsequent  to  the  application 
of  continuity  constraints  was  112.  The  approximate  systems  were  made  up  of  up  to  10  normal  modes  for  the 
bare  ;ing  and  flap,  6 discrete-load  modes  of  the  wing  . rimary  structure,  with  the  coefficients  modified 
to  take  account  of  the  flaps,  and  6 pylon  distortion  modes,  ai'ding  up  to  a maximum  of  22  modes.  The 

frequencies  of  the  first  10  modes  of  these  systems  were  compar'.d  with  those  of  the  true  system.  For  the 

biggest  approximate  system  the  maximum  error  was  found  to  be  0.044.  The  errors  were  not  all  of  the  same 
sign  which  suggested  that  the  theoretical  accuracy  of  the  approximation  was  greater  than  the  numerical 
accuracy  with  which  the  modes  were  calculated.  Were  the  differences  in  frequencies  larger  it  might  have 
been  necessary  to  formulate  other  criteria  of  accuracy.  A comparison  of  the  frequencies  obtained  when 
fewer  normal  and  discrete-load  modes-  «ere  used  showed  that  the  modes  were  dominated  by  the  effect  of  the 
heavy  stores . A good  indication  of  tl  is  was  that  the  approximations  with  6 discrete-load  modes  only  were 
better  than  those  with  10  normal  modes  only  apart  from  those  which  had  a large  flap  or  beam  content.  Also 
the  approximations  were  significantly  worse  when  the  discrete-load  modes  for  the  outer  pylon  fixing  points 
were  excluded . There  is  no  reason  to  believe  however  that  the  method  would  not  give  good  approximations 

with  lighter  stores  and  it  certainly  works  in  the  no-stores  case.  It  would  appear  that,  in  the  case  of 

the-  structure  considered  here  at  least,  accurate  results  can  be  achieved  by  including  all  the  discrete-load 
and  pylon  modes  and  normal  codes  up  to  and  includin.j  the  first  whose  frequency  is  above  the  maximum 
frequency  of  interest.  Further  details  of  the  calculations  are  given  in  Ref.l. 

The  approximate  method  used  here  has  wider  applications  in  the  representation  of  substructures.  For 
instance,  from  some  further  calculations  that  have  been  made*,  it  can  be  surmised  that  fewer  modes  of  Che 
wing  primary  structure  and  flaps  would  have  been  needed  to  calculate  the  wing-wich-flaps  normal  modes  if 
some  of  the  modes  for  the  wing  primiry  structure  had  been  the  modes  under  discrete  loads  at  the  flap  beam 
junctions  and  some  of  the  modes  for  the  flaps  had  been  discrete-constraint  or  -load  modes  for  the  beam 
junctions . 

3 CONTOURS  OF  CONSTANT  FLUTTER  SPEED 

A representation  of  the  wing  structure  such  as  the  above  is  useful  for  providing  data  for  computer 
programs  which  calculate  contours  of  constant  flutter  speed  against  two  score  mass  or  pylon  stiffness 
parameters  for  the  aerodynamic  coefficients  for  one  set  of  modes  will  cover  an  infinite  variety  of 
store/pylon  conditions.  The  contour  program  used  at  RAE  is  based  on  inverse  iteration  and  circular 
extrapolation . 

The  form  of  the  flutter  equations  generally  used  in  the  UK  can  be  written 


j^AX^  + (B  + Dv"’)x  + C + Ev"^jq  = 0 


(12) 


where  A,  D and  E are  reel  matrices  associated  with  the  structural  properties  of  the  system,  B and  C 
are  real  matrices  associated  with  the  aerodynamic  properties,  X is  the  (imaginary)  frequency  parameter, 

V the  speed  parameter,  and  q the  latent  vector. 


An  abbreviated  form  of  Eq.(l2)  is 


j^AX^  + BX  + cjq 


which  can  be  written  in  linear  form  as 
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(14) 


where  p = Xq . 


Iteration  according  to  the  equation 
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results  to  convergence  to  the  latent  root  nearest  n . 
Eq.(I5)  can  be  replaced  by  two  equations 


“ Ps+1 


(15) 


(16) 


and 

(rA  + 5)0^^,  + Cq^^,  - -(X  - . (17) 

Substituting  for  q^^l  in  Eq.(17)  from  Eq.(l6)  and  multiplying  by  n gives 

(v“A  + vB  + C)Pg+,  “ (X  - i')(Cqg  - rAq^)  (18) 

Once  one  point  on  the  contour  has  been  found  the  stability  at  any  point  near  to  it  can  be  found  by 
making  the  appropriate  modifications  to  the  coefficient  matrices  and  iterating  using  Eq.(I8)  and  Eq.(l6) 
with  the  frequency  parameter  at  the  known  point  on  the  contour  as  the  pole  n . The  sign  of  the  real  part 
of  (X  - it)  determines  which  side  of  the  contour  the  trial  point  lies. 

The  choice  of  where  to  search  for  the  contour  is  based  on  the  assumption  that  an  arc  of  a circle 

will  be  a good  approximation  over  a limited  range.  In  general  three  previous  points  on  the  contour  will 
be  available  (say  A,  B and  C of  Fig. 2) . A circle  is  drawn  through  these  points  and  the  positions  of 
the  last  two  points  are  defined  in  terms  of  the  angles  (a,B)  between  the  radii  through  them  and  the  first 
point  and  the  next  point  is  searched  for,  using  Muller  interpolation,  along  a radius  which  is  at  the  same 
angle  to  the  radius  through  the  last  point  as  that  is  to  the  radius  through  the  point  before  (y  “ B - a) . 

The  search  is  commenced  on  the  circle  (at  Dl)  and  continues  until  (X  - a)  has  a small  enough  real  part  or 
until  it  is  predicted  that  the  contour  lies  more  than  a limited  distance  from  Dl.  If  the  latter  is  the 
case  the  value  of  y is  halved  and  this  can  happen  a number  of  times  if  the  contour  proves  difficult  to 

find.  On  the  other  hand  if  the  contour  proves  easy  to  find  the  value  of  y used  for  the  next  point  may 

be  made  greater  subject  to  maximum  limits,  one  of  which  is  part  of  the  data  and  the  other  rt  of  the 
program.  Thus  the  intervals  between  points  on  the- contours  vary  constantly  with  the  tortuosity  of  the 
contour.  One  of  the  main  objects  in  writing  the  programs  has  been  to  reduce  the  awunt  of  human  interven- 
tion to  the  minimum.  Further  details  of  the  program  are  given  in  Refs. 2 and  3. 

The  method  does  not  cater  for  the  exact  lining-up  of  assumed  and  calculated  frequency  parameter  but 
if  it  is  so  wished  the  aerodynamic  coefficient  matrices  can  be  made  consistent  with  n . 

4 INTERPRETATION  OF  CONTOURS  OF  CONSTANT  FLUTTER  SPEED 

Fig. 3 shows  the  diagram  from  a typical  output  of  the  contour  program.  The  store  mass  parameters  are 
Che  mass  and  the  square  of  the  radius  of  gyration  of  a store  on  an  inboard  pylon.  This  output  was 
preceded  by  conventional  solutions  for  the  system  roots  versus  airspeed  at  the  mass  parameter  points  shown 
on  Fig. 3,  in  order  to  find  critical  flutter  speeds,  which  vi're  used  as  starting  values  for  the  loci  of 
flutter  speeds  versus  one  mass  parameter  with  the  other  mas.--  parameter  held  constant.  These  solutions 
were  obtained  by  one  of  the  options  built  in  to  the  contour  program,  and  they  themselves  are  used  for 
starting  values  for  finding  the  contours. 

It  may  be  seen  from  Fig. 3 that  there  is  a sharply  critical  combination  of  store  mass  and  radius  of 
gyration.  The  physical  basis  for  this  critical  combination  is  of  obvious  interest,  particularly  for 
assessing  the  likely  change  in  the  critical  mass  parameters  for  a change  in  another  parameter,  such  as 
pylon  stiffness.  Diagnostic  procedures  are  available'*  for  interpreting  flutter  conditions,  but  the 
number  of  parameters  in  a wirg-with-s tores  calculation  is  so  large,  and  some  of  the  stiffness  data  so 
uncertain  in  the  early  stages  ''f  design,  that  a very  rapid  and  easily  obtained  indicator  is  desirable. 

The  frequencies  at  zero  airspeed  of  the  modes  involved  in  the  flutter  solutions  are  easily  found,  and 
variations  of  tnese  with  one  mass  parameter  can  be  helpful.  Fig. 4 shows  the  zero  airspeed  frequencies 
verjus  store  mass  for  a value  of  store  radius  of  gyration  close  to  that  of  the  most  critical  store  combina- 
tion. Comparing  Fig. 3 and  Fig. 4,  it  may  be  seen  that  the  most  critical  store  mass  combination  is  associated 
with  the  proximity  of  the  second  and  third  normal  mode  frequencies.  Alternatively,  it  may  be  possible  to 
divide  the  system  of  equations  into  a set  of  predominantly  'bending'  modes  (e.g.  wing  vertical  bending, 
store  pylon  lateral  bending)  and  a set  of  predominantly  'torsion'  modes  (e.g.  wing  torsion,  pylon  pitching) 
and  the  variation  of  the  zero  airspeed  frequencies  from  these  sets  separately  a ly  be  plotted  against  a 
store  mass  parameter.  Fig. 5 is  an  example  of  this  presentation  for  the  same  system  as  in  Figs. 3 and  4. 

It  may  be  seen  that  on  Fig. 5,  the  critical  store  mass  co.nbination  is  associated  ’’th  the  proximity  of  the 
frequencies  of  the  second  'bending'  mode  and  the  first  'torsion'  mode. 

Obviously,  a serious  study  of  the  flutter  coupling  should  involve  the  shapes  of  the  modes  of  the  system, 
but  the  large  nunber  of  specific  configurations  of  potential  interest  in  a wing-with-s tores  analysis  can 
easily  lead  to  bulk  of  mode  shape  results  that  is  difficult  to  assimilate.  Modal  frequency  results  are 
more  compact  and  can  provide  a useful  quantitative  guide  to  the  effects  of  parameters,  so  that  the  plotting 
of  mode  shapes  can  be  confined  to  a bare  minimum  of  the  most  important  specific  configuration  cases. 
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The  paper  describes  the  flutter  clearance  of  U.K.  Jaguar  using  a 
combination  of  mathematical  modelling,  grc^nd  resonance  end  flight 
testing  leading  ultimately  to  clearance  of  a wide  range  of  under 
wing  stores.  Some  improvements  in  modal  modelling  techniques  are 
outlined.  These  should  enable  reductions  in  future  ground  and 
flight  testing  times  to  be  made. 


1.  unsopncnoH 


Kodern  vilitaiy  aircraft  are  aapected  to  carry  a wide  range  of  external  stores,  cany  of 
these  on  flexible  underwing  pylons.  For  a relatirely  saall  aircraft  such  as  Jaguar  where 
external  stores  are  carried  on  two  undarwing  pylons  and  an  underfUselage  pylon  (Figure  1) 
and  the  masses  of  the  heavier  stores  are  of  the  ease  order  as  the  wing  itself,  it  is 
evident  that  the  structural  modes  and  hence  the  flutter  status  of  the  aircraft  will 
change  considerably  free  one  store  combination  to  another.  On  Jaguar,  further  complications 
arise  because  multiple  stores  can  be  carried  on  a "tandem  beam"  and  side  by  side  carriers, 
which  are  fitted  to  the  underwing  pylons,  and  these  introduce  additional  flexibility. 

The  resulting  problem  of  clearing  for  flutter  a wide  range  of  "different"  etircraft  has 
been  epproached  by  a ccabination  of  theoretical  calciUationa,  ground  and  fli^t  tests. 

This  report  describes  the  procedures  followed  concentrating  on  the  techniques  employed 
and  discusses  soese  of  the  results  obtained. 


2.  CLEAHAMCE  ESOCSDOHE 

It  was  decided  that  the  moat  efficient  and  economical  method  of  demonstrating  the  clearance 
was  by  utilising  as  extensively  as  possible  a theoretical  modal  of  the  aircraft  structural 
modes  and  their  associated  aerodynamics.  The  calculated  modes  could  be  checked  against 
the  iqiproprlate  ground  resonance  test  data  and  if  necessary  modified  to  achieve  a successful 
natch.  With  the  structural  representation  thus  validated,  flutter  behaviour  could  be 
rechecked  by  calculation  and  the  critical  configurations  identified  for  investigation  in 
flight  tests.  Hopefully,  these  flight  tests  would  : 

(a)  achieve  a fli^t  flutter  dearsnee  for  these  critical  configurations 

(b)  provide  reliable  measured  data  on  critical  structural  mode  behaviour 
to  corroborate  predictions  and  hence, 

(c)  justify  an  analogous  clearance  by  calculation  for  all  other  nen- 
critical  stores  without  the  need  for  further  specific  instrumented 
flight  tests. 

This  approach  aims  to  derive  maximum  benefit  from  extensive  but  relatively  cheap  calcu- 
lations in  order  to  confine  flight  testing  by  a few  key  configurations.  It  hinges  crucially 
upon  : 


(a)  a confident  choice  of  the  key  configurations  for  flight  tests, 
and 

(b)  Measuring  modal  data,  in  flight,  which  is  extensive  and  reliable 
enough  to  allow  meaningful  comparison  with  predictions  in  order 
to  support  the  analogous  clearance. 

Its  final  success  depends  upon  a good  correlation  from  (b) 


3.  THEORETICAL  MODELLING 


Thus  the  underwing  stores  flutter  clearance  program  relied  fundamentally  on  the  quality 
and  efficiency  of  the  theoretical  calculations.  It  was  recognised  that  these  must  not 
only  predict  potentially  critical  store  combinations  with  sufficient  accuracy  but  must  also: 

(a)  allow  the  key  structural  and  aerodynamic  parameters  to  be  Identified 

(b)  allow  flutter  status  of  any  new  stores  to  be  established  quickly 
A theoretical  model  has  been  set  up  using  a combination  of  branch  modes  for  : 

(i)  clesui  aircraft,  and 
(ii)  pylon/storos. 
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This  cosfponest  aet  baa  eereral  adrantaseo  oTer  th*  normal  modea  of  the  aircraft  including 
underving  atores  : 

(a)  the  clean  aircraft  aodea  and  flutter  behaviour  ia  veil  underetood 
and  validated  froa  clean  aircraft  ground  resonance  tests  and 
flight  flutter  clearance  data. 

(b)  Sie  pjlon  Bodes  have  been  chosen  so  that  oodel  stiffnesses  can  be 
iuertad  easily  i^oa  static  stiffness  test  data.  Ihese  in  turn  can 
be  verified  by  calculating  the  associated  normal  modes  ard  comparing 
these  with  rig  oounted  pylon-stores  resonance  test  data. 

(c)  Ihe  modal  generalised  inertias  are  computed  once  for  all  for  each 
store  on  each  pylon  and  combinations  constructed  by  appropriate 
additions. 

(d)  As  Bany  pylon  modes  as  required  can  be  introduced,  to  represent 
E.B.n.  effects,  tandem  beass,  side  by  side  carriers  etc. 

(e)  The  nodel  generalised  aerodynamics  ne^ect  contributions  firom 
pylon  and  stores,  and  hence  can  be  computed  once  for  all  and  are 
Independent  of  store  configuration.  So  far  as  can  be  Judged  from 
calculation  and  test  data,  this  simplification  introduces  errors 
within  about  5^  on  flutter  speeds. 

Cieoretically  derived  modes  were  used  for  the  basic  aircraft  representation,  lliese  were 
computed  from  wing,  tail  and  fuselage  branch  nodes  and  in  the  initial,  pre  G.B.T.  stages 
were  confined  to  fixed  root  wing  mo^s.  Although  this  sacrifices  some  accuracy  in  flutter 
speed  prediction,  the  basic  store  flutter  couplings  are  unchanged  and  the  inich  reduced 
number  of  modes  allowed  a significant  economy  in  computation  and  eased  interpretation  of 
the  calculation  output,  dhis  allows  a better  understanding  of  the  fundamental  couplings 
and  sensitivities.  Finally  the  full  aircraft  free-free  theoretical  modes  were  incorporated 
into  the  model  and  have  been  used  to  match  aircraft  G.R.T. , to  define  the  fli^t  flutter 
test  configurations,  and  provide  comparison  with  flight  data. 


4.  FLOma  CALCPLATIONS 


It  was  clear  at  the  outset  that  fuel  tank  sequencing  would  need  particular  attention. 

Two  possible  flutter  couplings  were  identified  with  the  prototype  fuel  sequencing  system, 
using  fixed  root  wing  modes  (Figure  3)«  and  these  survive  when  outboard  stoi'es  are  intro- 
duced. (Figure  4).  Corresponding  calculations  using  complete  aircraft  modes  showed  that 
althou^  these  wing-store  flutter  speeds  changed,  the  flutters  still  occurred  at  similar 
frequency,  with  the  same  store  mass  properties  but  were  milder  in  type.  (Figure  3)* 

Detailed  inspection  of  the  flutter  mode  shapes  revealed  that  these  wing  flutters  arise 
from  couplings  between  modes  which  resemble  fundamental  wing  bending  and  torsion. 


Further,  it  can  be  shown  that  : 

(a)  componenv  modes  of  wing  bending,  inboard  and  outboard  pylon  lateral 
bending  combine  to  produce  three  "wing  bending"  modes 

(b)  component  modes  of  wing  torsion,  inboard  and  outboard  pylon  pitch 
combine  to  produce  three  "wing  torsion"  modes. 

The  wing  store  flutters  which  have  been  found  for  a wide  range  of  stores  all  arise  from 
couplings  between  modes  of  these  types.  These  effects  can  be  illustrated  most  clearly 
for  the  case  with  inboard  tank  alone  (Figure  6)  when  the  two  types  of  flutter  are  both 
associated  with  the  first  "wing  torsion"  mode  frequency  coinciding  with  first  and  second 
"wing  bending"  modes. 

Similar  relationships  can  be  identified  in  the  more  complex  inboard  tank/outboard  bomb 
case  (Figure  ?)•  It  is  clear  that  pylon  mode  stiffnesses  are  the  key  structural  feature 
in  wing-store  flutter.  Hie  stores  which  are  predicted  to  be  most  critical  will  change  if 
pylon  stiffnesses  change,  and  some  further  evidence  of  this  is  given  in  Figure  8.  Note 
that  these  trends  are  characteristic  of  frequency  coincidence. 
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5.  OaODKD  HES(KUMCS  TESmi^ 


The  raliance  on  theorotical  calculations  and  their  sensitiTity  to  pjlon  stiffnesses 
dictated  rig  Mounted  stiffness  and  resonance  tests  on  inboard  and  outboard  pylons  for 
a representatiTe  range  of  store  Bass  properties.  Subsequently,  a few  conplete  aircraft 
resonance  tests  w-'re  collated  with  a typical  store, 

(a)  on  inboard  pylon  only 

(b)  on  outboard  pylon  only 

(c)  in  coabination 

to  Teritj  the  orerall  Bathesatical  model  and  in  particular  the  important  pylon  stiffnesses. 

All  the  normal  modes  were  measured,  up  to  and  including  the  highest  frequency  pylon  induced 
wing  bending/torsion  mode  which  might  be  a source  of  flutter  coupling,  based  on  the  theoretical 
review. 

The  measured  normal  modes  ware  compared  with  the  ^propriate  predictions  (Figure  9),  where 
labels  have  been  attached  which  describe  the  predominant  structural  motion  in  that  mode. 

Ihe  correspondence  is  good  when  only  the  inboard  or  only  the  outboard  store  is  fitted,  but 
relatively  poor  when  these  are  in  combination.  An  alternative  set  of  pylon  stiffness  values 
was  derived  (Calc. 2}  which  provide  a bettor  natch  for  all  three  configurations.  These 
adjustments  were  small  for  the  inboard  pylon  but  up  to  O.^x  the  values  established  by  rig 
tests  for  the  outboard  pylon. 

It  was  conjectured  that  the  6 wing  normal  modes  which  were  retained  in  the  clean  aircraft 
representation  might  be  insufficient  to  define  local  wing  distortion  near  to  the  pylon- 
wing attachment.  Failure  to  account  for  these  effects  would  be  consistent  with  compensation 
via  a reduction  in  pylon  stiffnesses,  as  evidenced  by  the  aircraft  ground  resonance  test 
matching.  Farther  evidence  for  this  is  presented  in  Figure  10.  This  demonstrates  the 
accuracy  of  calculating  tdng-store  modes  (for  a rigid  pylon)  from  a finite  number  of  clean 
Jaguar  wing  modes.  IVom  this  there  are  substantial  inaccuracies  at  low  frequencies  even 
with  twelve  clean  wing  modes.  These  can  however,  be  avoided  even  with  only  3 clean  wing 
modes,  provided  that  additional  wing  modes  are  introduced  which  are  the  statically  distorted 
wing  shapes  for  unit  vertical  shear,  bending  moment  and  pitching  moment  applied  at  each 
pylon-wing  attachment.  .Such  discrete  load  modes  clearly  retain  significant  local  distortion 
effects  which  are  discarded  in  the  hi(^er  order  clean  wing  modes,  and  they  retain  the  concept 
and  economy  of  the  component  mode  approach. 

For  subsequent  applications,  it  might  be  expected  that  the  increased  acoiracy  of  predictions 
will  allow  the  critical  store  configuration  to  be  confidently  defined  prior  to  ground 
resonance  testing,  given  some  test  data  on  pylons.  In  addition,  it  mi^t  be  considered 
sufficient  in  future  only  to  resonance  test  one  configuration,  eg.  with  both  pylons  loaned. 


6.  SELECnON  OF  HOST  CRITICAL  STORES 


The  basis  of  selection  on  Jaguar  is  juxtaposition  of  wing  bending  and  torsion  type  mode 
frequencies,  supported  by  representative  flutter  calculations  for  each  possible  region 
of  coincidence. 

Inboard  stores  only 

Figure  11  shows  the  frequency  of  the  normal  modes  of  bending  and  torsion  types  which  fall 
within  the  important  frequency  range  for  a wide  range  of  projected  stores.  It  is  evident 
that  for  all  these  stores,  the  bending  and  torsion  mode  frequencies  are  substantially 
independent  and  uniquely  related  to  pylon  lateral  bending  and  pitch  mode  generalised 
Inertias,  respectively;  that  is,  to  store  mass  and  to  store  pitching  inertia.  Super- 
impositon  of  Figures  11(a),  11(b)  reveals  two  potential  flutter  couplings: 

(1)  Bi/Tl  (10  Hz)  - for  heavy,  slender  stores,  eg.  full  fuel  tank 

(2)  B2/T1  (10  Hz)  - for  medium  wei^t,  short  stores,  eg.  bombs/part  full 

fuel  tank. 

Flutter  calculations  are  already  available  in  these  classifications;  the  mild  flutter  on 
Figure  3 is  B1/T1  and  the  severe  flutter  is  B2/Tl. 
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Updated  flutter  calculations  with  new  stiffnesses  confirm  that  the  latter  is  more 

important , and  reveal  that  the  fUel  tank  with  protot^e  fuel  sequencing  isjetem  is  still 
the  most  critical  store  for  this  t^rpe  of  flutter.  Similarly,  the  fuel  tank  is  most 
critical  for  the  altemat  we,  milder  type. 


Outboard  Stores  only 

Similar  nodal  frequency  characteristics  can  be  identified  for  outboard  stores  only  (Figure  12). 
Superimposition  is  this  case  reveals  : 

(1)  B1/D1  coupling  - medium  wei^t,  slender  stores,  outside  the  present 

range. 

(2)  B2/E1  coupling  - for  most  of  the  stores  in  the  current  inventory 

Thus,  flutter  calculations  were  necessary  for  nearly  all  the  projected  outboard  single  stores. 
These  revealed  that  only  the  B1/Tl  coupling  is  potentially  critical.  Very  mild  couplings  and 
no  flutter  resulted  from  the  B2y^  frequency  coincidence.  Thus,  none  of  the  stores  in  the 
current  inventory  are  critical  tdien  carried  alone  on  outboard  pylon. 

Inboard  and  Outboard  Stores 


The  problem  becomes  more  cocsplicated  now,  as  there  are  three  bending  modes  and  two  torsion 
modes  to  consider.  Figure  13  shows  the  possible  bending  mode  frequencies  and  Figure  14  the 
possible  torsion  made  frequencies  as  a function  of  inboard  and  outboard  store  inertia 
properties.  Observe  that  there  are  six  possible  sources  of  coupling  : 

(a)  B1/T1  flutter  is  only  possible  when  heavy  slender  stores  are  carried  on  the  Inboard 
pylon.  The  most  critical  combination  is  with  a small  light  outboard  store,  (eg.  full 
fuel  tank  in  Figure  6). 

(b)  B1/T2  flutter  would  require  store  mass  and  inertia  outside  the  range  considered  at 
present. 

(c)  B2/T1  is  most  likely  for  store  combixiationB  with  moderately  heavy  stores  on  the  inboard 
pylon  and  li^t  stores  on  the  0Uv«,  ,rd  pylon,  (the  505^  full  fuel  tank  case  in  Figure  6). 

(d)  B2/T2  flutter  would  require  an  unrealistic  store  mass  and  inertia  combination. 

(a)  B3/tl  is  most  critical  for  a combination  of  inboard  and  outboard  stores  which  are 
similar  in  weight  (eg.  the  30%  fuel  tank  plus  outboard  bomb  flutter  in  Figure  7)> 

(f)  B3/T2  is  most  critical  for  heavy  outboard  stores.  (The  full  fuel  tank  plus  outboard 
bomb  flutter  in  Figiire  ?)♦ 

The  most  critical  store  combinations  were  established  for  flutter  calculations  in  each 
of  these  regions.  From  these  the  potentially  critical  flutters  occur  when  the  fuel  tank 
is  carried  on  the  inboard  pylon  both  with  and  without  outboard  pylon  stores.  T^es  B1A1 
and  B2/T1  havo  higher  flutter  speeds  than  the  inboard  only  stores  case,  and  B3/T1  and  B3/TZ 
have  the  lowest  flutter  speeds  of  all. 


SDHHARI 


From  this  review,  the  following  projected  store  configurations  emerge  as  critical,  in  order 
of  concern  : 

(1)  tank  Inboard  plus  bomb  outboard 

(2)  tank  inboard  only 

(3)  tandem  beam  stores  inboard  plus  outboard  bomb 

(4)  tandem  beam  stores  inboard  only 
and  these  were  selected  for  flight  flutter  testing. 

The  multiple  carried  stores  trere  predicted  to  be  less  critical  for  flutter  than  the  fuel 
tank,  but  there  was  a possibility  of  significant  structural  response  in  several  poorly 
damped,  local  beam  modes. 
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7.  njGgT  TEssma 


Strattgr: 

Calcttlatioss  hare  ehr«n  the  inboard  fuel  tank  wcs  the  ooat  critical  etore  for  flutter. 

It  vas  decided  to  flutter  teat  thin  alone  initially,  anticipacing  that  the  experience  and 
data  gained  would  help  considerably  in  interpreting  the  sore  cooplicated  and  most  critical 
configuration,  according  to  prediction,  with  outboard  bomb. 

Ciis  approach  was  also  adopted  with  tandem  beams  inboard.  In  each  case  (4  configurations) 
specific  excitation  and  instrumentation  vas  fitted. 

Finally,  on  successfully  completing  the  above  programme,  single  flights  wore  planned  with 
less  critical  stores,  which  were  confidently  expected  to  be  flutter  free.  In  these  flights, 
no  specific  excitation  uas  attempted,  other  than  pilot  opinion  of  vibration  response  to 
turbulence,  manoeuvring  etc.,  en  route  to  the  design  clearance  speed. 

Sxcitation  was  effected  via  electrical  Inputs  to  the  main  taileron  actuators  (inherited 
from  the  clean  aircraft  clearance ) covering  an  appropriate  frequency  range.  This  proved 
to  be  more  effective  than  bonking  the  wing  or  providing  electrical  inputs  via  wing  spoilers. 

Acceleration  transducers  ^ere  mounted  on  the  stores  and  on  fuselage,  wing  and  tail. 

Signals  from  these  could  be  recorded  on  magnetic  tape  during  the  fli^t. 


Analysis; 

Baw  time  histories  of  the  response  signals  allowed  a watching  brief  to  be  kept  on  real  events. 
Computer  analyeie  of  the  signals  yielded  power  spectra  which  assisted  in  the  subsequent 
interpretation  of  these  events.  Final  quantitative  analysis  provided  modal  data  for  justifying 
fll^t  envelope  expansion  and  for  correlation  with  predictions. 


BFSOLTS 


Seme  results  are  included  here  from  flight  teets  on  the  most  critical  inboard  taidc/outboard 
boob  combination. 

A typical  spectral  variation  is  shown  in  Figure  15»  For  this,  a large  number  of  aamples, 
each  of  short  duration  were  taken.  Oieae  particular  results  are  taken  from  the  tank  noae 
vertical  transducer  and  demonstrate  in  a striking  way  how  the  two  wing  torsion  mode  frequencies 
develop  as  fuel  is  used  from  the  tank.  Because  of  short  sampling  times,  these  spectra  cannot 
be  used  to  obtain  accurate  dampings  but  taken  together  they  give  a very  clear  picture  of 
frequency  trenda  eg.  Figure  16.  Note  that  all  the  potential  flutter  couplings  are  excited 
and  that  the  critical  condition  will  bo  close  to  tank  fuel. 

More  accurate  estimates  of  frequency  and  damping  were  obtained  by  further  analysis  of  the 
impulse  response,  derived  by  spectral  division  and  inverse  fourier  transformation  using 
spectra  taken  from  longer  data  samples. 

IVequencies  and  dampings  obtained  in  this  way  are  shown  in  Figure  17  for  the  potential 
flutter  mode  of  the  tank/outhoard  bomb  store  combination.  There  is  some  scatter  in  the 
results  hut  the  overall  trend  in  frequency  and  damping  is  clear.  Also  included  in  the 
plot  are  the  predicted  trends  for  this  mode.  If  critical  damping  increase  is  introduced 
from  the  Btruettu-e  (where  about  could  be  expected  from  GSIT  results)  then  a good  match 
follows.  Observe  that  60^  fuel  state  is  indeed  the  most  critical  condition. 

In  general,  the  comparison  with  predictions  was  acceptable  for  the  form  configurations 
tested,  and  no  flutter  limitations  were  encountered. 


8.  CONCLnSIONS 


At  this  stags  : 

(a)  a theorstical  oodel  has  been  established,  and 

(b)  validated  by  ground  resonance  teat  an  key  configurations. 

(c)  the  most  critical  store  configurations  have  been  successfully  cleared  by  flight  test 
to  the  full  flight  envelope  of  H = 0.95/Sea  Level 

(d)  good  correlation  between  flight  tests  and  predictions  has  been  achieved  over  a wide 
range  of  store  configurations. 


Ve  therefore  claim  that  clearance  can  be  justified,  by  analogy,  for  all  other  projected 
stores,  without  a need  for  specific  flight  flutter  tests.  In  the  course  of  subsequent 
flights  to  assess  handling/performance  of  non-critical  stores,  there  have  been  no  indications 
of  any  flutter  problems. 

Thus,  from  flutter  tests  on  4 configurations,  a full  flight  envelope  clearance  has  been 
achieved  for  85  possible  combinations  of  known  stores. 
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FLOTEaSOT  D»  AXLES  EPPSES  DE  JOTEURS  EH  HACELIE 
par  R.  Be3t>iyi!4er 

Office  Hatiocal  d'Etudes  et  de  Rechorches  Adrospatiales  (oiSRA) 
92320  ChfitiUoT.  (France) 


Ri^Eiirog 

Le  colcul  et  lee  cesures  dee  forces  aerodynaaiquce  inetationnaires  effectuee  en  eooulencnt  eubsonique  sur 
une  naquette  dquipce  d'un  reacteur  en  nacelle  contrent  quo  I'interaction  entre  3 reaoteur  et  I'aile  reste 
negligeable.  Par  contre  lee  forces  aiirodynaniquee  induites  sur  le  reacteur  lui-cfine  par  son  propre  souveixnt 
sont  inportantes  et  jouent  un  rdle  non  negli^able  dans  le  bilen  des  forces  gcneralisees. 

Le  rdacteur  en  nacelle  a done  pu  fitre  traite  sdpar&ent.  II  a dte  assiaild  a im  cylindre  a paroi  mince 
avec  ecouleaent  interne  et  eztemc. 

Une  application  k un  cas  de  flottemsnt  a etc  faite  nontrant  1' importance  des  forces  sur  le  rdacteur  ; un 
bon  recoupement  a dtd  dtabli  entre  theoric  et  expdrience  a Each  = 0,80. 


FLUTTER  OP  HIHGS  EQUIPPED  HIIH  LARGE  ElIGIHES  TO  POD 

Calculations  and  measurements  of  unsteady  aerodynamic  forces  performed  in  subsonic  flovi  on  a model  equipped 
with  an  engine  in  pod  showed  that  the  interference  between  engine  and  wing  remains  negligible. 

It  was  also  shown  that  the  aerodynamic  forces  induced  on  the  engine  itself  by  its  o'm  oscillation  are 
important  and  give  a significant  contribution  to  the  generalized  forces. 

Account  was  taken  of  these  two  remarks  and  the  aerodynamic  forces  were  calculated  separately  on  the  engine 
which  was  assimilated  to  a thin  walled  cylinder  with  internal  and  external  flow  in  the  axial  airection. 

An  application  to  a flutter  case  shows  the  importance  of  the  contribution  of  the  forces  on  the  engine.  A 
good  agreement  was  obtained  between  theory  and  experiment  at  Ikich  number  K = C.CC. 


I.-  HtTROEPCTIOlI 

A la  suite  de  uesaccords  observes  dans  I'^volution  des  purametres  de  frequences  et  u'aaortisscmc.  ts  theo- 
riques  et  experinentaux,  sur  une  maquette  adro^lastique  de  voiliirc  equipee  d*un  reaoteur  en  nacelle,  il  e paru 
interessant  d'dtudier  I'influence  de  ce  dernier  sur  les  forces  aerodynamiques  instationnairos, 

Les  resultats  ont  nontrd  quo  les  pressions  instationnaires  induites  sur  I'ailo  pnr  un  nouvenent  du  rdacteur 
restaient  negligeables,  bien  que  mesurables. 

Par  ailleurs  les  forces  adrodynamiques  instationnOiTes  induites  sur  le  rdactexur  lui-mfine  par  oa  propre 
oscillation  se  sont  revdldes  trds  inportantes.  C'est  ainsi  quo  poui  une  oscillation  Je  1'enseu.blc  ailc  rdacteur 
autour  d’un  axe  passant  au  voisinage  des  mi-cordes,  le  moment  aea  forces  adrodyna-doues  ajissant  sur  la  nacelle, 
par  rapport  a I'axe  d' oscillation  est  dgal  au  quart  du  moment  total. 

On  a pu  en  conclure  que  les  calculs  d'aerodlasticite  ne  pourraient  pas  ndgliger  l'in."'ucnce  uu  rdacteur, 
mats  pourraiEnt  Stre  effectuds  sans  tenir  compte  du  oouplagc  avec  I'ailo,  e'eot-a-dire  que  les  forces  adroi.V>TTa- 
miques  agissant  sur  le  rdacteur  pourraient  6tre  dvaludcs  sdpardment. 

A la  suite  de  cette  remerque  des  calculs  dans  lesquels  le  rdacteur  en  nacelle  est  asoinild  L un  ou  ueva 
cylindres  circulaires  ooncentriques  a paroi  mince,  ont  dtd  ddveloppds  d'une  part  par  une  ndtnode  uo  doublet 
(calcul  developpe  par  le  H.L.R.)  d'autro  part  par  une  mdthode  de  surface  portanto  (oalc;d  dlabord  prsr  l'c;ri.?a'i) . 

Les  rdsultats  foumis  par  les  thdories  se  recoupent  tres  bien  entre  eax  et  correctenent  avec  I'ox:  Jrionoo 
et  traduisent  bien,  entre  autres,  I'effet  du  rdasteur  sur  les  Unites  de  flottement. 

L'etude  dtait  d'autant  plus  juotifide  quo  cet  effet  s'est  rc-vdld  ddf avorablo . 


II.-  DETEKJHATIOK  EXPERPSHTALE  DSS  FORCES  AEROD-n:/'’.i;;USS  inav.i’Tomi.'id:: 


La  maquetto,  qui  comprend  une  demi-voiliu'e  rigide,  munie  ou  non  de  son  rdact?ur  en  pod,  dtai;  fixe  sur  un 
support  qui  permettait  de  la  faire  osciller  sinusoldalement , 

L'incidenoe  noyenne  pouvait  fitre  rdglde,  ce  qui  a pormis  de  rdalioor  deux  sdrics  d'essais  : 

- une  sdrie  dans  laquello  I'dcoulemont  stationnaire  dtait  supercritique  (a  !;  = 0,88)  reference  |_1_, 

- une  sdrie  dans  laquelle  1' incidence  :.iOyonne  (-  45  cn)  et  lo  nombre  de  Each  maxiimim  (id  ^ 0,S)  avaient  dtd 
choisis  de  telle  faqon  que  I'dcoulemont  reste  subcritique  sur  tout  le  profil, 

Cet  article  conceme  uniqueinent  la  deuxibme  sdrie  d'essais. 

Le  rdacteur  en  nacelle  dtait  lid  rigidemont  a I'aile  et  solidaire  de  oello-ci  dans  son  nouvement  d'oscilla- 
tion.  L'axe  d'osoillation  dtait  situd  a 52  ^ des  cordes. 

Des  capteurs  de  pression  rdpartis  sui  I'aile  o'-  sur  la  nacelle  ont  pernio  de  detci-iner  la  rdpaitition  ^,co 
coefficients  de  pression  instationnaire  et  a'obtcnir,  par  integration,  Its  coefficients  sans  dimsnsxons  qui 
expriment  les  forces  ot  les  moments  pour  le  mouvement  d'osoillation  oonsiddrd. 

La  oomparaison  des  rdsultats  obtonus  avec  et  sans  rdacteur  a porr.io  de  od>arer  I'influence  uc  ce  demioi'. 


11,1  - Descri'rtion  du  pontage 

Ii'aile  dont  la  foito  en  plan  eat  donnde  par  lea  figuroa  1 et  2 dtait  equipdo  de  capteura  do  preaaion  h 
court  tenpa  de  reponse  rdpartis  sur  cinq  coidea  diatribudca  en  envergure  aroivant  une  loi  de  I-Julthopp, 

Lea  onplitudea  angulairej  du  nouvenent  de  tangagc  dtaient  rd^abloa  entre  1C  et  45  minutea.  L'ailc  ctait 
conaideree  conse  rigide,  lea  prenierea  frdquencea  proprea  se  situant  a 66  Hz  en  flexion  et  185  Hz  on  toraion, 

Le  rappel  dlaat'ique  dloit  foumi  par  une  lace  de  rigiditd  rcglable  par  valcura  diacretes,  nontde  aur  un 
arbre  portd  par  dea  laces  en  croix  dana  le  prolongeacnt  de  I'axe  de  tangage,  L'excitation  dtait  rdalisde  a I'aide 
de  4 pota  dlectrodynasicuea,  portda  par  le  eSne  arbre  et  ddlivrant  un  noaent  naximun  de  20  Hn.  Afin  de  rdduire  Ics 
mouveaents  paraaitea  de  flexion,  I'axe  de  rotation  dtait  inpoad,  a I'extrdnitd  de  l'ailc,  par  dea  cdblea  d'acier 
tendus  i travera  la  veine.  Lea  frdquencea  d'eaaai,  a divers  noabro  de  Kach,  dtaient  de  35  a 48  IIz  dana  lea  deux 
configurations  (aveo  ou  aons  rdacteur),  Lea  ceaurea  ont  dtd  realisdea  dgalecent  en  stationnaire  (frdquence  nullo). 

Lora  dea  casaia  dynaniquea  on  ceaurait  en  outre  le  conent  adrodvna’aique  instationnaire  global  autour  de 
I'exe  d' oscillation. 

L'dquipecent  de  noaure  dea  pressiona  instationnaiioa  pemettait  I'analyse,  en  tenpa  rdel,  dea  signaux 
delivrda  aicultandsent  par  dlx  capteura  do  preaaion.  Aorda  analyse  ces  signaux  dtaient  coyennds  pendant  dix 
seoondea  pour  rdduire  lea  off eta  parasites  dua  a la  turbulence  de  la  soufflerie. 

La  chalne  do  nesure  entierenarit  autonatique,  denandait  coins  d'une  minute  pour  1' acquisition  et  le  traite- 
cent  de  I'ensenble  dea  rdaultats  (70  paranetres  coaplexes). 


11,2  - Bffei:  du  reaeteur  en  nacelle  aur  la  voilurc 

Lea  pressiona  instationnaires  sont  nrdsentdos^ soua  la  force  de  coefficients  sans  dinension,  Cp.  ddfinis  : 

cTp=.Ap  C p + i Cp 

A 10  ^ 

avec  Ap=  affixe  de  la  preaaion  locale  , q = preaaion  dynaaique 

0 = affixe  de  I'anglc  d'oscillation,  pris  parallelecent  a I'dcoulccent. 

La  frdquence  rdduite  eat  ddfinie  parO^=.i^(  Sdtant  la  deni-envorgure  do  I'ailo).  A la  frdquence  nulle 

(stationnaire)  la  distribution  de  preaaion  cxpcrinontale  a dtd  idduito  dea  c...icps  de  preaaion  statique  nesurds 
pour  lea  incidencea  de  zdro  et  coins  45  ninutes. 

Lea  figures  3,  4,  5 et  6 preaertont  lea  rdpartitiona  dedlpaux  nonbres  de  I-iach  M = 0,6C  et  11  = 0,00,  pour 
lea  configurations  avec  et  sana  notour  ot  pour  los  sections  les  plus  caractenatiquPs  (sections  2 et  3). 

Ces  sections  sont  loo  aoulco  concemdes  par  I'effet  a'interaction  entre  rdacteur  et  voilure  ainsi  qu'on 
pout  le  constater,  .rr  oxeaple,  sur  la  figure  7 qui  reprdsente  la  distribution  deCp=  K = 0,80,  avec  et  sans 
rdacteur,  pour  la  aeotior.  n°  5,  i 88  do  I'cnvorgJro  de  I'aile. 


On  pout  observer  quo  I'effet  i' interaction  .’roste  toujouro  faiblo  ot  cu'il  eat  localise  sui-  Ic  bord  a'attaque 
do  l'ailc.  II  n'affeotc  pratiquecont  quo  la  parlie  rdellc  cea  coefficients  Cp.  La  oonparaison  avee  la  distribu- 
tion de  Cp  stationnaire  aans  mo tour  centre  que  les  coefficients  sont  trda  peu  sonsibles  a la  frequence  rdduite 
dans  la  gacae  de  frdquence  dtudide. 

V.  caloul  d'intcraotion  adrodynu’.'d.nuo,  e.'feotud  par  le  uational  Luoht  en  .lul'.tovaartlaboratoriun  d'Anoterdae, 
par  la  nethode  dcs  doublets,  rdfdror.ee  f2.i,  confirae  bien  cettc  evolution. 

Pour  les  cordes  situdes  entre  I'ecplanture  de  I'aile  et  le  rdacteur,  I'effet  du  rdacteur  aui’  I'aile  erde 
une  Idgbre  ddportance  tandis  qua  i'effet  inverse  eat  observe  pour  ieo  sections  au-uelu  du  rdacteur, 

Calcul  et  experience  .sont  en  'w  . accord  et  .nontrent  que  I'interaction  du  rdacteur  sur  i'aile  est  trbs 
locale  et  d'intenoite  suffiscnEcnt  faible  pour  pouvoir  6trc  negligee  dans  les  oalculs  de  flottenent. 


II. 3  - Cffet  du  rdacteur  sur  lui-:*,3ne 

L'etudo  a ".entre  I'importuncc  relative  dcs  efforts  inotationn.ares  sur  Ic  rdacteur  li.i-r.3ne  compares  aux 
efforts  exereds  sur  !•'.  vcilurc  dans  un  nouvenent  de  tarq;agc  de  I'ensenble. 

les  efforts  aui‘  lo  rdacteur  ent  dtd  deterninds  de  doia  facono  : 

- par  difftrsnoo  ues  ..’.ocents  uc  tangage  ddduits  do  I'dvolution  tie  la  frequence  propre  sur  la  naquotte  avec 
ot  sans  rdacteur 

- par  peace  du  rdacteur. 


rl.3.1  - Ddtoraination  a pertir  dea  frequences 

Les  forces  adroa^-oiariiqucs  ddper  ent  li.neairorent  de  l'.on;,lo  d'oscillation  et  ont  jeur  effet  de  faire  varicr 
la  frequence  propre  et  I'anortissencnt  eVu  code  d'osoxlla,tion  do  la  raquette,  Le  i..ontage  coaportant  un  soul  uegr 
de  libertd,  los  forr.iulcs  qui  ox,.rincnt  la  variation  uc  frequenoo  ct  lu  variation  u'amortissenent  en  fonction  du 
coefficient  oomploxe  do  .noaent  sont  oinplcs  et  pouvent  6tre  exploitees  pour  udtoKaincr  co  coefficient  a partir 
des  frdquencea  et  dos  anortisseacn'.s  nesurds  en  soufflorie. 

mg  = ^0  +im0  etant  le  co(,rficient  do  .noaent  autour  de  I'axe  d'oscillation,  pour  une  anplitv'e  angulairo 
de  1 radian,  on  a : Po  ) 

Ps  V 

oil  fo  ot  P sont  los  frdquonces  propres  d'osoillation  sans  vent  ot  avoc  vent 

k un  coeffioiort  incluant  la  masse  gdndralisdo  du  -.ontago  ot  lo  voluac  do  rdference 

Psla  pression  staticuo  uc  i'dooulcaent,  ^ lo  rapport  dos  clialeurs  speoifiquos,  M lo  nombre  do  Imtch. 

Sur  la  figure  8 a dtd  portd  le  coefficient m'grelatif  1 I'ailo  avec  et  sans  rdacteur,  en  fonction  du  .noabro 
de  i.'ach,  pour  differentes  frequences  reduites  d'essai.  La  coui'be  do  la  figure  5 .nontre  I'iaportancc  rol„tivc  des 
forces  sur  le  rdacteur. 
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le  coefficient  k est  donnti  par 

k = m 0 global  - e aile  seule 
mo  globai 

La  valeur  de  K est  de  I'ordre  de  25  et  dei>end  peu  du  noaibre  do  Mach  et  de  la  frequence  rdduite, 

Dans  le  cas  d*un  node  do  vibration  de  "tangage  noteur"  d*une  aile  flexible,  1' amplitude  angulaire  dn  noteur 
serait  plus  inportante  que  celle  de  I'aile  et  la  valeur  relative  du  nonent  dfl  a la  nacelle  s*en  trouverait  accrue, 

Ce  conent  sgit  a I’inverse  des  forces  elastiques,  qui  tendent  toujours  a ramener  la  structure  vers  3a  posi- 
tion d’dquilibre,  et  il  a en  general  un  effet  defavorable  sur  le  flotterent  (effet  le  raideur  negative) , 

Son  importance  s'explique  lar  I'ecoulecent  h I'lntdrieur  de  la  r^elle.  Celle-ci,  en  effet,  se  conporto 
comne  uno  aile  annulaire  dont  lea  oscillations  ongendrent  des  variations  de  pression  inportantes  au  voisinage  de 
la  section  d'entrde, 

II.5.2  - Detemination  par  pesde 

Du  fjiit  de  la  faible  variation  du  coeffioien*  K avec  la  frequence  reduite  et  le  nonbre  de  Ilach  des  nesures 
comparatives  par  pesde  ont  pu  6tre  effectuees  en  slationnaire  et  h.  faible  Vitesse  (incompressible),  dans  une 
petite  soufflerie  de  Laboratoire,  en  faisant  varier  I'incidenoe  entre  - 3®  et  + 3®. 

Ces  nesures  oat  nontre  qua  le  foyer  dtait  situd  trfes  en  avaat  du  rdacteur  et  que  1'  effet  portant  dtait 
localisd  au  voisinage  de  la  section  d'sntree  d'air. 

Le  coefficient  de  nonent  qui  a etd  trouvd  repreaente  la  diffdrenoe  m'o  global  - in’o  aile  seule  et  permet 
d'dvaluer  le  rapport  de  moment  k ddfini  dans  le  paragraphe 

La  valeur  obtenue,  qui  est  dgale  a 23,2  Jo,  recoupe  bicn  la  valeur  ddteminee  a partir  des  variations  de 
frequence  (paragraphe  11.3,1 ). 

Ill  - CALCUL  DES  FORCES  AERODYKAMICUKS  SUR  LS  R5ACTSJR  EM  HACELIg 

L' etude  expdrinentalc  de  I'aile  avec  reaeteur  en  nacelle  a nontre  que  la  plupart  des  problfenes  d’adnoelas- 
ticitd  pouvaient  Stre  traites  en  ca.’.culant  les  forces  aerodynaniques  sur  I'aile  et  sur  le  rdacteur  sdpareaent, 
sans  tenir  conpte  des  interactions. 

Le  calcul  des  forces  sur  I'aile  ou  le  rdacteur  par  la  thdorie  des  doublets  ou  de  la  surface  portante  est 
donnd  dans  les  rdferences  [2J  et  [6]. 

Hous  nous  bomons  ici  a presenter  les  rdsultats. 

ni.1  - 

Le  calcul  foumit  les  coefficients  K. , K_,  K,  et  K.  qui  expriment  la  portance  et  le  nonent  adrodynaniques 
sur  la  nacelle,  par  rapport  a un  axe  de  ’ 34  choisi  h 50  Jo  de  la  longueur  du  cylindre  q\a 

la  schdmatise  dans  le  calcul.  Ces  coefficients  sont  ddfinis  par  1' equation  matrioielle  : 


FL 

n 


= 1 
2 V 


'ki  k2  ■ 

[s/q 

kq. 

, . 

ou  F est  I'affixe  de  la  portance 

Most  I'affixe  du  nonent  autour  de  I'axe  de  rdfdrence  oy  normal  a I'enplanture 

5 et  e6  definissont  le  nouvenent  d'osoillation  : S est  I'affixe  de  la  translation  au  niveau  de  l'a;:e  de 
rdfdrence,  et  oC  I'affixe  de  1' angle  d'osoillation. 

Les  coefficients  K, , ...  K.  sont  complexes  et  ils  ddoendent  du  nonbre  de  Mach  et  de  la  frequence  reduite. 

1 4 

Kj  ddtomine  la  portance  due  a une  oscillation  do  translation  verticalo. 

Kg  ddtemine  la  portance  due  a une  oscillation  de  tangago  autour  de  I'axe  oy. 

Kj  determine  le  nonent  dCl  a uno  oscillation  de  translation  vertioale. 

determine  le  nonent  dd  a uno  oscillation  de  tangage  autour  de  oy. 

III. 1.1  - Conparaison  avec  la  thdorie  de  I'dcoulciiient  bidinonsionnel 

Sn  considdrant  m cylincre  dont  le  dianetre  2R  est  grand  devant  la  longueur  2L  cn  doit  trouver  des  rdsultats 
comparables  aux  cooffioionts  donnes  par  la  thdorie  do  I'dcoulcment  bidimensionnel . Pour  le  vdrifier,  on  a offeotud 
une  application  a un  cylindre  de  24  metros  de  rayon  et  de  1,2  mfetro  de  longueur,  h.  (Jr  faible  (0,05)  et  a Maoh 
0,5.  Les  rosultats  sont  prdsentes  dans  le  tableau  ci-apres.  Leo  coefficients  domes  par  la  thdorie  du  cylindre  no 
sont  pas  rigoureusonent  dgaux  aux  coefficients  de  l'dcoulemf;nt  bidimensionnel,  ce  qui  est  normal  puisque  le  rap- 
port R/L  n'est  pas  «0  , nais  les  diffdrencos  sont  relativemcnt  faiblcs. 


Bidisensionncl 

Theorie  du  cylindre 

0,023  + 3 0,1 

0,016  + 3 0,11 

2,02  - 3 0,29 

2,2  - 3 0,27 

-0,012  - 3 0,05 

- 0,009  - 3 0,C6 

mm 

- 1,OJ6  + 3 0,205 

- 0,995  + 3 0,160 

Pour  cette  coaparaison  Ics  coefficients  K|f  ••♦X.  ont  ete  roacnes  aux  conventions  de  3'ecouleaent  bidi- 
censionnel.  ^ 

ni.J  .2  - Coaparaison  avcc  la  nethodc  des  doublets  de  la  reference  [2] 

Un  calcul  a ete  fait  h Hach  0,8  svj  im  cylindre  ce  rayon  0,0315  n et  de  lonsueijr  0,113  n aniae  d’un  couve- 
cent  de  tangacc  autcmr  du  centre  du  cylindre.  Ics  r^sultats  ont  ete  coapares  avec  ceux  obtenus  par  2.  Roos,  au 
!iia  en  eaployant  In  ndthode  de  la  reference  [2J. 

La  concordance  eot  tres  bonne.  La  seule  divergence,  tres  faible  d'ailleurs,  se  canifeste  sur  Ic  aosont,  la 
position  du  centre  de  poussee  etant  legereaent  differente. 


fethode  des  doublet:. 

Theorie  du  cylindre 

CJr=  0,045 

IM 

3,837  + 0,42  3 

3,915  + 0,387  3 

n 

- 3,23  + 0,144  3 

- 3,48  + 0,138  3 



GJr=  0,0675 

wm 

3,837  + 0,633  3 

3,920  + 0,600  3 

wm 

- 3,243  + 0,210  3 

- 3,48  + 0,210  3 

in.  1.3  - Influence  de  la  pulsation  et  du  Kach 

Sur  un  cylindre  de  rayon  0,054  c et  de  longueur  0,113  n sinulant  la  feme  en  plan  du  reactew  utilise  dans 
les  essais  des  oalculs  ont  ete  effectues  a diffdrents  noabres  de  Kach,  A I'exaoen  des  resultats  on  s'aperfoit 
quo  I'evolution  en  Kach  est  tres  faible,  ainsi  d'ailleurs  quo  I'evolution  en  frequence  reduite,  sauf  pour  certains 
coefficients  nuls  a CJ^=  0 coane  R (K,  ) ou  R (k.,) 


- 5,682 

- 3,809 

- 3,80 

- 3,984 

- 4,09 


- 0,060 

- 0,062 

- 0,0645 

- 0,0648 

- 0,068 


R (K2) 

R (K4) 

5,473 

5,644 

5,924 

- 4,534 

- 4,;iS5 

- 5,306 

,732 

10,31 

2,56 

,803 

10,80 

3,10 

,853 

11,44 

3,80 
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III. 1.4  - Cooparaiaon  avec  I'exp^rleneo 

Pour  valider  la  ach&atisation  du  reactcur  cor  un  cylinarc  crcin.  on  a coapare  on  incoacrcaai'U.e  atatiorr- 
nairc  (la  faiblc  variation  theoiiquo  cea  coePficionta  cn  fonction  do  (Ja  ct  du  Kach  pcraettant  de  ae  iiaitcr  a 
cette  cooparaiaor.)  la  valour  dca  coofficienta  do  portance  et  de  nocent. 


R (Kj) 

R (K  ) 
4 

Calcul 

4,57 

-3.63 

Essai 

4,15 

- 1,80 

lo  coefficient  de  portance  eat  catisc  corrccteaent. 

Une  difference  relative  plua  irportantc  ac  canifcate  aur  le  coefficient  de  rcoent  du  fait  quo  I'axe  de 
reference  eat  relative=ent  prfea  du  foyer  acrodyneaique.  Cette  inprecision  n'apportc  pan  d'erreur  notable  dons  le 
calcul  du  paragraphc  2.3.1  parce  oue,  dans  ce  caa,  I'axe  d'osciUation  eat  situe  tres  en  arriere  du  reactcur. 

La  variation  de  la  partio  inaginaire  du  Bonent  induit  aur  1*  aile  par  le  reacteur  etait  suffiasnsent  faible 
pour  reater  dans  la  diaperaion  dea  nesurcs.  Si  I'on  effectue  la  nfce  coaparaison  a pertir  dea  coefficients  theo- 
riques  on  trouve  wie  pour  une  frequerce  de  50  Ha  a H = 0,8  la  partie  inaginaire  du  nocent  aur  le  reacteur  n’apporte 
qu'vaie  variation  do  2 a 3 5^  sur  le  terae  de  noocnt  global  aile  + reacteur. 

L' influence  dea  tercea  dua  au  reacteur,  aur  lea  forces  e^obales  de  I'aile  a egalecent  ete  calculce. 

Un  calcul  a etc  fait  aur  une  aile,  pour  un  nouveaent  de  tangage  autour  de  I'axe  de  rotation  de  I’aile  a une 
frequence  de  36  Hz.  Lea  valcurs  conpareea  sont  cclles  de  la  partie  reelle  du  coefficient  TIL  dons.,  pir 

M = fV' 


Aile  seule 

Aile 

+ Rdacteur 

H = 

0,6 

Calcul 

- 10,4  10“  ^ 

- 13,4  10“5 

Sssoi 

-9,698  10“^ 

-12,682  10"^ 

Aile  seule 


Aile 

+ Rdacteur 


K = 0,8 

0-3  I _ „ 


Si  I'on  fait  le  rapport 


Hooent  rdysteur 

Kontant  global  (Aile  + reacteur) 


oonne  il  a ete  fait  au  paragrapho  2.3.1,  on  trouve  une  valeur  identique  en  calcul  et  cn  essai 
22,4  a partir  dea  valeurs  theoilques 
23  /J  a partir  dea  valeisra  cxperinentalea, 

Cea  rdsultata  aontrent  la  bonne  concordance  dea  cdthodes  theoruques  de  calcul  dea  foioea  aerodynaciques 
sur  lea  ailes  et  sur  le  cylindso  avec  lea  valeurs  experijnentales,  pour  un  terme  diroct  de  conent  de  tangage  a une 
frequence  ddterndnee. 


Pour  Bontrer  1' influence  defavorable  du  reacteur  sur  lea  limites  de  flotteaent,  la  naquette  a ete  aontde  a 
I'emplanturo  sur  un  support  a deux  degrda  de  libertd  (tangage  et  roulis)  qui  pemet  d'obtonir  un  flotteaent  dans 
le  doaaine  de  la  aoufflerie. 

Deux  configurations  ont  dtd  essaydes  s 

- une  configuration  dans  laquellc  I'aile  o''t  dquipde  du  rdacteur  en  nacelle,  avec  dcoulenenl  intense, 

- une  configuration  dans  laquellc  le  rdac  eui'  en  isacelle  est  ronplacd  par  ssn  corps  fuseld  de  mfinso  masse, 
mfine  inertie  et  mfime  centre  de  gravite,  de  faqon  h.  no  pas  modifier  les  frequences  propreo,  les  masses 
gdneralisdes  et  les  ddformdos  des  deux  codes  propres. 

Leo  dvolutions  dans  le  vent  des  dessx  frequences  propres  et  des  amortissoments  ont  dtd  ddtermineos  a partir 
des  enregistremento  de  rdponse  de  la  maquette  h la  turbsslenco,  par  la  mdthodo  de  Cole  (Rdf.  [4]).  Cette  methode 
pcrtset  d'obtonir  les  rdponres  inpsslaionnelles  de  la  maquette  (fig.  10). 

Los  limites  de  flottemont  ont  dtd  determindeo  par  extrapolation,  b partir  des  frequences  ot  des  amortis- 
sements  nessurdo  aloro  quo  la  maqsiette  etait  encore  loin  du  domaine  d'instabllitd  par  la  mdthode  do  Zimmerman 
[5].  Dans  cette  mdthode,  on  determine  un  paran&tro  de  stabilitd,  P,  qui  s'onnsile  k I'entrde  en  flottemont  et  qssi 


3^ 


erolue  suivant  uno  parabolc,  en  fonction  de  la  presaion  gdadratrico  pi  , pour  uae  evolution  a Hach  constant. 
Le  para:^tre  est  donnd  par  : 


(3l  - + -V  (3,(3i  (OLtlskM  + a 

2m  2 J 


(Jz  - CJ?  +2  ^ 


ou  (J'l  et  Oz  sont  les  deux  pulsations  propres,  et  Pg  = Oz  t.9(^tet  °^2  dtani  les  anortissenents 

rdduits. 

Les  courbes  F ( ) deteminces  a partir  des  frequences  et  des  anortissenents  thdoricues  et  ezocrisentaux 

ont  4tS  repr^sentdes  sur  les  figures  (ll)  et  (12)  pour  les  noabres  de  llach  0,60  et  0,60. 


On  voit  que  I'influence  relative  du  reacteur  est  plus  Icportante  h Kech  0,60  qu*a  0,80,  ce  qui  est  logique, 
puisque  le  cooent  dfl  an  rdacteur  est  peu  sensible  an  Kach,  tondis  que  le  cocent  dO  a I'aile  subit  un  accroisse- 
Dent  considdrable  entre  Kach  0,60  et  Kach  0,80. 

Les  courbes  expdrinentales  ont  etc  lissees  par  unc  sethode  de  coindres  carrds  autour  d'une  parabole. 

Dans  I'ensenble,  les  courbes  thcoilques  et  expdriacntales  presentent  un  decalage  qui  correspond  a une 
erreur  de  I'ordre  do  10  jS  sur  la  pression  generatrice. 

On  constate  que  la  pression  d'entree  en  flottencnt,  qui  correspond  a F = 0,  est  plus  faiblc  aver  le  reac- 
teur qu'avec  le  corps  fusele  sans  ccouleDent  interne,  d'environ  15  Cette  influence  ddfavorable  est  correctenent 
prevue  par  la  tlieorle. 


V - COKCmSIOH 

La  cesure  des  forces  adrodynaaiques  instatioranaires  sur  une  naquette  dquipde  d'un  rdacteur  on  nacelle,  en 
subsonique,  a nontrd  que  la  portance  et  le  noaent  induits  sur  le  rdacteur  l-i-cfine,  par  ses  propres  rscillations, 
dtaient  tree  t’jportants. 

Ces  forces  ne  doivent  pas  fitre  ndgjigdcs  dans  les  calculs  de  flottenent  des  avions,  car  elles  pouvent 
abaisser  sensibleaent  les  liaites  de  flottezxsnt  par  effet  de  "rigiditd  negative”  sur  les  codes  de  torsion  de 
voilure  et  surtout  de  "tangage  de  coteur". 

On  a pu  observer  dgaleaent  ; 

- que  les  forces  adrodynaaiques  sur  le  rdacteur  sont  pea  sensibles  au  Ilaoh  et  h la  frequence  rddxiite, 

- que  1' interaction  adrodynanique  rdacteur-yoilure  pout  6tre  negligee. 

En  se  fondant  sur  cette  demiere  recarque  on  a utilise  une  theorie  qui  pcrcet  de  trailer  sdpareDent  le 
coteur  en  nacelle  en  le  schematisant  par  un  tube  cylindrique  a poroi  mince  avec  dooulement  subsonique  de  r.6me 
Vitesse  a I’intdrieur  et  a I'oxtdrieur. 

La  coaparaison  des  rdsultats  theoriques  et  expdriiaentaux  est  satisfaisante. 
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CALCUIATION  OP  AEROUniAMIC  LOADS  OH  OSCILUTIHO  HIHG/STORE  COHBIHATIOKS  IN  SUBSONIC  FLOH*^^ 


B.  Beimekers,  R.  Roos  and  R.J.  Zwaan 
National  Aerospace  Laboratory  NLR 
Anthony  Fokkerweg  2 
Ansterdaa-1017 
The  17etherlands 

SUMMARY 

A method  for  the  calculation  of  aerodynamic  loads  on  wing-store  configurations  oscillating  in  subsonic 
flow  is  presented.  In  this  method  the  linearized  ecpiation  for  subsonic  conpressibl . flow  is  treinsformed  into 
two  sets  of  integral  equations  for  the  steady  and  a superimposed  tmsteady  flow  field. 

The  wing  loads  are  represented  hy  dipole  distributions  (wing  thickness  is  neglected)  and  the  store 
loads  by  source  distributions.  Discretizing  these  distributions  into  lifting  lines  and  source  panels  of 
constant  strength  results  into  a set  of  ailgebraic  equations.  These  are  solved  for  the  unknown  distributions 
by  forcing  the  flow  to  be  tangential  to  the  surfaces  of  the  jscillating  wings  and  bodies  in  - set  of  control 
points. 

The  solution  enables  the  calculation  of  pressure  distributions  on  the  wings  and  storuti  and  of  general- 
ized aerodynamic  coefficients.  Calculated  results  are  presented  and  compared  with  experiments. 
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LIST  OP  SYMBOLS 


influence  coefficient 
velocity  of  soiind 
pressure  coefficient 
reduced  frequency 
wave  number 

components  of  the  normal  vector  n 

Mach  number 

number  of  panels 

normal  vector 

velocity  vector 

position  vector 

disteuice  between  collocation  point  i and  the  centroid  of  panel  j 
displacement  vector 
surface  of  wing  or  body 


t longest  panel  diagonal 

U,  V,  W components  of  the  velocity  vector 
X,  y,  z right-hainded  cartesian  co-ordinate  system 
X , y , z components  of  the  displacement  vector 

O'  0 6 

ACp  pressure  jump  between  lower  and  upper  wing  surface 

ACp/A  normalized  pressure  jump  over  the  wing  surface 

a angle  of  attack 

* (1  - 

Y specific  heat  ratio 

6 arigular  co-ordinate  on  tip  tank 

U dipole  strength  on  the  wing  camber  surface 


(degrees) 


(degrees) 

(lVt) 


■'^)This  investigation  was  carried  out  under  contract  for  the 
Scientific  Research  Branch,  Air  Materiel  Directorate,  RNLAF. 
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right-handed  cartesian  co-ordinate  system 

source  strength 

velocity  potential 

disturbance  velocity  potential 

oscillation  frequen;^ 

dipole  strength  on  the  wake  surface 


Subscripts 


referring  to  the  body 


" " camber  surface 
" " lifting  surface 

" steady  (n  = O)  and  unsteady  Gi  = l)  flow  field 
" the  W3u:e 

" the  steady  reference  position 
" " unsteady  flow  field 
" " free-stream  condition 


Superscripts 


B referring  to  the  body 

D " " " lifting  surface 

**  " " " time-independent  quantity 

1 lUTROrUCTION 

Many  modern  military  aircraft  are  equipped  with  ever  larger  stores  under  the  wing  and  at  the  wing 
tips.  A problem  which  arises  in  the  investigation  of  the  aeroelastic  characteristics  of  such  aircraft 
(flutter,  gust  response  and  manoeuvring  loads)  is  how  to  obtain  reliable  estimates  of  the  aerodjoiamio  loads 
introduced  by  these  stores.  Important  questions  in  this  context  are;  What  does  the  presence  of  a store  do 
to  the  unsteady  aerodynamic  loads  on  the  wing  and  what  are  the  aerodynamic  locuis  on  the  store  itself  in  the 
presence  of  the  wing? 

Until  some  time  ago  the  aerodynamic  loads  on  these  stores  and  their  influence  on  the  wing  were  neglect- 
ed or  were  approximated  for  instance  by  using  slendir  body  theory  only.  However,  to  the  knowledge  of  the 
authors  the  value  of  this  application  of  slender  body  theory  has  not  yet  been  tested  conclusively  in  a com- 
parison between  theory  and  experiment. 

A more  thorough  attempt  to  determine  wing-body  interference  effects  on  harmonically  oscillating  con- 
figurations was  recently  made  by  Kalman,  Hodden  and  Giesing  (Ref.l).  Using  the  doublet  lattice  method  they 
assumed  the  body  to  be  a cylindrical  ring  wing  eind  thus  were  able  to  apprcximate  the  interference  effects 
on  the  wing  due  to  the  proximity  of  the  body.  Although  such  calculations  might  give  a reasonable  prediction 
of  these  interference  effects,  it  is  clear  that  the  diversion  of  the  flow  by  the  body  is  not  accounted  for 
and  that  forces  on  the  body  cannot  be  determined  in  this  way.  To  alleviate  this  problem  they  incorporated 
slender  body  theory  (Ref. 2),  and  to  overcome  the  problem  of  large  systems  of  equations  they  introduced  an 
image  system  in  the  body  representing  the  interference  effects  (Ref.3).  This  way  a method  was  obtained, 
capable  of  handling  configurations  consisting  of  lifting  surfa'^cs  and  bodies,  which  gives  satisfactory  re- 
sults in  the  calculation  of  aerodynamic  coefficients. 

A few  years  ago  an  extensive  research  program  was  started  at  the  NLR  to  investigate  airloads  on 
oscillating  wing-store  configurations.  Recently  detailed  pressure  measurements  were  completed  on  a wing- 
tiptank-store  configuration,  as  reported  by  Renirie  (Ref.4).  Another  part  of  the  research  program  was  the 
development  of  a method  for  the  calculation  of  aerodynamic  loads  on  oscillating  wing-store  configurations  in 
subsonic  flow.  The  purp-.se  of  such  a method  was  twofold; 

1 calculation  of  air  forces  on  wing-store  configurations  for  aeroelastic  analysis; 

2 support  of  wind  tunnel  measurements  through  the  calculation  of  detailed  pressure  distributions. 

~ The  method  presented  here  was  set  up  in  an  effort  to  meet  the  above  requirements.  It  combines  the 
elements  cf  four  existing  methods  into  one  general  method  for  the  calculation  of  steady  emd  unsteady  press- 
ure distributions  on  wing-body  configurations.  The  steady  flow  field  is  described  with  the  vortex  lattice 
method  and  the  panel  method  of  Hess  and  Smith  (Ref.5).  The  unsteady  field  is  modelled  with  the  doublet 
lattice  method  and  a new  panel  method,  which  is  partially  based  on  the  work  of  Hess  (Ref.6)  on  the  calcu- 
lation of  acoustic  fields  about  arbitrary  bodies.  In  the  following  the  part  of  the  method  describing  the 
steady  flow  field  will  be  called  the  KLRS-method,  while  the  part  in  which  the  steady  and  unsteady  fields  are 
combined  will  be  referred  to  as  the  NLRI-method. 

This  paper  gives  a description  of  the  method  and  presents  some  preliminary  calculated  results  on  a wing- 
tiptank  configuration  in  comparison  with  experimental  data. 

2 OOTLIHE  OP  THE  THEORY 

The  present  description  of  unsteady  airloads  on  wing-store  combinations  involves  airforces  on  rather 
thin  lifting  surfaces  (wing,  pylons,  tail  surfaces)  and  on  streamlined  closed  stores,  both  influencing 
each  other  through  interference  effects. 

For  the  calculation  of  aerodynamic  forces  on  lifting  surfaces  several  effective  methods  are  known.  In 
the  theory  presented  here  the  lifting  surfaces  are  treated  with  the  doublet  lattice  method  (DL-method),  since 
it  is  very  flexible  in  handling  different  configurations  (Ref.?).  Recently  the  reliability  of  the  DL-method 
was  shown  again  by  Roos  and  Zwaan  (Ref. 8)  in  calculations  of  the  unsteady  aerodynamic  loads  used  in  the  de- 
termination of  the  flutter  behaviour  of  a wing-pylon-naoelle  wind  tunnel  model. 

Considering  the  large  experience  with  the  HLR-peinel  method  for  steady  flow  (see  for  example  reference  J), 
the  unsteady  flow  field  about  the  stores  was  chosen  to  be  described  with  an  unsteady  source  panel  method. 
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This  ioplies  that  the  surface  of  the  hody  is  discretized  into  panels  each  containing  a time-varying  source 
distribution  of  constant  yet  unknown  strength  over  the  pemel.  The  wake  behind  the  store  is  assumed  to  be 
represented  with  an  extension  of  the  lifting  surface  into  the  store.  This  is  a rather  crude  simplification, 
but  appropriate  experimental  evidence  is  lacking.  However,  as  in  steady  calculation  methods  this  simplifi- 
cation leads  to  good  results,  it  seems  to  be  a reasonable  assumption.  The  boundary  conditions  are  enforced 
in  the  control  point  of  each  panel  on  lifting  surfaces  and  stores. 

Let  the  inviscid  flow  around  an  oscillating  arbitrary  wing-store  configuration  be  described  with  a 
velocity  potential 

<3CiytZi't>  = + HoZ  + P -ix,y,z,t>  (l) 

in  which  U„,  eind  W„  are  the  components  of  the  free-stream  velocity  vector  cuid  9 is  the  disturbance 
velocity  potential.  The  use  of  a small  disturbance  potential  implies  the  assumption  that  the  body  does  not 
disturb  the  free  stream  significantly.  For  the  region  of  the  stagnation  points  the  validity  of  this  assump- 
tion will  have  to  be  verified  afterwards.  Assuming  that  V„  and  Wo,  eire  small  as  compared  to  Ua>  and  that  <j>  is 
small  in  comparison  with  ILx  the  general  equation  f"'r  ^ reduces  to  a linearized  eqxiation  for  9 

2 1 

P <^xx  ’yy  ^ ’zz  " “ <^tt  ' ^ 1>xt  = ° (2) 


/ 2 2 

with  = Ug/ajo  and  p = 1 - The  Mach  number  is  taken  to  be  constant  throughout  the  flow  field,  which 
consequently  makes  the  value  of  the  method  doubtful  in  the  near-treinsonic  Mach  range. 

The  surface  of  an  oscillating  arbitrary  configuration  can  be  described  in  general  with  the  functional; 

S 


<x,y,z,t>  = 0 (3) 

The  boundary  condition  requiring  the  flow  to  be  tangential  to  the  moving  surface  is  then  given  by; 

Dt  ’ at 


= ■1^+  q .VS  = 0 


on  S = 0 


(4) 


In  Apperidix  1 the  general  solution  of  equation  (2)  is  derived  for  a harmonic  motion  of  small  amplitude. 
In  principle  this  solution  includes  the  effects  due  to  the  thickness  and  the  wake  of  the  wing  in  about  the 
same  way  as  does  the  theory  presented  recently  by  Morino  and  Kuo  (Refs. 10,  11).  However,  it  is  also  agued 
there  that  in  practice  an  accurate  consideration  of  the  wake  effects  complicates  the  calculation  substan- 
tially. Recalling  the  purpose  of  the  method  and  realizing  that  lifting  surface  theories,  although  they  ne- 
glect thickness  effects,  are  quite  accurate  already  as  demonstrated  in  references  ^ and  12,  the  wing  is 
assimed  to  be  infinitely  thin.  Thus,  a method  is  obtained  of  which  the  set-up  is  described  earlier  in  this 
section.  In  Appendix  1 this  set-up  is  formulated;  the  pertinent  expressions  for  the  velocity  potential  9 
are  repeated  here. 
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in  which 


and  for  the  lifting  surface 
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in  which 
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and 
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(10) 


In  the  above  expressions  the  disturbance  potential  <p  is  written  in  terms  of  a source  distribution  of 
strength  cr^  on  the  body  surface  and  a pressure  dipole  Sistribution  of  strength  ACp^  on  the  thin  lifting 
surface  S^.  The  subscript  n can  have  the  values  0 and  1,  indicating  the  steady  and  \msteady  potential. 

The  desorip)tion  of  the  surface  of  an  oscillating  arbitrary  body  is  not  as  simple  as  it  might  look  from 
expression  (j),  because  the  functional  form  of  S is  not  known  in  general.  However,  often  the  time-dependent 
position  of  the  individual  points  of  the  surface  can  be  expressed  in  terms  of  a reference  position  and  a 
time-dependent  amplitude  r^.  Assuming  that  this  amplitude  is  small  such  that  cp^«l^x,  the  boundary  con- 
dition on  the  surface,  as  given  by  (4),  can  be  divided  into  a ste^y  and  an  unsteady  part,  as  shown  in  Appen- 
dix 2.  For  the  steady  cp^  this  results  in  the  classical  condition*' 


(q„  + V<P„) 
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< x,y,z  > = 0 


(11) 


the  normal  vector  n is  directed  outside  the  body. 


44 


fe'-: 


and  for  the  unsteady 

V<?i  I .n  - - (q„+ V<P^).(t  V;  + mv^y“+  + n.[(r“7)V9^  | = 0 


on  S <.  x,y,z,t  > = 0 


(12) 


with  r|J  being  the  amplitude  of  the  displacement  vector  having  components  i",  y**  and  z“.  The  oomponents  of 
the  normal  vector  n on  the  tody  in  its  stationary  position  are  given  by  reap. t®,  m and  n.  For  the  thin 
lifting  surfaces  the  unsteady  condition  reduces  to  a simpler  form: 


Vfi  . n - iuTg.n  - ® 


(13) 


on  < x,y,z,t  > = 0 


From  (12)  It  is  clear  that,  while  the  linearized  equation  (2)  allows  for  a disconnection  of  the  steady  and 
unsteady  fields,  as  implemented  in  (5),  the  boundary  condition  restores  this  connection  again  througji 
second  order  terms  involving  first  and  second  order  derivatives  of  the  steady  disturbance  potential  qp  . In- 
cluding these  terms  is  necessary  to  obtain  a proper  representation  of  the  effects  due  to  diversion  orthe 
flow  as  a result  of  the  thickness  of  the  body.  Morino  and  Kuo  (Refs.  10,  ll)  do  not  indicate  this  connection 
in  their  work  on  oscillating  wings  with  thickness. 

Substitution  of  the  expressions  (6)  throu^  (lO)  for  <p  in  these  boundary  conditions  leads  to  the 
following  sets  of  integral  eqpiatioiB  for  resp.  the  steady  and  unsteady  flow  field: 
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(17) 


Equations  (I4)  and  (15) determine  the  steady  flow  field  on  the  body  sind  the  lifting  surface,  while 
eqxiations  (16)  and  (l?)  concern  the  unsteady  flow  field  superimposed  on  this  sxeady  field.  The  interference 
effects  between  the  body  and  the  lifting  surface  are  represented  in  the  secord  integral  of  equations  (14) 
and  (16)  and  in  the  first  integral  of  equations  (15)  and  (I7). 

To  solve  these  integral  equations  for  the  unknown  source  and  dipole  distributions,  these  distributions 
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are  discretized.  The  surface  of  the  body  and  the  lifting  surface  are  respectively  divided  into  Mg  and  Mg 
panels.  On  the  body  each  panel  is  assiuned  to  contain  a source  distribution  of  constant,  yet  uijaiown  strength. 
The  panels  on  the  lifting  surface  are  taken  to  contain  each  a lifting  line  of  constant,  also  unknown, 
strength,  along  the  I/4  -oh'^rd  line  of  the  panel.  Further,  each  panel  possesses  one  collocation  point  in 
w.uch  the  boundary  condition  is  imposed.  That  is,  in  these  points  the  flow  is  forced  to  be  tangential  to  the 
surface.  With  this  discretization  the  sets  of  integral  equations  are  transformed  into  two  sets  of  (Ng+Ng)  al- 
gebraic equations  for  the  steady  and  unsteady  part  of  the  flow: 
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in  which C''  and  ACp  are  the  unknown  source  and  lifting  line  strength  of  each  panel  of  the  body  and  the 
lifting  surface.  The  term  is  the  prescribed  normal  velocity  in  the  collocation  point  of  the  panel  of 

either  the  body  or  the  lifting  sui-face.  The  influence  coefficient  A^J  represent  . the  normal  velocity  induced 
by  the  i^*^  panel  with  a source  distribution  or  lifting  line  of  unit  strength  -..1  the  collocation  point  of  the 
panel.  Their  methods  of  calculation  are  indicated  in  section  3.  The  numerical  solution  of  the  sets  of 
algebraic  equations  is  given  in  section  4 of  this  paper. 

After  the  steady  source  and  lifting  line  strength  have  been  calculated,  the  steady  part  of  the  velocity 
in  the  collocation  points  can  be  determined.  The  second  derivatives  of  9^  are  obtained  by  differentiation  of 
the  expressions  for  the  velocity.  Next,  these  values  are  used  to  determine  the  right-hand  side  of  the  unstea- 
dy set  of  equations,  after  whicli  the  unsteady  flow  field  can  be  solved  as  well. 

For  the  lifting  surface  the  strength  of  the  lifting  l.ne  is  taken  to  be  equal  to  the  pressure  jump  over 
the  surface  in  the  middle  of  the  1/4  -chord  line  of  the  panel.  The  steady  and  unsteady  pressure  distributions 
on  the  body  follow  from  the  following  expressions,  derived  in  Appendix  3: 
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in  which 


(21) 


be  divided  into  two 


Figure  1 gives  a schematic  outline  of  all  necessary  calculations. 

3 CALCULATION  OP  THE  INFLUENCE  COEFFICIENTS 

For  both  the  steady  and  unsteady  flow  field  the  influence  coefficients  A^'-’  .... 

groups:  one  A^^  in  which  the  pa.iel  lies  on  a lifting  surface  and  a second'^A^^®  having  the  panel  on 
the  body.  For  the  steady  flow  field  the  calculation  follows  exactly  the  lines  o)"  two  well-known  methods.  The 
coefficients  Aj®  are  obtained  with  the  expressions  of  the  panel  method  of  Hess  and  Smith  (Ref. 5),  after  the 
compressibility  factor  d has  been  removed  temporarily  from  (?)  with  a transformation  of  the  Prandtl-Glauert 
type.  The  coefficients  AJ®  are  evaluated  according  to  the  formula+ions  for  the  vortex  lattice  method  given 
in  reference  13.  The  calculation  of  the  unsteady  coefficients  A^®  does  not  introduce  any  new  pi'cbloms  either. 
They  are  simply  obtained  from  the  formulations  of  the  doublet  lattice  method  given  in  reference  2. 

For  the  unsteady  coeffi  'lents  Aj®  represented  by  the  integral 
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a new  calculation  method  is  developed  which  in  part  is  based  on  the  panel  method  of  Hesc  (Ref. 6)  for  the 
determiriation  of  acoustic  fields  around  arbitrarily  shaped  bodies.  r 

To  evaluate  the  above  integral  five  different  regions  are  disting  ashed  based  on  the  ratio  ^ being  the 

distance  between  collocation  point  1 and  the  centroid  of  the  panel  j o'  ei  the  loigest  diagonal  ^ of  the 

panel.  In  these  regions  the  integral  is  calculated  as  follows: 
r 

1.  — ^ 4,  the  point  1 lays  far  away  from  panel  j;  the  source  distribution  on  the  panel  is  approximated  with 
a point-source  at  the  centroid  of  the  panel. 

r 

2.  4 > “ ^3-9,  fhe  point  1 lays  at  fairly  large  distance  from  panel  j;,.the  source  distribution  is  approx- 
imated with  a two  term  multi-pole  expansion  in  powers  of  (kj^t)  and  (— '•)  around  the  centroid. 

r c 

3.  1.9^  “ ^ 1,  fl’-e  point  i lays  at  fairly  small  distance  from  panel  j;  the  sajne  approximation  as  in  2., 
with  three  terms  teiken  into  account  instead  of  two. 
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4.  1 > ^ > 0,  the  point  i lays  close  to  the  centroid  of  panel  j;  the  integrand  is  approxinated  with  a 
one-dinensional  Taylor  series  expansion  in  poweM  of  (ikj^r). 

r 

5.  ^ = 0,  the  point  i lays  at  tne  centroid  of  panel  j;  the  integrand  is  expanded  similarly  as  in  4. i tut 
now  the  expressions  are  simpler. 

A detailed  description  of  the  above  approximations  and  the  resulting  expressions  is  given  in 
reference  14. 

In  their  work  Horino  and  Kuo  (Refs.  10,  11)  have  not  introduced  such  a ccmplicated  sclieme  to  calculate 
the  unsteady  influence  coefficients.  For  all  Vcilues  of  r^.  they  take  the  exponent  as  a constant  in  front  of 
the  integral  and  then  are  left  with  an  integration  which  is  identical  to  the  one  belonging  to  a steady 
source  eind  doublet  distribution.  Althou^  this  gives  a good  approximation  for  larger  distances,  it  can  be 

expected  that  for  small  r—  the  result  will  not  be  very  accurate.  On  the  other  hand  the  possibility  exists 
that  this  inaccuracy  might  be  compensated  hy  the  use  of  curved  panels  as  recently  proposed  by  Ci.en,  Suciu 
and  Horino  (Ref. 15). 


4 ITOMERICAL  SOLOTIOH 

After  the  influence  coefficients  a’’^  have  been  calcnilated,  the  two  sets  of  equations  given  by  (l8)  will 
have  to  be  solved.  Making  a distinction  in  these  equations  between  the  body  and  the  lifting  surface,  (l8) 
can  be  rewritten  in  matrix  form; 
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in  which  the  matrix  is  peirtioned  into  four  submatrices,  eacdi  containing  the  influence  coefficients  of  a pair- 
ticular  type  of  influence  between  nanels  on  the  body  and/or  the  lifting  surface.  The  advantage  cf  this  par- 
titioning iS  that  the  suhmatrix  [ A“  1,  describing  the  influence  between  the  panels  on  the  body,  is  strongly 
diagonal-domin^t.  The  sitoatrix  [ possesses  a similar  structure,  although  less  pronounced.  The  two 
submatrices  [A'^^3  a"^[  not  have  a specific  structure,  but  in  general  their  elements  will  be  small 

compared  to  the  diagonal "terms  in  the  other  two  submatrices. 

Baaed  on  these  expected  properties  with  regard  to  the  general  structure  of  the  matrix  equations  the  follow- 
ing iterative  procedure  is  used  to  solve  .hem. 

First  the  set  containing  [ A®®  ] is  solved  with  an  iterati;^  Causs-Seidel  process.  Then  the  residue  is 
18  determined,  which  is  used  to  solve  the  set  containing!]  A"  J through  a direct  Crout-process.  The  residue 
IS  determined  again  and  the  process  is  started  from  the  be^nning.  This  iteration  process  is  repeated  several 
times  until  the  increment,  which  during  each  run  is  added  to  the  solution  of and  dCp  has  become  smaller 
than  a certain  value.  To  overcome  the  problem  of  numerical  under-  or  overshoots  in  the  first  few  iteration 
steps,  the  calculation  of  the  residues  is  performed  with  a relaxation  factor  on  the  correction  term. 


5 CALCUUTED  RESULTS 


To  test  the  method  a configuration  is  chosen  on  which  extensive  steady  and  unstead.,  pressure  measure- 
ments have  been  performed  at  HLR  (see  reference  4).  The  configuration  which  is  shown  in  figure  2,  consists 
of  a tapered  wing  with  a tip  tank  and  a removable  pylon-store  attachment. 

The  panel  distribution  used  for  the  calculations  on  the  wing-tip  tank  configuration  is  designed  as 
follows.  The  wing  is  divided  into  10  chordwise  strips  of  10  panels  each.  The  last  strip  is  assumed  to  con- 
tinue up  to  the  tip  tank  axis,  such  that  part  of  these  pemels  fall  within  the  tip  tank.  This  strip  does  not 
contain  collocation  points  and  the  strength  of  the  singularities  is  taken  to  be  equal  to  the  once  of  the 
strip  next  to  it.  The  tip  tank  is  panelled  into  27  octeigonal  sections  perpendicular  to  the  axis,  totalling 
a number  of  216  panels. 

5a  Steady  results 

Figures  3 and  4 show  some  chorswise  distributioiB  of  the  steady  pressure  jump  ACp^  across  the  wing  of 
the  wing-tip  tank  configuration  for  M,,  = 0.45  and  a = O”.  For  both  sections  the  agreement  between  the  results 
of  the  ULRS-method  auid  the  measurements  is  satisfactory.  The  leading  edge  peak  visible  in  both  figures  is 
caused  by  the  fact  that  the  wing  profile  has  a drooped  nose  which  extends  over  the  first  40  % of  the  chord. 
Corresponding  pressure  distributions  on  the  tip  tank  are  given  in  figures  5 and  6.  Figure  5 shows  a compa- 
rison between  calculated  and  measured  pressure  distributions  in  the  direction  along  the  tip  tank  axis  at  a 
station  just  above  the  wing  reference  plane  (6  = 202.5  ). 

The  leading  edge  peak  induced  by  the  v;ing  is  clearly  visible  in  the  theoretical  results  but  is  too  local 
to  be  discernable  in  the  measurenients.  Further  p.way  from  the  wing  this  peak  has  disappeared  also  in  the  cal- 
culated results  as  can  be  seen  in  figure  6,  which  gives  the  pressure  on  the  outside  of  the  tip  tank. 

5b  Quasi-steady  and  imsteady  results 

Unsteady  results  have  been  calculated  for  a pitching  oscillation  around  an  axis  at  I5  / of  the  root- 
chord,  with  a reduced  frequency  k = O.305  based  on  the  wing  serai-span.  In  figures  7 and  8 a comparison  is 
given  for  the  same  wing  sections  as  used  in  thf  steady  case.  The  chordwise  distributions  of  the  unsteady 
pressure  jump  AC_^/a  calculated  with  the  NLRI-method  are  in  good  agreement  with  the  experimental  data  for 
both  the  real  ana  imaginary  part.  In  this  comparison  all  distributions  are  normalized  with  the  same  amplitude 
A.  Also  shown  in  these  figures  is  the  quasi-steady  distribution  based  on  the  difference  of  the  results  of  the 
HLRS-method  for  a = 2°  and  0°.  Apparently  the  reduced  frequency  of  k = O.3O5  is  still  too  small  to  give  a 
significant  difference  between  the  quasi-steady  results  and  the  real  part  of  the  unsteady  distribution. 

A signifiCcint  difference  is  found  indeed  in  these  figures  if  the  theoretical  results  of  the  wing  with  and 
without  tip  tank  are  compared.  Clearly,  the  presence  of  the  tip  tank  causes  the  pressure  level  on  the  wing 
to  rise.  This  interference  effect  is  most  prominent  in  the  vicinity  of  the  tip  tanx,  but  is  still  noticable 
in  section  3. 
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Contrary  to  the  theory  1;hese  unsteady  calculations  with  the  IHJ)I-aethod  are  perfomed  under  CQission 
of  the  second  derivatives  of  9^  in  the  boundary  condition  (12)  and  the  foraula  (20)  for  the  pressure  dis- 
tribution on  the  body  (tip  tank).  During  the  calculations  it  was  found  that,  as  is  assuned  in  the  theory, 
the  coaputer  progran  of  the  KLBS-aethod  indeed  generates  saall  perturbation  velocities  on  the  body.  However, 
the  second  derivatives  of  9^  do  not  coae  out  saall,  but  instead  are  so  large  that  they  doainate  the  solution 
of  the  flow  field  coapletely.  The  reason  for  this  is  not  yet  coapletely  understood,  but  it  is  believed  to  be 
inherent  to  the  coarseness  of  the  panelling.  That  these  second  derivatives  of  9^  should  be  saall  indeed  is 
partially  confiraed  by  the  fact  that  oaitting  then  leads  to  the  excellent  coaparison  with  the  aeasureaents 
on  the  wing  where  they  enter  the  calculation  throu^  interference  with  the  tip  tank.  Preliainary  calculations 
of  the  unsxeedy  pressure  on  the  tip  tank  indicate  that  the  second  derivatives  can  not  be  neglected  completely. 

6 CONCLUSION 

A method  is  presented  which  enables  the  calculation  of  unsteady  pressure  distributions  on  wing-body 
configurations.  Since  the  unsteady  flow  field  depends  on  the  steady  flow  field  on  which  it  is  superimposed, 
this  steady  flow  field  is  calculated  as  a part  of  the  nethod.  The  results  of  the  steady  method  seem  to  be 
in  agreement  with  the  measiurements  except  for  the  calculation  of  the  second  derivatives  of  the  disturbance 
potential.  This  problem,  which  is  probably  caused  by  a too  coarse  discretization  of  the  body,  is  being  in- 
vestigated at  the  moment. 

^eoretical  unsteady  pressure  distributions  on  the  wing  show  a very  good  agreement  with  the  experimental 
ones,  while  froa  calculations  for  the  wing  only  it  is  found  that  the  interference  effects  of  the  tip  tcink  on 
the  Wing  are  not  of  a negligibli  magnitude.  Preliminairy  results  on  the  body  seem  to  indicate  the  importance 
of  an  accurate  calculation  of  the  second  derivatives  of  the  disturbance  potential. 
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APPEHDIX  1 

THS  GEHERAL  SOLUTIOH 

The  disturhance  velocity  potential  ^ <ix,y,z,t  > satisfiee  the  eq'oation 

2 ' 1 / \ 

p 9 +9  +9 5'  9 - 2 — 9 * = 0 (Al.l) 

^yy  ^zz  2 '^tt  a_  xt  ' ' 

^to 

in  which  » Ho/®®  = 1 - 1^.  Assnning  harmonic  variations  in  tine  9 can  he  written  as; 

9 < *iy»z«‘t>  = 9q  < *ty»2  > + 9j  < x,y,2  > e^*^  (A1.2) 


with  9 being  the  steady  part  and  9j^e^  the  unsteady  part* 
Using  also  the  substitutions: 


9„  < x,y,z  > = 9„  < Xfy.z  > e 


ik„M^ 


(n  = 0.1) 


X*  = X , y*  = Py  , 
equation  (Al.l)  can  he  transforaed  into 


+ 9 “ + 9„**  + k%>!!  = 0 


(n  = 0.1) 


(AI.4) 


(AI.5) 


A general  solution  9 of  this  equation  can  he  written  in  terns  of  a source  distribution  over  the  body 
and  a dipole  distribution  over  the  wake.  To  avoid  numerical  instabilities  near  the  trailing  edge,  as  indic- 
ated by  Rubbert  et  al  (Ref.l6),  the  dipole  distribution  over  the  wake  should  be  extended  to  the  cawber  sur- 
face inside  the  body.  Transforming  the  result  back  to  an  expression  for  9 , one  obtains: 


< x,y,z  > 


Sr 


4tr  ji^n  <^C'”C*^C 


^ // 


ik  M (x-?,.)  , 
n <D'  ’if'  6 


(A1.6) 


for  n = 0.1  and 


= [(x-5)^  + P^(y-ri)^  + p^(z-c)^J 


In  (A1.6)  C represents  the  strength  of  the  sowce  distribution  on  the  body  stirface,  p and  A9  the 
dipole  distribution  on  the  camber  surface  and  the  surface  of  the  wake. 

The  calculation  of  the  integral  over  the  wake  surface,  is  rather  difficult  in  practice.  Obtaining  ana- 
lytical expressions  is  hampered  by  the  fact  that  the  stationary  position  of  the  wake  is  not  known  a priori 
and  thus  vn.ll  have  to  be  determined  first  with  some  kind  of  iteration  procedure.  Assuming  that  the  unsteeidy 
wake  coincides  vnth  the  steady  one,  which  in  general  is  not  a plane  surface  behind  the  wing,  introduces  the 
problem  of  calculating  integrals  vnth  oscillating  integrands  over  curved  surfaces,  practically  excluding 
analytical  formulations.  The  alternative  ’’f  discretizing  the  wake  sheet  and  calculating  the  contribution  of 
the  individual  elements  ?s  propcsea  by  Morino  and  Kuo  (Refs.  10,  ll)  leads  to  the  question  of  truncation. 

Besides  the  effects  of  the  wake,  the  general  expression  (A1.6)  for  9 does  also  include  the  effects  due 
to  wing  thickness.  However,  in  this  context  it  is  important  to  realize  that  lifting  surface  theory  gives  al- 
ready a sufficiently  accurate  prediction  of  the  aerodynamic  forces  on  the  viing,  while  at  the  same  time  the 
wsdce  effects  are  incorporated. 

Taking  into  account  the  above  arguments  eind  coneidering  the  goal  of  the  method,  being  the  calculation 
of  pressure  distributions  on  stores,  of  which  the  wake  is  not  as  well  defined  as  that  of  the  vring,  the  method 
is  simplified.  The  body  is  assumed  to  produce  no  wake,  while  the  lifting  surfaces  are  taken  to  be  infinitely 
thin.  With  these  assumptions  the  expression  for  the  velocity  potential  due  to  the  body  reduces  to; 


9„  <x,y.z>  = - 


jf^n 


Cb>  e 


ik  M (x-£„) 

n 00'  B e 


(A1.8) 


The  velocity  potential  due  to  the  lifting  surfaces  can  now  be  given  in  terms  of  a pressure  dipole  as  is 
common  in  existing  lifting  surface  theories. 
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APPEIIDIX  2 

THE  BOUHDAKf  COHDITIOIB  OK  AH  OSCILUTIHG  SURFACE 

The  surface  of  an  oscillating  body  is  described  the  fbnctional 

S <r>  = 0 


(A2.1) 


The  vector  r,  which  determines  the  position  of  points  on  this  surface,  can  be  written  in  the  form: 


r = r<r,t>=r  +r  <rt> 
o’  o e ^ o 


(A2.2) 


where  "r  represents  th_^  stationary  reference  position  and  the  tine  dependent  anplitude  of  the  notion, 

which  a?so  depends  on"? . The  vector  7 satisfies  the  functional 
o o 

S^<7>  = 0 (A2. 

describing  the  surface  of  the  body  in  the  reference  position. 

The  boundary  condition  requiring  the  flow  to  be  tangential  to  the  moving  surface  is  given  by: 


^ - as  . -oc 
Dt  " at 


on  S < r > = 0 


(A2.4) 


Except,  perhaps  for  some  very  sinple  bodies  the  functionals  S <7  > and  < 7^  > are  veiy  difficult  to 
determine.  On  the  other  hand  the  individual  point  on  the  oscillating  surface  can  be  given  in  terms  of  their 
reference  position  and  the  tine  dependent  anplitude  r as  is  done  in  (A2.2).  Also  the  normal  in  these  points 
on  their  reference  position,  cam  be  calculated  without  any  problen.  Therefore,  in  the  formulation  of  the 
boundary  condition  on  the  oscillating  body  the  derivatives  of  S will  have  to  be  expressed  in  terms  of  the 
anplitude  7^  and  the  normal  7.  This  can  be  done  as  follows. 

The  components  of  the  vector  7 are 

= X <x^,y^,z^,t> 


y = y < vyo’^o’*^ 

Interchanging  dependent  and  independent  variables  results  in: 

^0  = *0  < > 

Yq  = y^  < x,y,2,t  > 

< x,y,2,t>  (A2.6) 

With  (A2.6)  a point  <x,y,z>  can  be  related  at  any  tine  to  its  reference  position  <x^,y^,z^>  . The  func- 
tional for  the  moving  sxrrfaco  becomesj 

<r^>  = <x^  <Xiy,z,t  > ,y^  < x,y,z,t  > ,2^  < x,y,z,t» 

• S < x,y,z,t  > = 0 (A2.7) 

Assuming  that  I r I « I r 1 the  following  expressions  for  the  derivatives  of  S can  be  obtained  from 
(A2.7) 

(^2.8) 


_ fas  _ as  as 

7S  =-v  s - ■rrVx  + ~v  y + -^V  2 

0 0 ax  o e ay  o e oz  'o  e 

L o o o 


with  Xg,  y^  and  z^  being  the  components  of  the  displacement  vector 

7 and  V = (ar~i 
e ''o  '3x  ’ ay  az  ' 

o o o 

The  disturbjince  potential  in  a point  on  the  moving  body  can  be 
in  the  correspondir.g  point  of  the  surface  in  the  reference  position* 

T <x,y,z>  = 9 * 


(A2.9) 


in  the  value  of  the  poten+ial 


(A2.10) 


Similarly  one  can  write  for  the  disturbance  velocities; 

7<p  =V<p  I + Cfg  .V  ) Vv  j + 


(A2.11) 


Substitution  of  (A2.8),  (A2.9)  and  (A'.ll)  in  (A2.4)  and  retaining  only  terms  of  first  and  second  order 
in  the  amplitude  and  its  derivatives  gives 


dr  as  as 

-VS  .^+(q+v<p  ) . (vs  --^Vx  - — 

»o  0 9t  ' ' 0 0 dx  '0  e 9y 

0 0 

0 


Ve  - S;  ^0%)  + ^0%  • [ I ] = ° (A2.12) 


The  symbol  I indicates  that  the  term  preceding  it  will  have  to  be  evaluated  on  the  body  in  the  reference 
position.  ° 
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Using  the  bamonic  expansion  (A1.2)  for  9,  a time  harmonic  amplitude  of  the  form 


^e  <=^o>yo**o«*  > “ >^e  <*o»yo''o>  ® 


and  neglecting  prod  jots  of  r^  and  9^  leads  to  the  following  two  conditions 
n » 0: 

f5_+V<P„  1 ) . n = 0 


(A2.13) 


n = 1: 

V?>i|  . n - n ) • (tVx^  + nVzp  + . j(?e'^^^''’o  1 | " ° (A2.15) 

*0  *0  L O'* 

where  { , m and  n are  the  components  of  the  normal  vector  on  the  body  in  the  reference  position 
V S 

. . — ''00 


defined  hy  n 
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SUMMART 


An  analysis  is  given  of  aerodynaaic  loads  neasured  with  an  oscillating  wind  tunnel  model  represent- 
ing a wing  with  a tip  tank  and  a reaovable  pylon  with  store.  Attention  is  paid  to  the  interference 
effects  OB  the  wing  load  and  to  the  pylon  store  load  in  low  and  high  subsonic  flow. 
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1 nraoDncncsr 

■ili'taiT’  aircmft  configurationB  involve  frequently  wing-attached  atorea  of  conaiderahle  aize. 
When  dealing  with  aeroelaatic  prohleae  for  auch  oonfigurationa  two  queationa  ariae,  naaelyt 

- what  are  the  unateady  aerodynaaie  Icada  on  the  atorea 

- what  are  the  interference  loada  induced  hy  the  atorea  on  the  wing. 

Alao  the  Xm  wan  faced  with  tbia  queation  when  it  got  involved  in  the  flutter  inveatigation  for  a 
fighter-type  aircraft  with  a nuaher  of  wing  atorea,  differing  in  aise  and  poaition.  Appropriate  general 
data  concerning  aerodynamic  loada  on  oacillating  wing-atore  coahinationa  waa  not  available,  ao  that  a 
reaearch  prograaae  was  atarted  at  the  XLR  involving  wind  tunnel  aeaaureaenta  on  a model  of  a represent- 
ative wini^atore  coahination  and  the  development  of  a method  to  calculate  aerodynamic  loads  on  oacilla- 
ting wing-atore  combinations  in  subsonic  flow.  The  latter  subject  ia  reported  on  in  ref.l.  In  this  way 
a compariacQ  of  measured  and  calculated  reaults  ia  possible  and,  because  of  the  fact  that  the  wind  tun- 
nel model  has  been  scaled  to  the  wing  of  the  aircraft  in  question,  the  results  can  be  used  as  input  to 
practical  aeroelaatic  calculations. 

In  this  paper  a brief  description  is  given  of  the  wind  tunnel  measurementa  and  the  analysed  results 
now  available  are  diacuaael. 

2 KODEL  BESCRIPTIOir  ATO  TEST  COHDITIOKS 

A picture  of  the  model  Installed  in  the  wind  t\mnel  is  presented  in  fig.l.  It  is  a semi-span  model 
equipped  with  a tip  tank  and  a pylon  mounted  store.  (The  full-scale  wing  can  also  carry  an  inboard  pylon, 
but  in  the  wind  tunnel  model  this  pylon  has  been  left  out  because  of  the  expected  timnel  wall  intei> 
ference).  Some  geometrical  data  concerning  the  wing  and  the  stores  are  given  in  fig.  2.  The  wing  has  a 
moderate  sweep  angle  and  ia  strongly  teqiered.  Wing  thickness  is  4,8  ^ and  the  profile,  being  constant 
along  the  span,  shows  a droop  nose. 

During  the  meaaxirements  the  model  was  able  to  perform  a pitching  oscillation  about  an  axis  perpen- 
dicular to  the  wall.  The  pitching  axis  could  be  shifted  to  locations  at  I3  ^ and  50  % of  the  root  chord. 
The  model  was  driven  at  a reaonance  frequency,  which  could  be  varied  by  changing  springs  in  the  support- 
ing structure  of  the  pitching  euds.  However,  the  pitching  frequency  was  always  well  below  the  frequen- 
cies at  which  elastic  deformations  of  the  mod<:i  could  be  observed.  The  vibration  mode  was  measured  by 
12  acceleration  pick-ups,  of  which  6 were  installed  in  the  wing,  4 in  the  pylon  store  (2  vertical  and 

2 lateral)  and  2 in  the  tip  tank.  The  vibration  mode  for  which  results  are  presented  in  this  paper, 
was  a nearly  pure  pitching  oscillation  ; lateral  motion  of  the  pylon  store  could  be  neglected. 

‘Po  gain  a detailed  pressure  distribution  over  the  model  surface,  the  model  was  provided  with  340 
pressure  taps  of  which  88  were  distributed  over  both  upper  and  lower  wing  surface,  86  over  the  pylon 
store  and  78  over  the  tip  tank.  The  tap  locations  on  the  wing  and  the  pylon  store  are  indicated  in 
fig. 3. 

The  pressure  measuring  technique  was  described  very  recently  in  ref.  2 in  an  application  at  high 
subsonic  speeds.  In  this  reference  it  was  demonstrated  that  the  accuracy  of  the  measuring  technique  is 
satisfactory.  The  essence  of  this  method  is  that  through  vinyl  tubes  .'.ll  pressure  taps  are  connected 
with  a small  number  of  scanning  valves  outside  the  test  section  of  the  wind  tunnel.  The  measured 
pressures  are  corrected  using  the  calibrated  transfer  functions  of  the  tubes.  In  the  present  wind 
tunnel  model  the  calibration  was  performed  by  means  of  a number  of  in  situ  micro  miniature  pressure 
transducers,  installed  in  one  of  the  measuring  sections.  To  obtain  unsteady  pressure  coefficients  the 
pressures  have  been  related  to  a displacement  s in  a reference  point  (x/c  » 0,871  ; y/s  - 0.143). 

Transition  of  the  boundary  layer  was  assured  by  carborundum  strips  on  the  wing  and  the  stores. 

Or.  the  wing  the  strips  were  located  at  20  ^ of  the  chord  on  the  upper  surface  and  at  40  ^ on  the 
lower  surface.  The  strips  on  the  tip  tank  and  pylon-store  were  located  at  8. 5 and  11  jC  of  their 
length,  respectively. 

The  model  tests  were  performed  in  the  large  treuisonic  wind  tunnel  HST.  Parameters  in  the  test 
conditions  were:  configuration  (pylon  store  and  pylon  were  removable),  pitching  axis  location,  fre- 
quency, )!ach  number,  stagnation  pressure  and  angle  of  attack  (only  in  the  steady  measurements;  in  the 
unsteady  measurements  wing  angle  of  attack  was  zero). 

An  analysis  of  the  influence  of  all  these  parameters  is  in  process,  A part  of  the  results  now 
available  is  discussed  in  this  paper.  They  involve  pressure  distributions  on  the  wing  and  the  pylon 
store  for  the  configuration  shown  in  fig.l,  a pitching  axis  at  1^  frequencies  of  0,11  and  19  Hz, 

Mach  muibers  of  O.49  and  0.8  and  a stagnation  pressure  of  1 atm.  In  the  steady  measurements  the  wing 
angle  of  attack  is  0°  and  3°  with  respect  to  the  wing  reference  plane. 

3 DISCUSSION  OP  RESULTS 

In  the  dlucussions  in  this  section  the  following  line  wiH  be  pursued.  First  steady  results  are 
p|^esented,  next  quasi-steady  refults  which  have  been  derived  from  the  steady  results  for  o » 0°  and 
3 and  finally  the  unsteady  results.  In  doing  so  attention  ca.i  be  paid  separately  to  the  interference 
effects  between  wing  and  pylon  store  which  are  essentially  of  the  same  kind  for  the  quasi-steady  and 
the  unsteady  results  and  to  the  influence  of  frequency. 

3.1  Results  for  the  wing 

3.1.1  Steady  and  quasi-steady  pressure  distributions 
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In  figure  4 some  steady  chordwise  pressure  distributions  are  presented,  which  are  typical  for 
this  wing  model  in  low  subsonic  flow  (X  <•  0.45)  at  zero  angle  of  attack.  The  distributions  on  both  the 
upper  and  lower  sxirfaoe  of  the  wing  are  given  for  the  configurations  with  and  without  the  pylon-store. 
On  the  lower  surface  near  the  leading  edge  they  show  a largo  suction  peak,  which  is  caused  by  the 
droop  nose  of  the  wing  profile  and  therefore  is  characteristics  for  all  steady  pressure  distributions 
shown. 

The  addition  of  the  pylon  and  store  under  the  wing  has  a large  effect  on  the  flow  characteristics 
around  the  wing  in  the  neighbourhood  of  the  pylon.  The  divei-sion  of  the  flow  due  to  the  store  causes 
the  local  velocity  under  the  wing  to  speed  up,  resulting  in  a considerable  change  in  the  pressure 
distribution  on  the  lower  surface.  On  the  upper  surface  this  interference  effect  is  an  order  of 
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■a^itud*  SMiller,  whila  it  ia  raatricted  aaialy  to  the  vicinity  of  the  pylon-atore  attachnent. 

Ihoreasini'  the  freeatreaa  Machnuaher  only  pronounoea  the  picture  described  above.  While  the  wing 
with  tip  tank  only  becoaes  srqjercritical  between  X • O.85  and  X > 0.9  on  the  lower  surface  near  the 
leading  edge,  the  presence  of  the  store  lowers  it  down  to  just  below  X - 0.8.  Figure  5 shows  a few 
chordwise  distributions  for  this  Xachnunber  at  sero  angle  of  attack.  Clearly  the  floir  over  the  lower 
surface  of  the  wing  has  already  becoae  supercritical  just  ln''oard  of  the  pylon  and  a weak  shock  night 
have  been  foiaed  already  in  that  region.  Figure  6 presents  soae  quasi-steady  chordwise  pressure  distri- 
butions for  low  subsonic  flow  (X  - 0.45)*  niese  quasi-steady  distributions  can  be  interpreted  as  the 
"unsteady"  airloads  for  infinitely  slow  pitching  oscillations. 

As  in  the  steady  distributions  the  influence  of  the  store  is  found  to  be  largest  in  the  vicinity 
of  the  pylon.  The  quasi-steady  distributions  show  interference  effects  on  both  the  tq>per  and  the  lower 
surfaces.  On  the  upper  surfue  a decrease  of  the  quasi-steady  pressure  is  found  and  on  the  lower  sur- 
face an  increase  with  the  largest  interference  on  the  lower  wing  surface  inboard  of  the  pylon.  Outboard 
of  the  pylon  the  pressure  change  due  to  the  store  addition  is  largest  on  the  iqiper  wing  surface. 

For  high  subsonic  flow  (X  - 0.8)  the  quasi-steady  distributions  are  the  ease  qualitatively.  The 
only  difference  is  the  appearance  of  soae  wiggles  caused  by  the  apparent  existence  of  a weak  shock  on 
the  lower  surface  near  the  pylon,  which  was  already  noticed  in  the  steady  distribution  (see  fig.7). 


3.1.2  Steady  tuid  quasi-steady  spanwise  load  distributions 


The  BteeuSy  and  quasi-nteady  spanwise  loa^  distribution,  the  latter  being  derived  from  the  measured 
steady  pressure  distributions  at  a • 0 and  3 , are  given  in  figure  8a  for  X - 0.43  and  0.8.  For  X • 0.8 
these  experimental  results  are  also  compared  with  theoretical  calculations  performed  with  the  KUt  panel 
method  (see  reference  3). 

From  the  figure  it  follows  that  for  all  flow  conditions  the  interference  between  wing,  pylon  and 
store  results  in  a considerable  decrease  in  the  steady  lift  along  the  span  of  the  wing,  which  is  maximal 
near  the  position  of  the  pylon.  Clearly  this  effect  corresponds  with  the  large  change  in  pressure  over 


the  lower  wing  surface  observed  in  the  steady  chordwise  pressure  distributions.  The  pylon  also  takes 
away  part  of  the  circulation  from  the  wing,  causing  at  that  loeatign  s.  jgip  in  the  lift  distribution. 
'Hiis  Jimp  changes  its  sign  if  the  angle  of  attack  increases  from  O'^  to  3 , indicating  u change  in  the 


direction  of  the  side  force  on  the  pylon  and  store.  AC 

'Bie  influence  of  the  pylon  and  store  on  the  quasi-steady  wing  load  r—  is  shown  in  figure  8b. 

ACj 

It  is  chaimcterised  by  an  increase  of  r—  inboard  ot  the  pylon  and  a decrease  on  the  outboard  side. 


This  can  be  traced  back  to  the  quasi-steady  pressure  distributions  of  figure  6,  which  show  a inter- 
ference effect  being  different  on  both  sides  of  the  pylon. 

Both  figures  8a  and  8b  show  that  the  theory  overestimates  the  lift  which  probably  is  caused  by  the 
neglect  of  the  viscous  effects.  However,  the  trend  is  predicted  rather  well. 


3.1.3  Unsteady  pressure  distsrlbutions 

Figures  9 through  11  show  a set  of  unsteady  chordwiee  pressure  distributions  which  are  representative 
for  the  model  oscillating  in  subsonic  flow.  In  figures  9 and  10  the  pressure  distribution  in  low  sub- 
sonic flow  (X  • 0.43)  on  the  iq>per  and  lower  wing  surface  are  given  for  frequencies  of  11  and  I9  Bs. 

As  far  as  the  real  part  of  the  unsteady  pressures  is  concerned,  the  influence  of  the  store  is  very 
similar  to  what  is  found  for  the  quasi-steady  distributlonB}  the  pressure  on  the  upper  surface  shows  a 
decrease,  while  on  the  lower  surface  an  increase  is  found.  The  effect  of  frequency  appears  to  be  very 
limited.  Hie  imaginary  part  of  the  pressure  distributions  remains  more  or  less  uneffected  for  both  fre- 
queuciee.  For  high  subsonic  flow  (K  >0.8)  as  shown  in  figure  11  the  picture  remains  the  same  quali- 
tatively, except  that  the  distributions  in  section  9 6 at  the  lower  surface  show  wiggles  which  are 

related  to  the  presence  of  a local  supercritical  flow. 


3.1.4  Unsteady  spanwise  load  distributions 

Hie  spanwise  unsteady  load  distribution  obtained  from  the  measured  pressure  distributions  are  pre- 
sented in  figure  12  for  both  low  and  high  subsonic  flow.  For  comparison  also  the  quasi-steady  distribu- 
tions have  been  Included. 

The  influence  of  the  pylon  and  store  is  ol early  very  similar  for  0,  11  and  I9  Hs.  For  the  frequency 
range  considered  the  jus^  in  the  circulation  across  the  pylon  tends  to  decrease  a little  with  Increasing 
frequency.  The  magnitude  of  the  interference  doss  not  exceed  the  I3  % of  the  local  lift  values  for  the 
wing  without  pylon  and  store.  The  total  lift  coefficient  which  is  obtained  through  spanwise  integration 
shows  an  interference  effect  of  less  than  3 The  imaginary  part  of  the  unsteady  King  load  is  hardly 
effected  by  the  pylon-store  addition. 

3.2  Results  for  the  store 
3.2.1  Pressure  distributions 


Pig. 13  shows  pressure  distributions  along  the  store  in  4 sections.  These  sections  are  defined  by  the 
angle  f.  Wing  angle  of  attack  is  zero.  On  the  left  side  the  results  are  given  for  X > 0.43,  the  right 
side  for  X > 0.8. 

Starting  the  discussion  with  X » 0.43  'the  first  observation  is  that  the  pressure  distributions  for 
different  f-angles  are  very  similar.  This  is  aue  to  the  fact  that  the  wing  angle  of  attack  is  sero,  so 
that  the  presnice  of  the  wing  is  hardly  felt.  The  store  has  a small  negative  angle  0^  attack,  which  is 
reflected  in  the  positive  pressures  at  the  nose  being  higher  for  f » 0°  than  for  180  . Two  minima  occur 
in  each  sect'*  on,  clearly  related  to  the  transition  of  the  varying  cross-section  of  the  nose  to  the 
cylindrical  part  and  to  the  transition  of  this  part  to  the  aft-body. 

Hie  observations  for  X > 0.8  are  the  same.  In  this  case  theoretical  values  have  been  included  cal- 
culated with  a panel  method  presented  in  roferenoee.  The  agreement  is  very  satisfactory.  At  the  end  of 
the  cylindrical  part  the  flow  becoaes  slightly  supercritical,  just  below  the  part  of  the 
wing  where  the  flow  is  also  supercritical.  Obviously  this  phenomenon  is  caused  by  the  geometry  of  wing- 
pylon-store,  which  admits  only  a narrow  passage  for  the  flow. 


In  fig.  14  quMi-ate«dy  and  unateady  raaolia  ar«  praaan-ted  for  -the  aua  4 aectlona  aa  in  fig.l3>  Mach 
nnabar  is  O.43.  Tba  qaani-ataadj  proaaore  diatributlona  and  the  real  parta  of  the  nnateadjr  presaura  diatri- 
'butiooa  ahow  asatotly  iha  aaaa  -tandenciaa,  Tna  Mgnitudaa  of  tha  iaaglnai7  parta  ara  ralativaly  vei7  as«ll. 

In  fig.  15  aaalosoua  raanlta  are  ahown  for  X • 0.8,  where  again  a co^Modaon  hatween  Maanrad  and 
calculated  raanlta  la  poaaihla.  The  agreeaent  ia  fairly  good  { the  diffaranoea  ara  largaat  orer  the  rear 
half  of  the  atore,  poaaihly  due  to  huundaiy  layer  and  separation  effects. 

3.2.2  load  diatrihntlone 

In  this  section  quasi-steady  and  unsteady  results  ara  presented  in  the  fori  of  C - and  C -diatri- 

•V  jTv 

hutions  along  the  store.  In  fig.  16  the  C -distrihation  ia  shown  for  the  X - 0.43  and  0.8.  The  reaarks 

as  for  fig.  I4  can  also  he  aode  here.  A q^litative  analysis  of  the  load  diatrlhution  aay  he  given  using 
the  result  of  slender  body  theory  that  the  load  disirihution  is  proportional  to  JL  (S(x)a(x))  , where 

S(x)  is  the  local  croso-sectional  area  and  s(x)  the  local  nomal  wash  as  functions  of  the  co-ordinate  x 
along  the  store.  The  nomal  wash  a(x)  is  a coahination  of  the  geoaetric  angle  of  attack  a and  the 
wing-induced  tqiwaah,  which  increasss  towards  the  wing  and  passes  on  to  a downwash  near  the  leading  edge. 

The  resulting  downwash  distrihutlon  along  the  store  thus  will  have  the  global  shape  as  sketched  in 
fig.l8.  The  derivative  of  the  product  S(x)o(x)  then  leads  to  a force  upwards  on  the  front  part  of  the 
store  and  downwardr  on  the  rear  part,  in  accordance  with  the  C -distributions  of  fig.l6. 

*b 

A final  reeiark  concerns  the  siaall  imaginary  parts  of  C in  fig.  16,  indicating  that  the  damping 

effect  of  the  store  loads  on  the  oscillatory  motion  is  negligible. 

Lastly,  fig.  17  presents  the  distribution  of  the  coefficient  C along  the  store  for  both  Kach  nua- 

^b 

bers.  Also  here  the  quasi-steady  and  the  real  parts  of  the  unsteady  results  are  almost  identical.  Again 
slender  body  theory  nay  be  used  for  a qualitative  explanation  of  the  load  distribution.  Due  to  the  pre- 
sence of  the  swept  wing  the  side  wash  along  the  store  increases  gradually  towards  the  tail  as  is  shown 
qualitatively  in  fig. 18.  The  resulting  side  load  ia  pointing  outwards  nearly  all  along  the  store.  It  is 
very  small  over  the  front  part,  shows  a maximum  near  the  wing  leading  edge  and  falls  off  towards  the  tail. 
Through  integration  of  the  aide  load  distribution  along  the  store  a side  force  coefficient  is  found.  The 
modulus  of  this  coefficient  ia  nearly  2 ^ of  the  modulus  of  the  lift  coefficient  if  both  are  made  dimen- 
sionless in  the  same  way. 


4 coNcuronfG  remarks 

In  this  paper  a discussion  is  given  of  the  aerodynamic  loads  on  an  oscillating  wing-store  combina- 
tion measured  at  subsonic  speeds  duidng  a wind  tunnel  experiment.  Qsphasis  is  put  upon  the  influence 
of  the  store  on  the  unsteady  airload  on  the  wing  and  on  the  loads  acting  on  the  store  itself. 

It  is  ahown  that  for  the  configuration  under  consideration  the  store  has  a limited  influence  on  the 
wing  loads.  The  maximum  effect  0.1  the  local  unsteady  lift  coefficient  amounts  to  about  I5  while  the 
maximum  effect  on  the  overall  unsteady  lift  coefficient  ia  about  5 

The  real  part  of  the  unsteady  load  on  the  store  Msembles  very  much  the  quasi-steady  load}  the 
imaginary  part  (damping)  appears  to  be  very  small.  Noteworthy  is  the  wing-induced  lateral  load  on  the 
store.  The  resulting  side  force  is  in  the  order  of  2 ^ of  the  unsteady  lift  force  on  the  wing. 
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Plff.3  Pooitlon  of  pressure  holes  on  the  ving  and  pylon  store. 
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Fig.4  Influence  of  pylon  and  store  on  the  steady  chordwisc  pressure  distributions 
on  the  wing  in  low  subsonic  flow. 
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Kig.5  Influence  of  pylon  and  store  on  the  steady  ehordwise  pressure  distributions 
on  the  wing  in  high  subsonie  flow. 
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Pig,  10a»  Influence  of  pylon  and  store  on  the  unstaeidy  chordwiee  pressure  distributions 
at  K • 0.45  and  f • 19  Rs  S wing  \:pper  aide. 


Pig, 10b  Influence  of  pylon  and  store  on  the  unsteady  chordwise  pressure  distributions 
at  N > 0.45  and  f > 19  Ha  s wing  lower  side. 
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Fig.  If)  Quasi -stendy  and  unstuady  pressure  distributions  along  the  store  at  high  subsonic  speed 
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Fig.  16  Quasi -steady  and  unsteady  normal  force  distribution  along  the  store 


Fig.  17  Quasi -steady  and  unsteady  side  force  distribution  along  the  store 


WING  WITH  STORES  FLUTTER  ON  VARIABLE 
SWEEP  WING  AIRCRAFT 


O.Sensburg,  A.Lotze  and  G.Haidl 
Messerschmitt-Bolkow-BIohm  GMBH 
Unternehmensbercich  Flugzeuge  ■ cntwicklung 
Ottobrunn  bci  Mfmchen 


1 . INTRODUCTION 

The  large  number  of  wing  movmted  stores  with  varying  mass  and  inertia  properties  for 
modern  fighter  concepts  especially  in  conjunction  with  a variable  wing  geometry  requires 
economical  procedures  for  investigation  of  total  airplane  dynamics*  There  are  two  basi- 
cally different  analytical  methods  which  are  commonly  used  to  explore  wing-store  flutter 
characteristics : 

A)  Determination  of  the  stiffness  matrix  of  the  coupled  system  by  application  of 
subtructure  techniques  and  calculation  of  normal  vibration  modes: 

This  method  has  the  advantage  that  every  degree  of  freedom  is  implemented  (no  trun- 
cation effects).  There  ore  two  major  disadvantages,  namely  that  for  total  airplane 
analyses  one  must  work  with  very  large  matrices  and  that  for  each  mass  variation 
leading  to  different  normal  modes  the  unsteady  aerodynamic  forces  must  be  produced 
anew.  For  the  latter  reason  the  method  is  restricted  to  the  application  of  two-di- 
mensional unsteady  aerodynamic  forces  with  all  known  Inaccuracies  at  high  surface 
sweep  angles.  Uce  of  correction  methods  to  the  unsteady  aerodynamic  forces  reduces 
the  error  margin  but  still  leaves  uncertainties  about  the  absolute  error  in  pre- 
diction of  flutter  speeds. 


Representation  of  the  total  aircraft  dynamics  by  superposition  of  clean  wing  canti- 
lever modes  and  cantilevered  external  store  modes. 

Here  wo  have  the  advantage  of  wing  branch  modes  and  unsteady  aerodynamic  forces 
not  changing  with  external  store  variation.  As  will  be  shown  later  convergence 
problems  concerning  the  number  of  necessary  wing  branch  modes  will  arise.  If  many 
higher  order  branch  modes  are  utilized  then  the  order  of  the  eigenvalue  problem 
increase  and  the  amount  of  work  which  is  needed  for  good  definition  of  aerodynamic 
forces  will  be  unacceptable. 


It  is  also  possible  to  use  wing  branch  modes  which  incorporate  the  external  mass  rigidly 
attached  to  a wing  point.  This  method  contains  the  disodvantages  of  methods  A and  D 
since  wing  mode  shapes  will  change  with  varying  external  stores  thus  rcquxring  recalcu- 
lation of  unsteady  aerodynamic  forces  and  the  convergence  problems  will  not  be  alleviated. 

Modified  branch  mode  techniques  are  found  to  bo  most  suitable  to  solve  the  problem.  By 
this  approach  the  frequencies  and  modeshapes  of  the  coupled  system  are  obtained  by  super- 
position of  a limited  number  of  clean  wing  normal  modes  and  the  so  colled  junction  modes, 
which  are  employed  in  order  to  improve  the  convergency  of  the  results  related  to  the 
number  of  cantilevered  bare  wing  normal  modes.  The  procedure,  using  modified  branch  mode 
techniques  for  the  total  airplane  representation  and  crosschecking  the  solution  of  ground 
vibration  and  flutter  calculations  with  total  airplane  model  test  results  is  considered 
to  be  the  most  efficient  approach  for  obtaining  reliable  flutter  results.  It  will  be  shown 
that  only  free-free  dynamically  scaled  total  aircraft  models  can  give  good  correlation 
when  tuning  effects  occur.  The  most  efficient  method  of  Investigating  flutter  behaviour  of 
asymmetrical  store  configuration  is  therefore  the  testing  of  this  above  mentioned  model. 


2,  ANALYTICAL  MODEL 


Structural  representation 

For  the  representation  of  total  airplane  dynamics  the  following  sets  of  finally  retained 
generalized  coordinates  are  considered  for  vibration  and  flutter  analyses: 


o Airplane  rigid  body  modes 

o Cantilevered  fuselage  normal  modes 
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Cantilevered  tailcron  normal  modes  and  talleron  attachment  modes 
Cantilevered  fin  normal  modes  and  fin  attachment  modes 
Cantilevered  ving  normal  modes  and  wing  attachment  modes 
Wing  - external  store  Junction  modes 

Cantilevered  external  store  normal  modes  and  pylon  attachment  modes 


During  early  design  stages  each  airplane  component  was  represented  by  a simple  beam 
structure  to  allow  simple  parameter  variation:-  For  initial  investigations  also  the 
taileron  and  the  fin  were  assumed  to  be  rigid.  After  refinement  of  the  design  the 
fuselage,  taileron  and  fin  idealisation  was  replaced  by  a finite  element  representa- 
tion, according  to  the  complexity  of  the  structure.  For  the  cantilevered  component 
modes  a distinction  is  made  between  primary  modes  for  a main  structure  os  for  the 
wing  and  secondary  modes  for  substructures  as  for  example  an  external  store.  Attach- 
ment modes,  which  ore  rigid  body  modes  for  components,  are  used  to  vary  the  connection 
stiffness  of  main  structures  with  substructures.  One  attachment  mode  for  each  external 
store  has  been  introduced  to  establish  a yaw  degree  of  freedom  corresponding  to  the 
flexibility  of  the  pylon  control  rod,  which  provides  constant  streamwise  direction  of 
the  pylon  for  all  wing  sweep  positrons.  It  should  be  mentioned  that  according  to  the 
single  point  pylon  attachment  with  free  motion  of  the  pylon  in  yaw  relatively  co  the 
wing,  the  wing  modes  do  not  contribute  to  the  pylon  yaw  displacements. 

Determination  of  Junction  modes  for  the  wing  - externa]  store  attachment 


Following  the  procedure,  which  has  been  developed  for  the  dynamical  coupling  of  sub- 
structures with  statically  undetermined  attachments  using  constraint  modes,  (Ref.  1; 
2;  3)  junction  modes  can  be  derived  from  the  stiffness  matrix  of  the  main  structure 
and  the  stiffness  matrix  of  the  attached  siibstructure  respectively 


kAA  k« 
I.IA  111 


where  the  superscripts  A and  I refer  to  the  freedoms  of  the  attachment  and  the  interior 
structure  respectively.  The  coordinates  6 are  physical  displacements  of  the  attachment 
points  and  interior  points.  F and  k describe  the  corresponding  force  vector  and  the 
stiffness  elements  respectively. 

For  the  determination  of  Junction  modes  the  forces  at  all  interior  points  arc  set  to 
zero  while  the  displacements  of  the  attachment  successively  obtain  unit  values. 

Using  equation  (l)  this  procedure  leads  to  a number  of  Junction  modes  corresponding  to 
the  number  of  attachment  freedoms. 


I = unity  matrix 

Each  column  of  refers  to  a Junction  mode 


- $3 

* , 


I 


defining  the  displacements  at  the  interior  structure  points  and  having  unit  displace- 
ment in  one  attachment  freedom. 


6-3 


For  the  case  in  consideration,  a wing  with  pylon  mounted  stores,  equation  (2)  con  be 
used  to  calculate  the  junction  modes  for  the  wing  structure.  Combined  with  all  modes, 
which  effect  displacements  on  the  wing  structure,  the  resulting  modal  matrix  produces 
a nondiagonal  generalized  mass  and  stiffness  matrix. 


A typical  example  of  Junction  modes  for  a wing  with  a store  connected  at  a wing  inboard 
station  in  three  attachment  freedoms  (roll,  pitch  and  z-translation)  is  shown  in  Fig.  1. 
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Fig.  1 Wing  Junction  Modes 


For  the  connected  substructure,  the  store  having  the  flexibility  of  the  pylon,  the 
Junction  modes  turn  out  to  be  rigid  body  modes,  if  the  attachment  between  the  main 
structure  and  the  substructure  is  statically  determined  and  the  unit  displacements  at 
the  attachment  points  do  not  produce  deformations  or.  the  substructure. 


Restricting  our  explanations  to  the  wing-external  store  combination,  the  following 
generalized  coordinates  hove  to  be  introduced  into  the  calculation: 


Rigid  wing  modes  (R,  Wg ) 


Primary  modes  (P) 


Wing  Junction  modes  (J) 


Rigid  store  modes  (R.  St) 


Secondary  modes  (S) 


Compatibility  conditions  for  tho  boundary 
compatibility  equation 


/ 1 fh* ' 
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3 ! IUST  I S . 
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9 ttor 

9=^ 


= 0 


(4) 


ensure  that  the  displacements  on  the  boundary  of*  the  substructure  match  those  of  the 
primary  structure  on  the  attachment  points.  Before  this  operation  can  be  performed, 
the  motions  of  the  store  at  the  attachment  points  are  to  be  expressed  in  terms  of  local 
wing  coordinates.  This  can  be  achieved  by  a suitable  transformation 


(j)^  = T,(b 

Tstwe  • 


A* 

SToee 


(5) 


where  represents  the  motion  of  the  store  at  the  attachment  related  to  local 

store  coordinates • 


The  compatibility  equation  (4)  is  used  to  elimii.ute  dependent  coordinates.  PnrtStionanc 
the  generalized  coordinates  according  to 


9 depend  , 
3 independ 


(6) 


pqnntion  (h)  can  be  rewritten 


[§dcp«n<i]  depend]  [4^indee*nd]{3’n<*b»nd|  ^ 


(7) 


^ 9 depend 
^independ 


"[$  depend  [^I'ndependil 

[11 


|9  indepond  | 


(8) 


The  transformation 


independ 


rnduces  the  generalized  coordinates  to  the  finally  retained  coordinates  q j j 
whore  T2  is  defined  by  equation  (8).  This  operation  also  relates  the  displacem?n?s  ’ 
of  the  substructure  to  the  generalized  coordinates  of  the  main  structure.  For  convenience, 
the  rigid  store  modes  with  unit  displacement  in  pitch,  roll  and  z-translation  at  tho 
boundary  can  bo  considered  os  dependent  modes.  For  total  airplane  calculations  a simi- 
lar treatment  has  to  be  performed  for  each  connected  component,  as  for  example  tho  rigil 
wing  modes  are  used  to  satisfy  tho  compatibility 'conditions  for  the  wing-fus' ^ ago  nttach- 
mont. 
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System  equation  for  vibration  and  flutter  analyaes 

Applying  the  transformation  matrix  T_, the  generalized  mass  and  stiffness  matrix  as 
used  for  vibration  and  flutter  analyses  of  the  complete  system  can  be  calculated  from 
initially  obtained  generalized  mass  and  stiffness  matrices  for  the  cantilevered  compo- 
nents t 
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Kqtnz 


'N 

(9) 

l[r.] 

(10) 

Kgenl 


where  the  indices  I and  II  denote  the  generalized  coefficients  of  the  main  structure 
find  the  secondary  structure. 

The  eigenvectors  and  frequencies  of  the  coupled  system  are  obtained  by  solving  the 
eigenvalue  problem 


[j^9«n]T^J{Sindep}  ' 0 


The  flutter  problem  is  described  by 


• ([^  9«n  + j f C )]  | 


(12) 


“ 0 


where  : m^,  oj  j. 

s 

A 

8 

V 

CO  ^ 


reference  mass  and  frequency 
semispan  of  roforenco  plane 
area  of  reference  plane 
structural  damping 
true  airspeed 
flutter  frequency 


a 


•If 


flutter  vector 


C 'gen.  C ' ' gen. 


real  and  imaginary  parts  of  non  dimensional  generalized 


airforces  defined  for  di scret  values  of 
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Using  the  Q-R  Algorlthnms  the  flutter  equation  is  solved  for  the  equivalent  amount 
of  structural  damping  g necessary  to  provide  harmonlcal  oscillation. 

The  generalized  airforces  which  are  to  be  introduced  into  equation  (12)  can  either 
be  calculated  directly  for  the  constrained  modes  defined  by  equation  (8)  or  be  copiputed 
for  the  initially  employed  unconstrained  modes  and  transformed  by 


Jtoo  “ [^z] 
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Unsteady  aerod>  .amic  forces 

3-dimensional  unsteady  aerodynamic  forces  are  calculated  According  to  the  theories  of 
('l),  (5),  (6),  (7).  Subsonic  airforces  including  wing- tailplane-f in  interference  effects 
and  supersonic  aerodynamic  forces  considering  wlng-tailplane  interferences  can  easily 
be  introduced  into  flutter  analyses  by  usage  of  branch  modes. 

Investigations  show  that  there  is  little  influence  on  flutter  behaviour  caused  by  un- 
steady aerodynamic  forces  on  external  stores  and  store-wing  interference  effects. 

For  comparison  with  flutter  model  test  results  in  this  paper  only  calculation  for  the 
subsonic  region  at  .Mach  number  0.2  arc  discussed. 


3.  SUBSONIC  WIND  TUNNEL  FLUTTER  MODEL 

A dynamically  scaled  free  flying  subsonic  wind  tunnel  flutter  model  was  bu  It  at  MDD 
and  tested  in  the  flutter  tunnel  of  the  Eidgenossisches  Flugzeugwerk  tn  Emmen,  Switzer- 
land. 

The  wings  an  ’ the  fuselage  of  this  model  are  represented  by  elastic  axis  lumped  mass 
systems.  The  fin  and  tailplane  struciure  is  simulated  with  a beam  laticc. 

A set  of  stores  with  varying  mess  and  inertia  properties  was  fabricated  and  tl.ts. 
stores  could  be  quickly  changed  on  the  model.  The  store  was  attached  to  the  wing  by 
a beam  representing  the  pylon  stiffness. 

The  flutter  speed  and  frequency  of  the  model  was  found  by  increasing  the  tunnel  speed 
up  to  the  flutter  speed  and  recording  tie  frequency  with  accelerometers. 


'l . RESULTS 


As  already  mentioned  there  is  a huge  variety  of  configurations,  arising  from  the  attach- 
ment of  stores-largely  differing  in  weight  and  inertia  - to  wings  of  variable  sweep  air- 
planes. 


Fig.  2 illustrates  a few  major  points.  It  shows  flutter  speed  as  a function  of  wing 
sweep  angle  for  a store  with  a constant  weight  and  variable  radius  of  gyration.  One  can 
see,  that  the  increments  of  flutter  speed  with  wing  sweep  change  with  store  radius  of 
gyration.  It  is  also  important  to  note,  that  the  minimum  flutter  speed  occuring  at  diffe- 
rent sweep  angles  changes  with  radius  of  gyration  ns  is  illustrated  by  curve  1 and  2 of 
Fig.  2.  It  is  also  possible  that  the  flutter  speed  Increment  with  wing  sweep  angle  might 
not  be  sufficient,  bearing  in  mind  that  the  flight  envelope  is  also  extend,.d  for  higher 
sweep  angles.  The  flutter  behaviour  may  look  completely  different  for  a store  having 
another  weight  than  that  in  Fig.  2. 
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FIG.  2 FLUTTER  SPEED  VERSUS  SWEEP  ANGLE  FOR  DIFFERENT  RADII  OF  GYRATION 


An  effective  method  is  therefore  necessary  to  cover  all  possible  cases.  A subsonic 
flutter  model  can  be  used  and  measured,  flutter  trends  must  be  correlated  with  analysis 
results.  Having  had  good  correlation  at  lov  subsonic  Mach  numbers,  the  whole  flight  en- 
velope with  varying  Mach  numbers  and  altitudes  can  bo  covered  analytically.  All  these 
analyses  can  be  performed  effectively  with  our  developed  branch  mode  system. 


Correlation  of  test  and  analysis  results 


In  order  to  prove  the  validity  of  a theoretical  method  it  must  be  tested  with  the  most 
severe  boundary  conditions.  This  is  the  case  when  a store  is  rigidly  attached  to  a wing. 
We  calculated  the  normal  modes  of  this  problem  applying  three  different  methods: 


A.  Exact  solution 

For  this  solution,  the  store  mass  properties  were  added  to  the  wing  mass  matrix  and 
the  eigenvalue  problem  was  solved. 


D.  Popular  wing  branch  mode  solution 

The  bare  wing  branch  modes  are  used  as  assumed  modes.  A generalized  mass  matlx  is 
formed  by  pre  and  post  multiplying  of  the  store  mass  matrix.  This  generalized  mass 
matrix  is  added  to  xhe  wing  generalized  mass  matrix.  The  eigenvalue  problem  is 
solved  and  the  solution  vectors  are  used  to  superimpose  the  clean  wing  modes,  thus 
forming  normal  modes  and  a diagonal  generalized  mass  matrix. 


C.  Combination  of  wing  branch  modes  and  Jiinctlon  modes. 


The  same  solution  method  as  for  B is  applied. 
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FIG,  3 TORSION  MODE  SHAPE  RESULTING  FROM  DIFFERENT  APPROACHES  FOR  A 
CANTILEVERED  WING  WITH  A RIGIDLY  MOUNTED  EXTERNAL  STORE 
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FIG.  4 TABLE  OF  FREQUENCIES  AND  GENERALIZED  MASSES 
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Results  of  these  calculations  are  shown  in  Fig.  3 as  deflections  hr,  bending  angles  ©a 
and  twist  angles  Ov  along  the  wing  elastic  axis  for  the  second  torsion  mode  shape. 
Solutions  with  6 branch  modes  and  3 branch  and  3 junction  modes  are  compared  with 
the  exact  solution.  Whereas  the  branch-mode  junction  mode  system  matches  the  exact 
solution  very  well,  there  is  bad  correlation  with  the  6 branch  mode  system  for  the 
twist  angle  which  is  very  important  because  it  forms  to  a larp  extent  the  genera- 
lized mass  of  this  mode.  From  Fig.  3 it  can  be  deduced  that  higher  wing  branch  modes 
arc  not  cs  well  suited  for  representing  a mode  shape  that  features  on  increasing 
twist  angle  up  to  the  store  attachment  point  and  has  a constant  angle  from  there  on, 
than  the  torsion  junction  mode  is,  which  has  exactly  those  characteristic  (see  Fig.  1). 
For  tL  other  normal  modes  correlation  of  the  three  methods  is  much  better  - 
ore  not  shown,  but  their  generalized  masses  and  frequencies  are  tabled  in  Fig.  -i.  This 
table  proves  again  that  only  the  branch  mode- junction  mode  solution  matches  the  exact 
solution  accurately  enough. 

In  this  table  it  is  also  shown  that  the  branch  modes- junction  mode  solution  converges 
rapidly.  Only  2 branch  modes  and  one  junction  mode  are  needed  to  get  a reasonable 
match. 
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FIG.  5 FLUTTER  SOLUTION  RESULTING  FROM  DIFFERENT  APPROACHES 


Results  of  flutter  calculations  with  the  two  different  approaches  are  presented  in 
Fig.  5.  The  classical  wing  bending  torsion  flutter  problem  occurs.  With  the  branch 
mooe  system  the  flutter  speed  is  unconservative  22  % for  structural  damping  g = 0 
and  25  % for  g = 2 %.  Again  the  branch  mode- junction  mode  and  the  exact  solution 
correlate  very  well. 

It  should  be  noted  that  we  tested  the  method  with  a rigorous  case  i.e.  rigid  store 
attachment  stiffness.  For  stores  elastically  attached  to  a wing  the  branch  mode  system 
should  inevitably  lead  to  better  results  because  for  the  zero  attachment  stiffness 
case  the  clean  wing  modes  arc  exact  solutions. 
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FIG.  6 FLUTTER  SPEED  VERSUS  C.G.  POSITION  OF  THE  OUTBOARD  STORE 


The  V-g  plot  Fig.  7 depicts  this  behaviour  at  O.5  c.g.  shift.  Bending  mode  and  store 
pitch  mode  are  producing  flutter  mode  1 and  bending  mode  and  store  yaw  mode  produce 
flutter  mode  2 of  Fig.  6. 


FIG.  8 FLUHER  SPEED  VERSUS  STORE  WEIGHT  FIG.  10  FLUTTER  SPEED  VERSUS  RADIUS  OF 


FOR  INBOARD  STORE  CONFIGURATIONS 


GYRATION  FOR  INBOARD  STORES 


Fig.  8 Illustrates,  that  there  is  little  change  of  flutter  speed  with  store  weight  for 
a certain  inboard  store  (constant  radius  of  gyration).  This  is  attributable  to  the 
fact,  that  the  inboard  store  weight  changes  the  wing  bending  frequency  very  little. 
Fig.  9 shows  the  V-g  plot  for  W = 0.886.  Here  it  con  be  seen  that  the  flutter 

^Rcf. 

mechanism  is  wing  bending  store  pitch  coupling. (Mode  2,  3) 
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FIG.  9 FLUTTER  SPEED  VERSUS  DAMPING  AND  FREQUENCY  FOR  TOTAL  AIRPLANE  WITH 
A STORE  ATTACHED  TO  THE  INBOARD  WING 


Fig.  10,  showing  flutter  speed  versus  Inboard  store  radius  of  gyration,  proves  the  well 
known  fact,  that  the  lowest  flutter  speeds  occur  when  the  store  pitch  frequency  matches 
the  wing  bending  frequency.  The  two  small  V-g  plots  in  Fig.  10  demonstrate  that  the 
flutter  cases  for  store  pitch  frequency  below  wing  bending  frequency  are  mild  whereas 
strong  flutter  occurs  for  the  opposite  frequency  ratio. 


Fig.  11  ehow8  • strong  dspsndency  of  the  flutter  spaed  of  two  stores  Inboard  and 
outboard  on  a wing  with  inboard  store  weight.  The  V-g  plot  for  a W = 0.886  is 

given  in  Fig.  12.  y n,.- 


FIG.  11  FLUHER  SPEED  VERSUS  INBOARD  STORE  WEIGHT 
FOR  VARYING  INBOARD  STORES  AND 
CONSTANT  OUTBOARD  STORE 


RG.  13  FLUTTER  SPEED  VERSUS  INBOARD  STORE 
RADIUS  OF  GYRATION  FOR  VARYING 
INBOARD  AND  CONST.  OUTBOARD  STORE 


In  Fig.  13  variation  of  flutter  speed  of  an  airplane  with  two  stores  mounted  on  each 
wing  as  a function  of  radius  of  gyration  of  lie  inboard  store  is  depicted. 
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FIG.  12  FLUTTER  SPEED  VERSUS  DAMPING  AND  FREQUENCY  FOR  TOTAL  AIRPLANE 
WITH  ONE  STORE  ON  THE  INDOATD  WING  AND  ONE  STORE  ON  THE  OUTBOARD 
WING 
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Influenc-  of  flutter  model  auspenalon  syatein 

Host  of  the  flutter  cases,  experienced  during  our  tests  were  symmetrical  cases,  but 
we  have  also  encountered  antisymmetrical  flutter  cases.  We  have  also  found,  that 
there  can  be  a strong  influence  of  fuselage  tuning  effects  upon  store  flutter  cases, 
fluttering  symmetrically.  For  these  two  reasons  we  recommend  to  use  a freely  sus- 
pended total  aircraft  wind  tunnel  model  for  wing-store  flutter  investigation.  In  the 
next  figures  the  effect  of  boundary  conditions  on  the  flutter  speeds  wing-store  con- 
figurations will  be  illustrated.  These  pictures  also  prove,  that  the  usage  of  uncorrec- 
ted two-dimensional  unsteady  air  forces  may  lead  to  large  errors  in  flutter  predictions. 


FIG.  14  FLUTTER  SPEED  VERSUS  PYLON  PITCH  STIFFNESS  FOR  FREE  AND  RESTRICTED 
AIRPLANE 


i- 


Fig.  shows  the  differences  in  flutter  speed  for  varying  boundary  conditions  as  a 
function  of  pylon  pitch  stiffness.  Curve  (2)  illustrates  how  the  flutter  speed  of  a 
cantilevered  wing  with  an  external  store  varies  with  pylon  pitch  stiffness.  The  im- 
portant store  pitch  mode  shape  is  presented  for  different  pylon  pitch  stiffness.  These 
mode  shapes  allow  a qualitative  assessment  of  flutter  behaviour.  The  mode  shape  for 
100  % pylon  pitch  stiffness  shows  a nodal  line  at  obout  throe  quarter  wing  chord.  Both 
nodal  lines  for  50  % end  200  % pylon  pitch  stiffness  are  further  forword  on  the  wing 
chord,  thus  lending  to  higher  flutter  speeds.  The  free-free  aircraft  flutter  model  has 
a different  flutter  behaviour  versus  pylon  pitch  stiffness  (Curve  l).  This  time  the 
nodal  line  is  forward  for  100  % pylon  pitch  stiffness.  When  the  heave  mode  of  the  total 
airplane  model  is  restricted  (Curve  3),  than  the  model  behaves  as  the  cantilevered  wing. 
From  Fig.  I'l  it  can  be  learned,  thot  the  flutter  speed  can  be  largely  underestimated 
with  cantilevered  wing  models  for  certain  pylon  pitch  stiffness. 
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FIG.  15 


FLUTTER  SPEED  VERSUS  DAMPING  FOR  FREE  AIR''LANE  AND  CANTILEVERED 
WING  FOR  TWO  DIFFERENT  PYLON  PITCH  STIFFNESSES 
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V-g  plots  for  different  boundary  conditions  and  two  pylon  itch  stiffnesses  ore  given 
in  Fig.  15.  The  large  differences  in  flutter  speed  for  cantilevered  wing  and  total  air- 
plane with  two-dimensional  and  three-dimensional  air  forces  are  striking  for  the  case 
with  100  % pylon  pitch  stiffness.  Much  smaller  differences  ore  apparent  for  200  % pylon 
pitch  stiffness. 


Flutter  of  asymmetrical  external  store  configurations 

I This  phenomenon  can  only  bo  investigated  with  total  airplane  flutter  models.  Asymmetri- 

' cal  store  configurations  ore  possible  for  an  airplane  carrying  bombs  and  missiles. 

Fig.  l6  shows,  that  it  is  possible  to  have  increases  or  decreases  of  flutter  speeds 
with  sweep  angles  for  asymmetrical  store  configurations  compared  with  symmetrical  store 
configuration.  Whereas  Fig.  l6  a shows  higher  flutter  speeds  for  a store  with  a high 
; radius  of  gyration  on  one  wing  than  for  the  same  store  on  two  wings,  Fig.  t6  c denon- 

• strates  the  opposite.  A different  behaviour  occurs  for  two  stores  on  the  wing  (Fig. 

;•  l6  d,  l6  e,  16  f). 

Couplings  of  three  different  vibration  modes  ore  feasible.  These  vibration  mode 
' shapes  are  shown  in  Fig.  17.  In  Fig.  17  a first  wing  bending,  in  Fig.  17  b store 

f pitch  and  in  Fig.  17  c store  roll  is  represented.  The  mode  shapes,  considering  only 

the  right  side  wing  change  little  for  symmetrical  and  asymmetrical  configurations, 
but  the  frequencies  vary. 
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An  explanation  of  the  flutter  mechanism  will  be  given  using  the  following  figures 
for  the  inboard  store  configuration.  In  Fig.  l8  a the  flutter  speeds  for  a symmetri- 
cal and  asymmetrical  store  configuration  are  drawn  versus  radius  of  gyration  of  this 
store  for  a constant  swoop  angle  .A..  = 25°  (they  are  taken  from  Fig.  l6  a,  l6  b, 

16  c). 

Fig.  l8  b presents  the  frequencies  of  the  different  modes  varying  with  radius  of 
gyration. 


one  can  see,  that  the  asymmetrical  store  pitch  frequency  is  below  the 
asymm^rlcal  wing  bending  frequency.  This  leads  to  a milder  flutter  case  than  that 
of  the  symmetrical  configuration  where  store  pitch  frequency  is  above  wing  bending 
frequency.  For  this  reason  we  ^ot  a higher  flutter  speed  for  the  asymmetrical  con- 
° - - - =0.9  the  wing  bending-store  pitch  flutter  is 


figuration.  For;^j^  = 0.7  to 
detuned  and  a mode  coupling  stor 
pitch  and  store  roll  mode  have  the  same  frequency  at 
speed  for  this  coupling  occurs. 


= 0.6  the  lowest  flutter 


The  different  behaviour  of  the  si-rametrlcal  and  a asymmetrical  configurations  may  be 
explained  by  the  different  mode  frequencies  of  these  conf igtu-atlons . 
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A PARAXSrai  J STUDY  0?  WIKS  STORE  FLUTra? 
by 

L.  CH2STA 

A2RITAUA  S.p.A.  - TOHIHO 


INTRODUCTIOll 

At  Aerxtali%  an  intensive  study  progran  based  on  tests  and  calculations  on  a sweepable  vdng  nodel 
has  been  carried  out  in  the  past  years  anc  is  still  in  progjv'ss, 

NiJst  of  the  work  has  been  devoted  to  the  stores  problem  to  study  the  influence  of  different  paramet- 
ers in  order  to  optimize  the  flutter  performance  of  a wing  with  stores* 

Sore  than  3000  configurations  have  been  tested  in  che  wind  tunnel  and  many  computations  have  been  do 
ne  to  extend  the  research  to  oarancters  ttiatware  difficult  to  simulate  on  the  models. 

Fig,  1 shows  the  geometry  of  one  of  the  models  used. 

The  model  has  been  built  in  the  classical  way  with  a spar  simulating  the  bending  and  torsional  stiff 
nesses  and  nine  sections  simulating  the  inertia  of  the  wia~. 

The  pivot  and  fuselage  flexibilities  have  been  simulated  by  a bean  clamped  to  a ri/-id  rig  fixed  to 
the  tunnel  roof. 

It  was  possible  to  change  the  sweep  -ingle  by  a screw  with  a spring  incorporatec  in  order  to  simulate 
the  rigid  wing  fore  ana  ait  mode. 

Two  pylon  stations  have  been  proviued  at  suitable  positions  and  the  two  pylons  have  been  simulated 
by  beans  of  different  pitch,  lateral  and  yaw  stiffnesses. 

The  store  C.O.  position,  weight  and  moment  of  inertia  were  si.mulated  by  two  variable  masses  mounted 
on  a beam  at  seven  standard  positions  without  any  aerodynamic  shape  simulation. 


TESTS  PROCEDU-HES 

A complete  investigation  has  been  perfor-med  on  the  following  parameters  : 

- Store  mass 

- Store  radius  of  inertia 

- Store  C.a.  longitudinal  position  referred  to  the  pylon  attachment  point 

- Pylon  pitch  stiffness 

- Wing  sweep  angli 

Moreover  a limited  analysis  has  been  done  to  evaluate  the  influence  of  parameters  such  as;  Pylon 
length,  pylon  lateral  stiffness,  fuel  in  the  wing,  fore  and  aft  rigid  wing  stiffness,  open  flaps 
and  slats,  pivot  backlash,  stores  aerodynamic  shape  and  the  size  of  the  stores  vanes. 


RESULTS 

tost  of  the  studied  configurations  showed  a flutter  induced  by  the  coupling  of  the  wing  fundamen- 
tal bending  and  the  store  pitch  modes. 

Pig.  2 left  shows  a typical  contour  diagram  obtaineu  for  a configuration  with  the  inboard  station 
loaded  and  a wing  sweepangle  of  25  degrees. 

Equal  speed  lines  are  drown  with  reference  to  the  max  wind  tunnel  speed  in  the  store  weight  -radius 
of  gyration,  plane. 

On  the  right  side  of  the  figure  a section  of  the  diagram  for  fixed  weight  is  shown  together  with 
the  frequency  behaviour  of  the  relevant  flutter  modes. 

It  can  be  readily  seen  that  the  flutter  speed  decreases  increasing  ti."  store  radius  of  inertia  un- 
til the  frequency  of  the  store  pitch  mode  is  higher  than  the  fundamental  lending  and  starts  again  to 
increase  as  soon  as  the  pitch  frequency  becomes  lower  than  the  fundamental. 

This  second  type  of  flutter  is  normally  very  mild;  it  is  very  difficult  during  ■'■he  test  to  define  ccr 

rectly  the  flutter  speed,  and  consequently  to  make  a correct  comparison  with  the  computation  without 
a careful  evaluation  of  the  model  damping. 

On  Figure  3 ■the  results  obtained  for  the  basic  configurations  with  only  one  store  at  the  inboard  p^ 
Ion  are  shown. 

The  figure  is  a oompxisition  of  many  diagrams  like  that  one  shown  on  Fig.  2. 

The  diagrcims  are  set  from  left  to  right,  in  order  of  increasing  values  of  the  C.G.  forward  snift  and 

from  top  to  bottom,  in  order  of  increasing  swee  angle. 

It  is  clear  that  for  all  the  sweep  angles  the  forward  C.G.  shifting  produces  a sligit  reduction  of 
the  value  of  the  minimum  flutter  speed  associated  with  a reduction  of  the  area  encircled  by  a defin- 
ed constant  speed  line.  Such  area  increasing  the  C.G.  shift  is  progressively  reduced  to  the  region 
of  stores  with  small  weight  and  large  radius  of  inertia. 


An  increase  in  tlis  sweep  angle  produces  a rise  of  the  value  of  the  minicuiii  flutter  speed  and 
the  reduction  of  the  effectiveness  of  the  C.Q.  shift. 

Pig.  4 shows  a siailar  diagraa  obtained  for  a configuration  with  a pylon  pitch  stiffness  four 
tioe  the  one  used  in  fig.j 

Comparing  the  two  figures  it  is  possible  to  see  that  an  increase  of  the  pylon  pitch  stiffness 
does  not  change  the  value  of  the  Bininia  flutter  speed  but  shifts  it  towards  stores  of  high- 
er weight  and  radius  of  inertia  and  wduces  the  elfectiveness  of  the  forward  shifting  of  the 
C.G.  - 

On  Fig.  5 6 the  effects  of  the  same  parameters  for  a configuration  with  a store  at  the 

outboard  station  are  presented.  The  trend  is  the  same  ac  before,  only  the  effect  of  the  pylon 
pitch  stiffness  IS  loss  important,  because  the  total  deformability  of  the  wing  pylon  system 
from  the  root  to  the  store  attachment  does  not  change  toomuch  multiplyirg  the  pylon  nominal 
pitch  stiffness  by  four. 

Pig.  7 and  8 show  the  effect  on  the  flutter  speed  of  a fixed  and  inertially  well  defined 
store  at  the  outboard  pylon  for  the  same  wing  sweep  angles  and  inboarc  store  paraveters  pre 
viously  analyzed. 

The  configuration  chosen  for  the  outboard  store  was  inherently  very  stable  and  therefore  the 
result  was  only  a shifting  of  the  minimum  flutter  speed  toward  the  region  of  the  large  store 
weights  and  radii  of  inertia. 

Tne  diagram  presented  before  could  be  considered  in  many  aspects  ideal  because  only  one 
type  of  flutter  occurred  in  the  configurations  chosen. 

Some  different  types  of  flutter  can  appear  changing  other  parameters. 

The  lateral  stiffness  of  the  pylons  comes  out  as  a very  important  parameter  for  heavy  sto- 
res with  a small  radius  of  inertia  like  bombs  mounted  on  double  or  triple  carriers. 

In  these  cases  a flutter  due  to  the  coupling  between  the  lateral  pylon  bending  and  the  pylon 
pitch  modes  can  appear. 

Fig.  9 shows  nine  coatoiir  plots,  obtained  for  the  same  wing  configuration  and  for  different 
pylon  stiffnesses  and  J.’i,  position,  in  which  the  distorsions  of  the  equal  speed  lines  due  to 
the  presence  of  different  flutter  modes  are  eviaont. 

On  top  of  fig.  10  three  sections  of  these  diagrams  for  a constant  store  weight  are  presented 
together  with  the  f^c^uency  trends  of  the  most  interest; ng  modes. 

Pylon  Cl  presents  the  normal  behaviour  because  the  lateral  pylon  mode  has  alway  a frequency  hi- 
gher than  the  pitch  mode. 

Pylon  B01  has  the  same  pitch  stiffness  but  a lower  lateral  stiffness  so  that  the  pitch  and  la 
teral  mode  frequencies  are  very  close  for  the  first  two  test  point  and  a flutter  due  to  the  ecu 
pling  of  these  modes  appears  at  speed  lower  than  espected. 

Pylon  30j2  has  about  tne  same  lateral  stifiness  as  ?yion  d01  but  a larger  pitcn  stiffness 
and  also  in  tnis  case  when  the  frequency  of  the  lateral  bending  mode  is  lower  or  close  to 
tr.at  one  of  the  pitch  mode  the  same  type  of  flutter  appears  at  speed  lower  than  expected. 

Also  when  bot!.  sta  ions  are  loaded  the  lateral  bending  stiffness  of  the  inboard  pylon  has 
been  found  import. int;  Fig.  11  show.,  a configuration  in  which  some  cases  of  flutter  cue  to 
the  coupling  of  the  lateral  bending  of  the  inboard  pylon  with  the  pitch  mode  of  the  outboard 
pylon  are  present. 

In  this  conf miration  the  outboard  pylon  yaw  mone  some  ti.'.ies  has  nearly  the  same  frequency  as 
the  flutter  mode  and  therefore  a rotating  motion  of  the  outboard  pylon  nay  appear  during  the 
flutter. 

The  most  peculiar  flutter  found  .s  however  that  one  origina*  d by  the  rigid  for'*  and  aft  mo 
tion  of  the  wing. 

This  Mode  couples  itself  very  stron-gly  with  the  wing  .orsion  when  a store  is  fitted  on  the  p^ 
Ion  because  of  the  large  vertical  distance  between  the  wing  and  the  store  C.C.  - Therefore,  if 
the  frequency  of  the  uncoupled  fore  and  aft  mode  .s  low  enough,  a new  mode  at  a lower  frequen- 
cy with  a large  amount  of  torsion  may  appear  inoucing  a flutter  at  a lower  speed. 

On  Pig.  12  bottom  the  frequency  trends  with  the  radius  of  inertia  of  the  relevant  modes  for 

two  cases  with  and  witnout  fore  and  aft  motion  are  presented  together  with  the  cc  responding 
flutter  trends  on  the  top  of  the  figare. 

It  IS  clear  that  for  heavy  stores  with  a small  radius  of  inertia  the  wing  fore  and  aft  mode  has 
a frequency  lower  than  the  store  pitch  and  couple..  it.^clf  with  the  fundamental  wing  bending  at 
a lower  speed.  This  is  made  evident  on  figs.  13  cuid  14  where  the  modes  nodal  lines  together  with 
the  V-g  plots  for  two  typical  cases  are  shovm. 

Apparently  the  aerodynamic  shape  is  not  .a  very  important  parameter,  only  if  the  store  is  ve- 
ry large  and  lignt  a small  i.screase  in  the  flutter  speed  has  been  found  (see  Fig.  15). 

Also  the  vanes  nay  increase  slightly  the  flutter  speec  when  the  flutter  is  mild  as  shown  in  fig.l6. 
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CONCUUSKKJS 

The  dia^ams  shovn  are  only  a few  santple  of  the  results  obtained  in.  the  worlc;  unfortunate- 
ly the  short  tioe  available  allows  only  the  presentation  of  some  highlights  and  their  discus 
Sion  in  a very  arid  way. 

It  is  anyway  possible  to  sunmarize  here  some  iriportant  general  conclusions  reached  during  the 
progress  of  the  study. 

If  an  aircraft  carrying  many  underwing  store  configurations  must  be  developed  a complete 
flutter  analysis  should  be  done  taking  into  accoui;t  the  effect  of  the  variation  of  as  many  pa 
rameters  as  possible  . 

In  fact,  with  an  aialysis  limited  to  the  icTiircd  storci  oonfigurations  or,  much  worse,  only  to 
the  configurations  considered  critical,  it  is  impossible  to  extrapolate  the  value  of  the 
flutter  speed  even  for  very  similar  configurations  and  to  make  any  forecast  of  the  critical 
speed  trends  if  some  parameter  is  changed. 

The  evaluation  of  the  effect  of  a change  of  the  value  of  some  parameters  is  the  .lormal  work, 
during  the  development  of  an  aircraft,  bit  this  is  very  cumbersome  and  difficult  to  do  if  cany 
underwing  store  configurations  have  to  be  cleared  together  . In  fact  ctn  increase  of  the  value 
of  one  parameter,  for  example  the  pitch  pylon  stiffness,  may  be  beneficial  for  one  store  and 
simultaneously  unfavorable  for  another  store  only  slightly  differing  in  the  i..ertia  characteris 
tics. 

To  be  able  to  judge  about  the  convenience  of  a modification  of  a parameter  it  is  therefore  ne- 
cecsary  to  have  well  i.i  mind  how  the  <^verall  picture,  and  not  only  the  single  points,  changes 
with  that  parameter. 

This  IS  even  more  complicate  for  an  aircraft  with  a variable  sweep  wing  where  it  is  necessaury 
to  get  a good  compromise  with  the  different  flight  envelope  requirements  at  each  sweep  angle. 

In  the  design  of  the  wind  tunnel  models  and  in  the  definition  of  the  mathematical  models 
for  tne  computations  it  is  important  to  Lake  into  account  all  the  degrees  of  freedom,  also 
those  normally  neglected  like  thi  in-plane  flexibility  of  the  wing,  the  lateral  and  yaw  flexi- 
bilities of  the  pylons. 

In  particular  for  the  computations,  it  is  necessaury  to  obtain  a good  accuracy  in  the  represen- 
tation of  the  wing  modes.  In  fact  most  of  the  discrepancies  found  in  the  cross  checks  bet- 
ween tests  and  co  iputations  were  due  to  inaccuracies  in  the  modes  definition  specially  in  the 
earlier  calculations  when  the  method  of  the  branch  modes,  without  junction  modes  was  used. 

For  this  reason  all  the  flutter  computations  have  been  done  using  wing-pylons  normal  modes  com- 
puted for  each  specific  configuration  teaching  in  this  way  a good  agreement  between  test  and 
calculation,  notwithstanding  some  further  difficulty  in  the  correlation  due  to  a str..,tural 
damping  of  ■’.he  model.  In  fact  if  the  flutter  is  mild,  and  normally  it  is,  a small  amotmt  of 
damping  is  sufficient  to  increase  very  much  the  value  of  the  flu"^ter  speed;  moreover  the  defi- 
nition of  the  flutter  speed  in  the  wind  tunnel  may  get,  in  these  cases,  inaccurate  since  it 
is  difficult  to  define  the  point  of  the  flutter  onset.  Anyway  a good  comparison  has  always  been 
found  introducing  the  measured  value  of  the  flutter  m de  damping. 

As  a last  sentence  it  has  been  found  very  convenient  and  powerful  to  use  a mixed  wind  tun- 
nel - oompat.tjon  flutter  analysis.  In  fact,  as  both  methods  b'.ve  been  proven  equivalent  in 
tne  accuracy  of  the  results  it  has  bean  possible  in  traa  wa'  to  reach  an  optimum  between  the 
time  necessary  to  got  the  results  and  the  difficulties  tc  represent  some  parameter  either  on 
wind  tunnel  or  mathe.matical  models. 
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IPL 
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Inboard  pylon  lateral  mode 

IPP 
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Inboard  pylon  pi .oh  mode 
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Inboard  pylon  yaw  mode 
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m 

Outboard  Pylon  lateral  mode 

OPP 

m 

Outboard  pylon  pitch  "ode 

opy 

m 

Outboard  pylon  yaw  mode 

V.-nax 

K 

I<ax  wind  t\umel  speed 

Vp 

m 

Flutter  speed 

VBII 

m 

Wing  second  bending 

WFA 

s 

Wing  fore  and  aft 

Wf'B 

ft 

Wing  fundamental  bending 

"ISa 

« 

Inboard  store  weight  (mouel) 

X 

«s 

Sweep  angle 

0).,, 

fit 

Flutter  frequency 

ft 

Model  frequencies 

Fig.  7 - Influence  of  fore  C.G.  position  and  wing  sweep  angle 

- Store  at  inboaid  and  outboard  station  outboard  store 
C.G.  position  .20  Cj.  fwd 

- Nominal  pylon  pitch  stiffness 


Fig.  8 - Influence  of  store  C.G.  position  and  wing  sweep  angle 

- Store  at  inboard  and  outboard  stations  outboard  store 
C.G.  position  .20  Cf  fwd 

- 4X  Nominal  pylon  pitch  stiffness 


RECENT  OBSERVATIONS  ON  EXTERNAL-STORE  FLUTTER 


Eugene  F.  Baird  and  William  B.  Clark 
Grumman  Aerospace  Corporation 
Bethpage,  New  York  11714 


INTRODUCTION 

The  problem  of  wing  flutter  with  external  stores  is  generally  regarded  as  one  of  the  most  complex 
flutter-prevention  considerations.  Tne  sheer  magnitude  of  the  number  of  configurations  to  be  studied. 
Figure  1,  all  too  often  staggers  the  imagination  of  the  flutter  engineer;  it  places  him  in  a state  of  near- 
exhaustion, contemplating  the  number  of  analyses  he  must  perform  to  fully  understand  the  problem. 

Often,  one  of  the  main  dilemmas  is  jost  where  to  start.  This  paper  attempts  to  show  how  ttis  problem 
can  be  approached  realistically  and  h<-w  its  magnitude  can  be  reduced  to  manageable  proportions.  The 
paper  first  discusses  the  current  app.'oach  tj  external-store  flutter,  and  then  touches  on  some  of  the 
recent  developments  being  made  to  better  cope  with  this  complex  problem. 

As  is  well  known,  classical  wing  flu  .ter  involves  the  coupling  of  at  least  two  iegrees  of  freedom,  one 
or  more  of  which  must  contain  motion  of  ,he  wing  hi  angle  of  attack.  Wing  external-  store  flutter  is  no  ex- 
ception to  this  rule.  The  basic  mech<.nism  of  external-store  flutter  can  be  thou^t  of  as  the  coupling  of 
wing  modes  with  store  modes  which  p<-oduce  wing  angle -of -attack  motion.  Most  external-store  flutter 
problems  occur  either  as  an  existing  b^tre-wing  flutter  mode,  modified  by  the  addition  of  store  mass  and 
inertia  on  a relatively  stiff  pylon,  or  as  a new  flutter  mode  introduced  by  the  addition  of  the  store  on  a 
relatively  flexible  pylon.  The  occurrence  and  characteristics  of  either  type  of  flutter  are  very  much  a 
function  oi  an  craft  configuration,  store  location  relative  to  the  aircraft,  and  system  mass  and  elastic 
properties.  It  is,  obviously,  impractical  to  provide  detailed  desigri  charts  for  any  and  all  possible  com- 
bmations  of  these  parameters.  This  paper  will  show  how  the  problem  can  be  approached  and  interpreted 
in  a manner  sufficient  to  gain  the  understanding  required  for  a realistic  problem  solution. 

EARLY  TRADEOFF  DESIGN  CONSIDERATIONS 

In  the  early  design  stages  of  an  aircraft,  the  flatter  analyst  is  generally  asked  to  provide  configura- 
tion tradeoff  information  necessary  to  arrive  at  an  optimum  ovoiall  design.  Analyses  are  usually  based 
on  very  preliminary  desigu  data  and  little  or  no  test  information.  This  problem  is  greatly  complicated 
when  the  aircraft  being  designed  is  required  to  carry  a large  variety  of  wi.ig-mounted  external  stores. 

If  the  design  happens,  for  example,  to  be  a strike  fighter,  then  the  tradeoffs  between  the  fighter  config- 
uration, which  is  required  to  carry  a few  air-to-air  missiles,  and  the  attack  version,  which  most  likely 
must  carry  a la;  ge  assortment  of  heavier  air-lo-ground  store.<?.  become  a major  design  consideration. 

How  much  weigl.t  penalty  should  be  imposed  on  the  fighter  version  to  have  a satisfactory'  attack  capability 
Early  design  flutter  studies  of  a reasonable  number  of  store  configurations  with  a strength -design  wing 
will  indicate  quickly  the  scope  of  the  problem.  It  is  not  unreasonable  to  find  a ratio  in  flutter  speeds  of 
two  Oi  more  between  the  bare  wing  and  the  wing  heavily  loaded  with  stores.  Such  a ratio  would  surely 
require  that  the  torsional  stiffness  of  the  wing  be  increased.  The  additional  stiffne  required  will  then 
depend  on  exactly  what  pe.  formance  is  needed  for  the  attack  mission.  An  early  check  on  the  aircraft 
performance  with  the  added  store  drag  is  certainly  a calculation  which  should  be  made.  These  results 
help  the  analyst  decide  just  what  level  of  stiffness  is  needed  to  satisfy  the  various  mission  profiles. 

The  most  gratifying  approach  to  carrying  external  stores  to  come  along  in  a long  while  is  the  con- 
cept of  conformal  weapon  carriage.  Figure  2.  This  concept  provides  for  the  mounting  of  stores  within 
a pallet  that  fits  along  the  bottom  of  the  aircraft  fuselage.  The  stores  are  submerged  and  held  in  place 
by  ejectc’’  racks  within  the  pallet.  Flight  tests  have  demonstrated  substantial  performance  improvements 
with  the  contormal  carriage  pallet.  Results  indicate  a very  dramatic  reduction  in  the  drag  associated 
with  external  stores.  However,  because  ihe  trends  in  weapon  dev.elopment  are  toward  very  large, 
heavy,  sophisticated  bombs,  the  use  of  the  pallet  -will  be  limited.  Because  sophisticated  weapons  have 
aerodynamic  fins  for  controllability  and  will  be  elongated  with  complex  guidance  units,  the  installation  of 
these  bombs  into  the  pallet  may  not  be  possible.  Therefore,  satisfactory  operation  of  sophisticated 
bombs  will  still  require  the  use  of  wing  pylon  mounting. 

Unfortunately,  there  is  little  hope  of  ever  having  enou^  quantitive  data  to  set  down  a complete  set 
of  guidelines  for  the  optimum  arrangement  of  external  stores  on  a wing.  Some  obvious  general  guide- 
lines to  follow  in  the  design  of  a new  aircraft  which  carry  many  external -store  configurations  are: 

• Keep  all  store  stations  as  close  to  the  fuselage  as  possible,  as  shown  in  Figure  3 

• If  store  stations  outboard  of  about  the  wing  mid-span  are  required,  minimize  the  possible  mass  ' 
inertia  variations  at  these  stations 

• Obtain  as  high  a bare-wing  flutter  speed  as  is  practical,  an  example  of  which  is  shown  in 
Figure  4 
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Once  the  configuration  has  been  established,  some  preliminary  analyses  must  be  performed  to 
establish  what  the  critical  flutter  modes  are  and  an  idea  of  what  the  critical  parameters  are.  The 
amount  of  work  necessary  in  this  preliminary  study  will  depend  to  a large  e^dent  on  the  corifiguration  to 
be  analyzed. 


CURRENT  APPROACH 


Mission -Loading  Definition 

The  mission  sf.- re-loading  rec.uirements  for  a new  aircraft  are  generaily  well  defined  in  the  speci- 
fications. Here,  the  various  roles  of  the  aircraft  are  detailed  and  the  armament  specified.  For 
example,  a strike  fighter  will,  typically,  perform  a filter  escort  mission  for  which  a few  relatively 
li^t  air-to-air  weapons  a.a  carried.  This  same  aircraft  might  well  be  required  to  also  perform  an 
interdiction  mission  carrying  a much  heavier  assortment  of  external  stores  and  a close -air-support  mis- 
sion with  a still  heavier  load  of  wapons.  If  all  the  possible  mission  and  weapon  requirements  of  the  air- 
craft are  set  forth,  one  soon  begins  to  realize  the  scope  of  the  problem  that  faces  the  flutter  analyst. 

The  enormous  number  of  possible  wing-store  combinations  is  the  central  problem  to  most  wing-store 
flutter  studies. 

The  first  job  of  the  flutter  analyst  is  to  reduce  this  seemingly  infinite  problem  to  one  of  manageable 
size.  Generally,  this  task  is  begun  by  discussing  mission  loadings  with  both  the  contracting  agency  and 
the  operational  squadrons.  In  the  case  of  tl  - . ,pr'thetical  strike  fighter,  it  is  not  unreasonable  to  expect 
that  out  of  these  discussions  will  emerge  si  eight  mission  loadings  which  must  be  thorou^ly  analyzed 
to  establish  stiffness  levels  necessary  for  proper  flutter  margins  throughout  the  defined  flight  profiles: 
that  is,  the  wing  should  be  designed  to  carry  these  basic  mission  loadings  and  their  subsets  to  specified 
limit  speeds. 

The  problem  then  becomes  one  of  minimizing  the  number  of  subsets  which  must  be  analyzed  to 
clear  an  entire  loading.  An  effective  way  of  accomplishing  this  reduction  is  to  make  use  of  inertial  sim- 
ilarity. This  procedure  is  based  on  the  establishment  of  a relationship  between  store  inertial  properties 
and  wing  flutter  characteristics.  This  is  accomplished  by  plotting  the  inertial  properties  of  all  possible 
stores  within  a mission  loading  at  a given  station,  as  shown  in  Figure  5.  Examination  of  this  plot  re- 
veals important  areas  where  analysis  of  a few  cases  is  representative  of  a larger  number.  Considera- 
tion of  plots  of  this  type  for  each  station  and  each  loading  v/ill  lead  to  a list  of  cases  which,  when  ana- 
lyzed, will  represent  the  entire  array.  This  list  becomes  the  basis  for  the  preliminary  analysis. 

Analysis 

A successful  theoretical  analysis  is  one  which  initially  identifies  problem  areas,  categorizes  them 
in  terms  of  impact  on  design  and  performance,  and  provides  detailed  relationships  between  stiffness, 
weight,  and  flutter  margin  to  the  designer.  The  accurate  identification  of  a wing-store  flutter  problem 
at  an  early  stage  of  development  is  essential  for  a minimum -impact  solution. 

A search  for  problem  areas  can  be  undertaken  as  soon  as  a tentative  strength  design  has  been 
determined.  The  clean  wing  is  analyzed,  of  course,  in  each  of  its  possible  modes:  cantilever,  sym- 
metric and  antisymmetric;  with  and  without  wing  fuel;  at  various  wing  sweeps,  if  movable;  and  at  several 
key  Mach  numbers.  The  wing  flutter  speed(s)  should  be  considered  in  terms  of  gradients.  If  flutter 
speed  changes  rapidly  with  respect  to  any  important  parameter  (fuel,  wing  sweep,  root  attachment  stiff- 
ness, outer  panel  stiffness,  altitude,  Mach  nuir.ber,  fuselage  inertia,  etc.),  the  phenomena  should  be 
studied  in  more  detail. 

As  a next  step,  candidate  stores  from  the  mission  loadings  described  above  are  introduced  into  the 
analysis,  one  station  at  a time.  The  effect  of  their  weight  and  inertia  on  the  clean  wing  should,  initially, 
be  determined  by  making  the  attachment  stiffness  infinite.  (A  typical  flutter  velocity  contour  is  shown  in 
Figure  6.)  The  modes  obtained  in  this  way  are  an  excellent  set  of  generalized  coordinates  for  later  use 
when  store  flexibilities  are  first  evaluated.  Combinations  of  rigidly  attached  stores  should  also  be  ana- 
lyzed using  ti  e mission  loadings  as  a guide.  Rather  extensive  modifications  to  the  clean-wing  flutter 
results  cai>  o expected,  especially  for  those  stores  most  outboard  on  the  surface.  As  above,  steep 
gradients  in  flutter  speed  as  a function  of  mass  or  inertia  at  a certain  wing  station  should  be  studied  and 
the  phenomena  properly  identified  and  understood. 

The  intvodiiction  of  pylon  flexibilities  is  next  undertaken.  It  is  suggested  that  cantilever  pylon 
modes  (velat've  to  the  fixed  wing)  be  coupled  to  the  modes  above  (wing  with  rigid  store  attachments) 
directly  in  the  flutter  determinant  during  these  early  design  studies  when  little  data  are  available  for 
precise  prediction  of  coupled  modes.  The  flv'icr  detern.mant  shown  in  Figure  7 can  be  configured  such 
that  it  is  necessary  only  to  change  a single  diagonal  element  to  change  a particular  flexibility;  if  store 
aerodynamics  are  ignored,  only  mass  coupling  exists  between  the  store  and  wing  coordinates.  Flutter 
trends  vs  pylon  flexibility  can  thus  be  obtained  very  efficiently  and  effectively.  The  hopeful  result  of 
these  trend  studies  would  bo  ranges  of  pylon  stiffness  where  acceptable  flutter  margins  exist.  If  satis- 
factory margins  cannot  be  predicted,  the  entire  analysis  should  form  a basis  for  reevaluation  of  t..e  sit- 
uation and  reconfiguration  of  the  design. 
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The  analyses  desciibed  to  this  paint  can  be  regarded  as  a stuify,  the  output  of  which  is  a descrip- 
tion of  the  flutter  behavior  of  a particular  aircraft  when  external  stores  are  introduced.  The  knowledge 
accumulated  to  this  point  should  be  sufficient  to  plan  the  remainder  of  the  theoretical  analysis  program. 

A relatively  small  number  of  cases  should  be  chosen  as  representative  and  a detailed  analysis  of  each 
performed.  These  cases  will  be  a basis  for  monitoring  the  design  as  it  progresses. 

Flutter-Model  Tests 

Having  completed  the  preliminary  theoretical  analysis,  a flutter-model  test  program  can  be 
planned.  An  early  goal  of  this  program  should  be  identification  of  the  transonic  flutter  characteristics 
due  to  the  presence  of  a store  shape  near  the  wing.  Component  models  tested  in  a small,  hi^-density 
blowdown  tunnel  are  useful  to  this  end.  (See  Figure  8.)  The  component  models  are  not  necessarily 
dynamically  similar,  but  should  have  node  lines  similar  to  full  scale  and  the  proper  thickness  distribu- 
tion. A vibration  test  is  performed  before  tunnel  testing  so  that  an  accurate  determination  of  generalized 
mass  and  stiffness  can  be  made  for  theoretical  analyses  of  the  models.  A clean-wing  flutter  boundary  is 
first  obtained  and  the  test  repeated  with  a very  lightwei^t  store  shape  added  to  the  model.  The  store  is 
kept  li^t  so  that  only  small  changes  in  the  wing  modes  occur,  ensuring  flutter  at  nearly  the  same  tunnel 
condition.  The  difference  between  the  two  boundaries  is  now  primarily  a measure  of  the  change  in  aero- 
dynamics due  to  the  presence  of  the  store. 

A large,  low-speed  flutter  model,  such  as  the  one  shown  in  Figure  9,  is  the  most  useful  experi- 
mental tool  to  determine  increments  in  flutter  speed  for  a large  variety  of  stores.  This  type  of  testing 
can  be  accelerated  greatly  ;?ith  the  use  of  a variable-inertia  store  whose  internal  wei^ts  can  be  post 
tioned  remotely  for  variations  in  store  center  of  gravity  and  inertia.  The  use  of  a low-speed  model  is 
also  an  excellent  way  of  assessing  the  influence  of  nonlinearities,  such  as  free  play,  on  store-flutter 
mechanisms.  Low-speed  models  should  be  cable-  or  rod-mounted  so  that  body  fre^oms  are  relatively 
uninhibited. 

Free-flying  transonic  flutter  models,  such  as  those  flown  on  cables  in  the  NASA  Langley  16 -foot 
Freon  Tunnel  (Figure  10),  are  best  suited  to  demonstrating  margins  in  the  ‘ransonic  regime  and  investi- 
gating special  problems  associated  with  Mach  number  and  altitude.  The  co.nplex  design  problems  and 
high  fabrication  and  test  costs  associated  with  these  models  do  not  make  them  attractive  for  tests  of 
more  than  a very  few  store  configurations.  t 

Ground  Vibration  Tests 

The  ground  vibration  survey  of  the  flight  article  is  traditionally  performed  just  before  first  fligjit. 

It  IS  costly,  time-consuming,  and  complex,  and  it  is  unreasonable  to  perform  complete  aircraft  vibration 
tests  for  each  external-store  configuration  of  interest. 

An  effective  way  of  alleviating  this  situation  is  to  perform  vibration  tests  of  isolated  pylon/store 
configurations  and  couple  these  measured  modes  analytically  to  either  calculated  or  measured  modes  of 
the  aircraft  without  stores.  A fixture  must  be  provided  which  properly  simulates  the  wing  pylon  attach- 
ment into  a fixed  base.  (See  F'.gure  11.)  These  measured  modes  are  comparable  to  those  used  in  the 
theoretical  study  described  above.  The  aircraft  modes  should  include  those  at  least  as  high  in  frequency 
as  wing  second  torsion  so  that  enough  wing  deformation  shapes  are  available  as  coordinates.  The  pylon 
store  modes  can  be  measured  whenever  the  stores  are  available  without  access  to  the  aircraft  itself. 

When  this  approach  is  employed , several  of  the  most  critical  cases  should  be  checked  by  performing 
ground  vibration  tests  on  the  fully  loaded  aircraft  to  demonstrate  the  accuracy  of  this  analytical  proc  ’- 
dure.  (Figure  12  shows  a heavily  loaded  aircraft  during  a ground  vibration  survey. ) 

Flight-Flutter  Tests 

It  IS  obvious  that  only  a relatively  small  number  of  store  configurations  can  be  flight -flutter  teote'’. 

To  this  end,  it  is  important  that  an  identification  be  made,  based  on  all  the  analyses  described  above,  of 
a few  configurations  which,  when  demonstrated  in  flight,  will  clear  all  remaining  configurations  of 
interest.  This  difficult  and  seemingly  idealistic  task  can  be  accomplished  if  the  flutter  mechanisms  are 
thorou^ly  understood  and  the  pertinent  parameters  have  been  monitored  throughout  the  design  process. 

Put  another  way , money  well  spent  to  support  the  engineering  analysis  is  money  saved  many  times  over 
in  flight  testing. 

To  make  the  fli^t  test  program  as  brief  as  possible,  it  is  recommended  that  a dedicated  control 
surface  be  installed  on  each  wing  to  act  as  an  exciter.  A surface  of  this  type  is  very  effective  in  the 
frequency  spectrum  of  interest.  Ample  excitation  of  the  stores  themselves  can  be  achieved  in  this 
manner.  Data  acquisition  methods  should  be  comparable  to  those  used  for  the  clean  aircraft,  e.g. , those 
described  in  Reference  1. 

RECENT  DEVELOPMENTS 

In  the  past  few  years  there  have  been  some  noteworthly  developments  in  several  areas  which  we 
believe  will  greatly  ease  the  burden  of  the  flutter  analyst  in  his  quest  for  a more  practical  solution  to  the 
wing-store  flutter  problem.  Generally,  these  developments  fall  into  three  main  categories:  active 
flutter  suppression,  fast  flutter  routines,  and  flutter  optimization. 
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The  most  promising  application  of  active  flutter  suppression  is  to  problems  of  wing  external-store 
flutter.  Reference  2.  A considerab  e amount  of  analytical  work,  wind  tunnel  testing,  and  flight  testing 
has  been  performed  in  the  United  States  on  this  approach.  Figure  13  is  a schematic  of  the  flutter- 
suppression  system  used  in  one  of  the  studies.  Active  flutter  suppression  is  attractive  for  store-flutter 
problems  for  three  main  reasons:  store  flutter  instabilities  are  commonly  very  mild,  which  suggests 
that  large  increases  in  flutter  speed  are  potentially  achievable;  a flutter-suppression  system  is  adaptable 
to  store  loadings  not  covered  in  the  original  aircr^t  specifications  ai?d,  accordingly,  not  covered  by  the 
original  flutter-prevention  program;  and  a wide  variety  of  store-flutter  mechanisms  can  conceivably  be 
suppressed  by  properly  adapting  the  electronically  compensated  feedba.tk  network. 

In  the  field  of  fast  flutter  routines,  two  promising  methods  have  been  developed;  one  by  McDonnell 
Douglas  (Reference  3)  and  one  by  Northrop  (Reference  4). 

McDonnell  Douglas  has  developed  a high-speed  wing-store  flutter  analysis  program  which  features 
a data  storage  and  retrieval  system  and  diagnostic  and  interpolation  subprograms.  A classical  vibration 
and  flutter  analysis  routine  which  has  been  configured  to  treat  wing-store  cases  in  an  optimum  fashion  is 
the  central  program  in  the  group.  Results  from  this  program  are  permanently  stored  for  later  use. 

When  a significant  amount  of  data  have  been  stored,  the  diagnostic  program  is  used  to  provide  information 
such  as  minimum  flutter  speed  or  flutter-speed  gradient  vs  some  pertinent  parameter.  When  a new  store 
is  to  be  analyzed,  the  interpolation  program  is  used  to  estimate  the  flutter  speed  instead  of  perJorming  a 
complete  analysis.  The  entire  system  is  tied  to  a cathode-tube  display  and  real-time  control  panel  pro- 
viding the  flutter  engineer  with  a powerful  analysis  tool. 

Northrop  engineers  Cross  and  Albano  have  shown  that  computing  time  can  be  reduced  considerably 
by  usuig  perturbation  theory  to  find  flutter  speeds  of  wing -store  combinations.  The  savings  are  accom- 
plished by  avoiding  the  repetitive  and  time-consuming  operations  of  extracting  eigen  data  from  real 
dynamic  matrices  (vibration)  and  complex  flutter  determinants.  A representative  case  is  first  solved  by 
classical  methods;  a new  store  mass  (or  flexibility)  representation  is  then  chosen,  which  is  expressed  as 
a'  incremental  change  in  the  dynamic  matrix;  the  roots  and  vectors  for  the  new  case  are  computed 
directly  in  terms  of  the  original  values  and  the  uicremental  dynamic  matrix  using  elementary  algebraic 
relationships;  an  incremental  flutter  determinant  is  then  derived  from  the  new  vectors,  etc. , and  the 
perturbation  equations  are  applied  again.  Northrop  has  been  able  to  solve  1000  analyses  per  hour  on  an 
IBM  370/165  computer  with  this  method. 

These  recent  developments  in  wing-store  flutter  prevention  have  been  directed  mainly  toward  in- 
creasing the  efficiency  of  hutter-speed  computations  for  variations  in  design  parameters,  the  so-called 
fast-flutter  techniques.  Altliough  these  techniques  aid  significantly  in  establishing  flutter-speed  ^store- 
parameter  trends,  definition  of  the  most  critical  design  cases  remains  based  largely  on  interpolation  of 
the  computed  data.  This  implicit  dependence  on  the  extensive  "mapping"  of  flutter-speed  "contours" 
leads  to  considering  the  feasibility  of  utilizing  a more  direct  numerical  search  procedure.  The  relatively 
recent  derivation  of  an  exact  expression  for  the  derivative  of  flutter  velocity  (Reference  5)  with  respect 
to  design  parameters  leads  to  the  possioility  of  employing  a flutter-velocity  gradient  search  to  determine 
the  critical  (minimum  flutter  speed)  configurations.  Pylon  stiffness,  store  mass,  and  store  inertia 
would  be  considered  as  continuous  design  variables  within  the  prescribed  envelope  of  permissible  values. 
Pylon  stiffness  parameters  would  be  constrained  not  to  fall  below  their  strength -required  values.  The 
gradient -search  approach  would  be  used  to  determine  directly  a minimum  flutter  speed  condition  without 
the  requirement  for  extensive  contour  mapping.  For  cases  where  more  than  one  critical  root  exists  for 
a given  flight  condition,  the  gradient  search  would  necessarily  be  applied  to  each  root  separately  to  find 
the  respective  worst  cases. 

Having  somehow  arrived  at  a critical  wing-store  configuration,  the  next  logical  step  in  the  design 
process  is  that  of  determining  a "fix"  to  achieve  a specified  flutter-speed  margin  for  selected  critical 
cases.  Up  to  new,  this  process  has  tended  to  rely  rather  heavily  on  engineering  judgment.  An  automated 
procedure  (Reference  6)  has  been  developed  recently  for  resizing  a structure  to  meet  both  flutter  and 
strength  requirements  for  minimum  weight  increase.  Although  the  reported  method  is  currently  being 
applied  to  structures  having  a single  critical  flight -flutter  condition,  it  could  be  extended  to  multiple 
conditions,  more  typical  of  the  store  flutter  problem,  by  holding  the  increased  gage  sizes  resulting  from 
redesign  for  one  critical  case  as  minimum  gage  for  redesign  for  any  subsequent  case. 

CONCLUDING  REMARKS 

We  have  attempted  to  show  how  the  problem  of  aircraft/external -store  flutter  can  be  handled  in  a 
realistic  and  practical  manner  and  how  the  magnitude  of  the  problem  can  be  reduced  to  manageable 
proportions.  We  have  also  touched  on  what  we  feel  is  the  current  approach  in  the  United  States  to  the 
design  of  aircraft  which  are  required  to  carry  a large  assortment  of  external  stores.  Finally,  some 
remarks  on  several  promising  new  developments  have  been  presented  which  provide  the  authors  with 
some  hope  that  in  the  future  aircraft  designers  will  be  able  to  cope  with  this  formidable  problem  in  a 
more  sophisticated  manner. 
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SUMMARY 

This  paper  first  presents  a sunmary  or  overview  gleaned  from  a brief  review  of  some  of  the  literature 
on  the  practiial  aspects  of  wing-store  flutter  prediction  and  prevention.  Brief  comments  are  given  on  the 
advantages  and  disadvantages  of  various  aspects  of  analytical  and  test  procedures.  Descriptions  of  improved 
analytical  procedures  developed  for  the  United  States  Air  Force  is  then  given.  Two  methods  are  described 
in  some  detail  and  the  results  of  the  investigators  are  outlined.  One  is  a rapid  special  purpose  wing- 
store  flutter  analysis  program  called  F'aCES.  It  has  data  storage  and  retrieval  capabilities  which  together 
with  a diagnostic  and  interpolation/extrapolation  procedure  estimate  the  flutter  speed  of  new,  similar 
stores.  The  system  can  be  coupled  to  a cathode  ray  tube  to  increase  man/machine  interaction  and  reduce 
decision  times.  The  other  analysis  metnud  described  is  based  on  the  perturbation  approach.  Computation 
times  can  be  reduced  90%  by  using  the  previously  available  data.  The  method  produces  good  results  when 
the  mass  or  stiffness  changes  are  small  so  that  in  turn,  eigenvalue  and  eigenvector  changes  are  small.  A 
graph  of  flutter  speed  versus  important  parameters  can  be  produced  in  one  minute  on  a modern  comruter. 
Limitations  of  the  method  are  discussed. 


INTRODUCTION 

Modern  capabilities  in  flutter  analyses  and  subsonic  and  transonic  wind  tunnel  flutter  model  techniques 
are  used  in  the  mid-term  design  period  to  determine  the  optimum  location  of  pylon  mounted  jet  engines  on 
transports  from  a fl.itter  avoidance  standpoint.  The  chordwise  and  spanwise  locations  of  the  propulsion 
systems,  and  natural  frequencies  of  their  suspension  systems,  are  selected  to  provide  higher  flutter 
speeds  and  avoid  valleys  of  low  flutter  speeds.  Such  action  is  possible  for  a select  critical  number  of 
"external  stores"  on  fighters.  However,  this  procedure  is  not  possible  for  the  vast  number  of  externally 
mounted  stores  or  weapons,  since  many  other  technical  considerations  force  placement  of  stores  on  fighter 
wings  at  positions  which  adversely  affect  flutter  characteristics.  Fighter  f. utter  speed  placards  on  many 
stores  are  frequently  in  the  speed  range  near  H = 0.8  and  thereby  restrict  and  constrain  the  speed-altitude 
mission  envelope  which  otherwise  might  be  available  to  a field  commander  (Reference  1).  Another  adver.e 
effect  of  the  performance  limiting  placard  is  a noticeable  reduction  in  survivability.  Some  relative 
evaluations  indicate  that  moderate  increases  (12-15S)  in  speed  near  M = 0.8  at  low  altitude  can  almost 
double  the  su> vivability  against  anti-aircraft  guns. 

Another  serious  consideratiun  for  the  engineer  concerned  with  the  prediction  and  prevention  of  flutter 
of  wings  carrying  external  stores  is  tne  potentially  overwhelming  magnitude  of  different  store  combinations 
that  could  be  carried.  It  would  not  be  unusual  to  be  concerned  with  several  thousand  (10,000)  realistic 
cases  In  analysis  which  are  further  evaluated  and  narrowed  down  by  subsonic  model  tests  to  several  hundred 
more  critical  cases.  These  ca,es  in  turn  are  reduced  further  to  the  most  critical  which  are  investigated 
by  transonic-low  supersonic  speed  model  tests  and  additional  final  analyses.  A limited  number  of 
representative  stores  exhibiting  the  main  highlights  are  selected  for  flight  flutter  checkouts  to  sub- 
stantiate safety  and  prove  out  the  reliability  of  the  process.  Store  inertia-flutter  mode  characteristic 
similarities  are  sought  out  to  reduce  the  number  of  individual  store  investigations.  However,  the  process 
is  extensive,  time  and  labor  consuming  and  expensive. 

The  above  briefly  summarizes  the  balanced  analytical-flutter  model  approaches  used  by  most  of  the  US 
contractors  for  wing-store  flutter  prediction  and  prevention.  The  spectrum  of  store  loadings  and  flutter 
critical  store  loadings  can  be  reduced  by  early  discussions  with  the  Air  Force  engineering  and  procurement 
agencies,  the  aircraft  manufacturer  and  the  using  Air  Force  command.  This  includes  a discussion  of  best 
oylon  and  suspension  for  critical  stores  and  store  release  sequences  to  minimize  speed  restrictions.  The 
dialogue  should  be  maintained  throughout  the  design  and  early  life  of  the  aircraft.  This  obvious  type  of 
action  must  be  given  priority  and  contacts  must  be  renewed  to  arrive  at  the  store  (mass-inertia)  loading 
spectrum  possible  on  each  oylon  and  rack. 

Another  essential  aspect  is  a complete  understanding  of  the  physical-mathematical  model  for  the  subject 
aircraft  with  stores  end  t!ie  cause  for  the  flutter  conditions  as  a function  of  the  store  variations. 
Exploratory  zero  airspeed  vibration  and  flutter  calculations  covering  a reasonable  part  of  the  store  mass- 
inertia  matrix  are  absolutely  necessary  to  locate  modes  which  could  couple  due  to  frequency  proximity 
(frequencies  versus  store  parameter)  and  to  obtain  experience  about  the  relative  magnitudes  of  mass  (or 
elastic)  couplings  which  induce  the  flutter  causing  wing  twist  and  also  the  subsequent  adve>'se  aero- 
dynamic couplings.  These  calculations  will  then  shed  light  on  the  more  critical  conditions  and  further 
calculations  needed.  They  also  can  confirm  or  refute  "rule  of  thumb"  statements  such  as  "locate  heavy 
stores  in-board  and  forward",  which  could  produce  more  critical  conditions  for  some  aircraft. 

This  introduction  ends  by  summarizing  the  highlights  of  a brief  but  not  encompassing  review  of  some 
literature  on  engineering  irediction  and  prevention  of  flutter  of  wings  with  stores.  The  overview  includes 
some  recent  comments  by  a few  specialists  and  is  given  in  the  following  Table  I.  One  important  new  develop- 
ment is  the  potential  use  of  act-ve  controls  to  suppress  wing-store  flutter  (Reference  1).  Comments  on 
recent  United  States  Air  force  research  on  wing-store  flutter  analyses  then  follow  in  subsequent  sections. 
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TABLE  I.  COMMENTS  ON  THE  HIHG-STORE  aUHER  PREDICTION  PROCESS 


Coiiiputer  Methods 


Inertia  Characteris- 
tics of  Inert  Bon4)s 


COMMENTS  1 REFERENCE 


I 

Host  analyses  are  conducted  with  digital  coitputers. 

General  Dynanics  uses  the  Direct  Analog  Electric  Analog  Computer 
also  in  a parametric  approach  to  obtain  flutter  boundaries 
rapidly  and  economically  in  terms  of  inertia  versus  weight  for 
various  c.g.  conditions  and  velocities 

LTV  employed  a data  storage  and  search  systein  to  3 

evaluate  different  stores  on  the  basis  of  inertial  similarity. 

A library  stores  wind  tunnel  and  flight  data.  Decisions  are 
based  on  the  evaluating  engineer's  experience  and  the  location 
of  the  store's  characteristics  within  an  acceptable  inertial 
similarity  field. 

Height  and  inertia  can  vary  18S.  Variations  should  be  4 and  5 

covered  in  analyses.  Flutter  speed  could  be  very  sensitive  to 
and  drop  sharply  with  a moderate  change  in  bomb  weight  for 
flutter  modes  involving  rack  bending,  store  pitch,  and  wing 
twist. 


Fuselage  Dynamics 


Pylon-Store  Dynamics 


The  fuselage  is  sometimes  relatively  stiff  so  coupling  is  Personal 
not  large.  However,  in  some  cases,  fuselage  frequencies  are  in  Comment 
the  range  of  wing-store  flutter  frequencies.  The  effective  wing 
bending-to-torsion  frequency  ratio  or  root  impedance  is  changed. 

Also,  fuselage  deformations  produce  a wing  angle  of  attack  at 
the  root  and  cause  aerodynamic  coupling.  This  is  very  important 
for  low  aspect  ratio  wings  with  long  root  chords  and  can  pro- 
duce fuselage  induced  wing  pitch  and  wing  bending  flutter. 

Active  flutter  control  could  be  a good  prevention  approach  in 
such  case.  Detailed  finite  element  methods  are  sometimes 
necessary  to  predict  fuselage  and  car^y-through  stiffness  and 
dynamics. 

This  is  the  prime  cause  of  flutter  through  production 
of  mass  coupling  and  adverse  frequency  ratios.  Store  motions 
including  pitch  produce  wing  twisting  and  aerodynam'c  coupling. 

Critical  frequency  ratios  exist  where  bending  type  wing  modes 
coalesce  with  higher  torsion  mode.  Pylon  or  rack  modes  change 
effective  critical  wing  bending  torsion  frequency  ratios. 

Flutter  modes  can  occur  which  seem  to  be  mostly  confined  6 

to  the  stores.  Modes  are  influenced  by  store  aerodynamics  and 
the  use  of  fins  is  beneficial  for  preventing  flutter. 

A good  definition  of  the  store-pylon  dynamics  including 
wing  connections  is  essential  but  can  be  complicated  for  higher 
complex  rack  modes.  Ground  vibration  and  stiffness  tests  and 
analytical  correlations  are  essential,  but  cause  a problem  in 
that  test  data  cannot  be  obtained  early  for  analyses  and  model 
tests.  Wing,  pylon,  fuselage  and  carry-through  stiffness 
coefficients  can  be  approximated  early  by  finite  element  models. 

Cantilevered  store-pylon  tests  can  be  conducted  to  7 

determine  resonant  frequencies  and  back-calculate  influence  3 and  8 

coefficients. 


Model  Tests 


Proper  dynamic  simulation  of  tip  missiles  and  launcher 
rails  through  similar  tests  are  required.  Tip  modal  deflections 
are  important  for  proper  estimation  of  aerodynamic  coupling. 

Cantilever  pylon-store  tests,  while  useful,  can  be 
influenced  by  too  much  flexibility  in  short  wing  section-sturdy 
jig  tests.  Centerline  pylon-store  modes  were  developed 
successfully  from  vibration  tests  with  wing  and  fuselage  on 
jacks,  however.  Pylon  and  rack  stiffness  characteristics  can 
be  non-linear  (buckling).  Effective  fuel  inertia  characteris- 
tics can  be  approximated  by  conducting  tests  on  empty  and 
partially  full  fuel  tanks  and  using  the  stiffness  data  determined 
from  the  empty  tank  vibration  tests. 

Both  subsonic  and  transonic  tests  are  employed  generally 
with  analyses  for  balanced  approach.  Lower  frequency  fuselage 
degrees  of  freedom  are  included  as  well  as  lower  mode  pylon  or 
rack-store  freedoms. 
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TABLE  I.  COHHEHTS  ON  THE  HING-STORE  FLUHER  PREDICTION  PROCESS  Cont=d) 


TOPIC COWHENTS REFERENCE 

Special  tanks  or  shells  are  used  together  with  variable 
or  irovable  masses  to  cover  inertia-mass  range  of  stores.  This 
may  not  be  a sufficiently  detailed  representation  of  bombs  on 
Multiple  Ejection  Racks.  Perisho  mentions  that  a shell  model  6 

with  coarse  variations  of  inertia  did  not  produce  acceptable 
results.  Use  of  fore  ard  aft  bombs  on  model  MER  rack  produces 
mild  store  flutter  mode  involving  S shaped  lateral  bending  ode 
of  r.^ck.  Bomb  nose  changes  and  bomb  fins  eliminate  flutter. 

Thus  higher  order  rack  modes  and  store  dynamics  can  be  important. 

Following  initial  model  program  with  extensive  analysis  to  3 and  8 
define  flutter  modes  which  have  gradual  loss  in  damping  and  the 
symmetry,  LTV  uses  half-span  symmetric  model  as  main  source  of 
flutter  data.  Full  span  tests  showed  flutter  velocity  is 
extremely  sensitive  (+40S  increase)  to  102  (+52)  asymmetrical  5 

(LH/RH)  difference  in  pylon  frequencies.  Also,  velocity  is  very 
sensitive  to  small  pylon  c.g.  (1/4  inch  on  aircraft)  asymmetry. 

Planform  symmetry  adjustments  improve  model  behavior  significant- 
ly also,  and  provide  some  different  results.  Pylon  dynamics  are 
critical  and  require  accurate  determination  from  stiffness  and 
ground  vibration  tests.  TER  flexibilities  modeled  from  ground 
vibration  test  results.  Rigid  TER  tests  are  conservative 
generally,  but  one  long  store  model  test  without  fin  shows 
unconservative  results  (no  flutter  on  rigid  TER).  Flexible  TER 
model  and  flight  tests  confirm  limited  rack-store  flutter.  LTV 
has  probably  now  evaluated  approximately  10,000  store  combina-  3 and  8 
tions  for  subject  aircraft.  Clean  wing  flutter  speed  well 
above  limit  dive  speed. 

LTV  later  conducts  half-span  subsonic  end  larger  transonic 
model  tests.  Due  to  sensitivity  of  flutter  to  store  rack  8 

parameters,  rack  influence  coefficients  are  obtained  from  rack 
vibration  data  rather  than  stiffness  tests.  Pylon  discrepancies 
on  previous  smaller  transonic  model  are  eliminated.  Low  speed 
and  transonic  model  test  results  agree  rather  well  with  flight 
data,  although  transonic  model  results  indicate  somewhat 
conservative  results. 

Store  Aerodynamics  Frequently,  store  aerodynamics  are  not  included  in  analyses  10 

and  lowly  damped  modes  are  evaluated  with  some  apprehension 
pending  test  data.  However,  subsonic  unsteady  aerodynamic 
methods  now  exist  for  reasonable  approximation  of  pylon-store 
aerodynamics.  Effects  are  expected  to  be  small  in  many  cases. 

Tip  missile  unsteady  aerodynamics  are  difficult  to  predict  7 

but  are  important.  Flow  separation  effects  at  transonic  speeds 
could  be  important  but  are  difficult  to  approximate  Calcula- 
tions for  various  wing  mounted  store  and  tip  missiles  do  not 
produce  consistent  results. 

Some  flutter  modes  occur  which  seem  to  be  confined  to  6 

motion  of  stores  and  are  influenced  by  store  aerodynamics.  3 

Use  of  fins  is  beneficial. 

Interference  was  suspected  between  fuselage  nacelle  ducts  9 

and  inboard  pylon-store.  Flow- through-nacell e-duct  might 
change  flutter  speed  of  inboard  pylon  on  transonic  model  and 
therefore  was  employed.  Proper  simulation  of  fuselage-wing 
cavity  seals  is  important. 

Store  Loadings  Individual  pylon  loadings  do  not  produce  independent 

flutter  effects.  Frequency  proximities,  mass  coupling  and 
aerodynamic  coupling  can  produce  widely  varying  results  for 
different  airplanes 

Some  comments  made  by  S.  Schwartz  for  one  aircraft  are:  7 

1.  Careful  attention  is  required  for  light  tip  stores  with 
heavy  inboard  stores.  Results  are  sensitive  to  the  centerline 
store  carried. 

2.  High  flutter  speeds  result  with  heavy  tip  stores  with 
forward  eg  and  with  inboard  stores. 

3.  Analyses  are  less  likely  to  accurately  predict  tip 
store  and  outboard  store  cases  than  other  cases. 
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TABLE  I.  COMMENTS  ON  THE  WING-STORE  FLUTTER  PREDICTION  PROCESS  (Confd) 


TOPIC 

COMMENTS 

REFERENCE 

4.  Flutter  speeds  with  all  stores  loaded  were  most 
difficult  to  predict  but  satisfactory  ’•arget  or  placard  speeds 
could  be  established  based  on  mission  performance  requirements 
and  lower  store  limitations  and  speeds. 

Coupling  with  Control 

Heavy  stores  or  moderate  to  light  tip  stores  reduce  the 

11 

System 

natural  frequencies  of  the  wing-store  system  significantly  so 
that  coupling  with  the  natural  frequencies  of  the  control 
system  are  possible.  Brock  covers  a static  or  mechanical 
coupling  problem.  However,  as  is  well  known,  dynamic  coupling 
is  possible  and  needs  to  be  predicted  and  prevented.  Aeroservo- 
elastic  prediction  methods  are  available  and  are  not  discussed 
here.  A closely  related  and  reverse  problem  to  aeroservoelasti- 

city  is  active  flutter  control  where  a high  frequency  SAS  system 
is  used  a priori  to  advantage. 

1 

THE  FACES  HING-STORE  DIGITAL  COMPUTER  - CRT  SYSTEM 


Fast  computer  programs  are  required  to  reduce  the  cost  and  time  of  extensive  analyses  for  very  large 
numbers  of  stores,  to  determine  the  more  critical  stores  and  their  flutter  speeds,  and  to  reduce  number 
and  costs  of  f'  tter  model  tests  and  flight  tests.  One  such  research-development  program  was  conducted 
for  the  Air  Force  Flight  Dynamics  Laboratory  (AFFDL)  by  Ferman  et  al.,  of  the  McDonnell  Aircraft  Company 
and  is  reported  fully  in  Reference  12  and  summarized  in  Reference  13.  Later  figures  are  excerpts  from 
these  documents. 

In  addition  to  the  development  of  a fast  flutter  routine,  the  AFFDL  considered  that  R&D  was  needed 
for  systematic  accumulation  and  retrieval  of  the  large  quantities  of  data  which  were  too  great  for  one 
or  two  individuals  to  recall  and  use  from  human  memory.  Furthermore,  difficulties  were  frequently 
caused  by  transfer  of  key  flutter  engineering  personnel  when  the  aircraft  was  in  main  use  leading  to 
delays  and  difficulties  in  flutter  decisions  and  considerable  repeating  of  previous  indoctrination  and 
training.  Other  R&D  requirements  specified  by  AFFDL  were  the  development  of  an  interpolation-extrapo- 
lation procedure  for  new  but  similar  stores  and  investigation  of  the  potential  of  man-machine  interfaces 
using  the  cathode  ray  tube  and  associated  controls.  The  data  storage  retrieval,  diagnostic-interpolation- 
extrapolation,  and  CRT  control  developments  were  of  moderate  scope  in  view  of  the  resource  limitations. 
However,  they  proved  successful.  The  CRT  aspect  proved  very  valuable  for  a priority  flutter  prevention 
study  on  a different  problem  for  a new  design. 

The  above  requirements  resulted  in  a computer  system  called  FACES  (Flutter  of  Aircraft  with  External 
Stores).  The  procedures  have  been  checked  out  and  the  computer  program  is  being  employed  extensively 
now  by  several  USAF  contractors  and  USAF  organizations.  Well  over  5000  store  configurations  have  been 
recently  investigated  using  this  system  for  several  aircraft.  This  program  or  system  is  outlined  on 
Figure  1. 

The  main  characteristics  of  the  Fast  Flutter  Routine  are  outlined  in  Table  II.  The  wing  simulation 
is  based  on  sectionalized  chordwise  rigid  strips.  This  beam  procedure  has  been  found  to  be  generally 
adequate  for  f ghters.  However,  the  program  can  also  accept  mode  shapes  and  frequencies  from  external 
inputs  in  case  chordwise  deformations  are  present  on  lower  aspect  ratio  wings.  The  incorporation  of  a 
mini-finite  element  stiuctural  analysis  module  is  a future  possibility.  The  Fast  Flutter  Routine  (FFR) 
can  use  subsonic  strip,  subsonic  lifting  surface,  or  supersonic  piston  theories.  Other  aerodynamic 
information  can  be  inserteo  from  other  external  sources.  Notice  that  the  program  ran  handle  five  pylons 
per  side.  Each  pylon  is  represented  as  a rigid  body  restrained  at  top  and  bottom  with  six  degrees  of 
freedom,  whereas  each  Multiple  Ejection  Rack  (MER)  has  six  degrees  of  freedom  to  represent  three  rigid 
sections  which  are  elastically  connected. 

Figure  2 shows  the  Data  Storage  and  Retrieval  System  flow.  It  accepts  data  internally  from  the  Fast 
Flutter  Routine  for  the  four  modes  with  the  lowest  flutter  speeds  or  from  other  external  data,  eithe" 
theoretical  or  experimental.  The  data  control  is  accomplished  by  KEYS  which  are  groupings  of  five  two- 
digit  words.  The  ten  digit  words  code  five  physical  properties.  The  number  of  words  required  increases 
with  the  number  of  stores  carried.  If  the  key  code  matches,  stored  data  information  is  listed  as  a 
result  of  the  retrieval  request.  Otherwise,  DSRS  can  be  used  in  a further  process  with  diagnostic  and 
interpolation-extrapolation  procedures  described  next. 

The  Diagnostic  Process  (DIAPR)  and  the  Calculation  Interpolation  Decision  Proces*’  (CID)  operate  as  a 
combined  process.  Table  III  shows  the  sensitivity  information  used  as  engineerino  reminders  and  for 
aiding  engineering  evaluations.  Sensitivity  criteria  are  built  into  the  program.  If  a new  store  satisfies 
these  criteria,  then  CID  uses  these  stored  data  to  estimate  the  similar  but  new  store's  flutter  speed. 

The  user  is  notified  regarding  satisfaction  of  acceptance  or  sensitivity  criteria.  Figure  3 is  a carpet 
type  diagram  and  outlines  a limited  but  succe.sful  application  of  the  data  storage  and  retrieval  system. 

The  various  cases  show  that  stored  data  are  ret. ’eved  and  listed  when  data  matches  exist.  Two  cases 

show  no  data  are  retrieved  (paths  4 or  6)  if  there  is  .no  parameter  match.  If  no  exact  match  results, 
the  combi.ned  DIAPR/CID  module  can  be  asked  to  delete  a variable  and  perform  a trend  study.  Thr^e  or  lore 
matches  from  a trend  study  are  required  to  estimate  flutter  speeds  of  new  stores.  Simple  (parabolic) 

curve  fittings  are  used  for  estimates.  Extrapolations  are  acceptable  if  the  parameter  is  within  20*  of 

maximum  or  minimum  values  of  the  parametric  stores  data  arid  if  flutter  velocity  derivative  data  or 


acceptance  criteria  are  within  specified  ranges  (velocity-damping  crossing  is  not  too  shallow  or  velocity 
does  not  vary  too  rapidly  with  the  parameter).  The  260  knot  estimated  flutter  speed  at  a pitch  inertia 
of  ? X 10®  lb  inches”  and  600  lbs  was  extrapolated  but  was  questioned  by  the  acceptance  criteria  because 
the  pitch  inertia  of  the  store  was  well  outside  the  20%  data  range  criterion. 
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FIGURE  1.  BLOCK  DIAGRAM  FOR  FACES  DIGITAL  COMPUTER  PROGRAM 


TABLE  II:  OUTLINE  OF  FAST  FLUTTER  ROUTINE  CHARACTERISTICS 


CANTILEVER,  SYMMETRIC,  ANTISYMMETRIC 

W WING  SECTIONS  - SERIALLY  KINKED  ELASTIC  AXIS 

STREAiiWISE  MASS,  ELASTIC  AXIS  MASS 

TWO  WING-ROOT  FREEDOMS,  STIFFNESS  OR  INFLUENCE  COEFFICIENTS 

WINS  El  AND  GJ  STIFFNESS  OR  INFLUENCE  COEFFICIENTS 

FIVE  PYLONS  PER  SIDE  - - MER,  TER,  SINGLE  STORES 

5 FREEDOM',  SPRINGS,  INFLUENCE  COEFFICIENTS  FOR  EACH  PYLON 

6 FREEDOMS,  SPRINGS,  INFLUENCE  COEFFICIENTS  FOR  EACH  MER  RACK 
SAVE  WING  MASS,  STIFFNESS,  EFFECTIVE  MASS 

STRIP  THEORY,  LIFTING  SURFACE  AND  STORE  AERO,  PISTON  THEORY 
MULTICASE  OPTIONS:  (FLUTTER) 
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FIGURE  2.  DATA  RETRIEVAL  FLOH 
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TABLE  III,  DESCRIPTION  OF  DIAGNOSTIC  PARAHETERS 


-SYMBOL 

NAME 

FUNCTION 

'Vm 

MINIMUM  FLUHER 
SPEED 

TO  ESTABLISH  THE  MINIMUM  FLUTTER  SPEED 
SPEED  FROM  A FLUTTER  SPEED  TREND  VS  A 
PARAMFTFR 

«f)r 

“REQUIRED  PLUTTER 
SPEED 

USED  TO  ESTABLISH  FLUTTER  CLEARANCE 

i^yF>A 

ASYMPTOTIC  FLUTTER 

TO  ESTABLISH  THE  FLUTTER  SPEED  ASYMPTOTE 

^^F^A 

ASYMPTOTIC  <; 

TO  ESTABLISH  THE  PARAMETER  VALUE 
WHERE  aV/ac- 

aV/ag 

FLUTTER  SENSITIVITY 

TO  ESTABLISH  FLUTTER  SPEED  GRADIENT  WITH 
DAMPING 

(aVAO* 

FLUTTER  GRA.DIENT 

TO  ESTABLISH  FLUTTER  SPEED  GRADIENT  WITH 
CHANGES  IN  A PARAMETER 

fc 

FLUTTER  CLEARANCE 
FACTOR 

TO  ESTABLISH  FLUTTER  SPEED  CLEARANCE  IN 

PERCENT.  EQUALS 

C(Vf)h  . (Vf)r]  (100)/(Vp)R 

ACT 

ACCEPTANCE  CRITERIA 

EVALUATES  FLUTTER  SENSITIVITY  BEHAVIOR  AND 
ESTIMATED  FLUTTER 

“o-Vj-P/Vj 


u 


i 


fi 

A3 
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FIGURE  3.  DSRS  RESULTS  FOR  FIGHTER  WITH  POD  AT  60*.  SEMISPAN 

One  version  of  FACES  is  coupled  to  a cathode  ray  tube  system  to  increase  interface  between  engineer  and 
the  computer,  control  computations,  accomplish  rapid  variations  of  parameters,  display  results  and 
accelerate  turn-around  time  end  engineering  decisions.  Flutter  data  from  the  CRT  mode  of  operation  are 
incorporated  into  the  storage  system. 

Ihe  CRT  system  can  display  a list  of  natural  frequencies  and  as  shown  in  Figure  4,  also  display  mode 
shape  information.  Such  data  should  be  valuable  in  indicating  the  likelihood  of  flutter  due  to  proximity 
of  frequencies  of  important  modes  and  the  amount  of  response  coupling  oetween  wing-pylon-store  components 
in  a mode  as  well  as  the  amount  of  relative  wing  twisting.  The  latter  could  be  an  indication  of  aero- 
dynamic coupling.  The  CRT  user  options  are;  Display  Eigenvectors,  Display  Deflections,  Plot  Deflections, 
Change  Input  Matrix  Data,  Change  HtR  Rack  Configuration,  Redisplay  Frequencies,  Rerun  Vibration,  and  Vibra- 
tion Analysis  Complete.  The  program  also  displays  flutter  mode  crossing  m terms  of  mode  number,  flutter 
speed  and  frequency.  Damping  versus  velocity  and  frequency  versus  velocity  trends  are  displayed.  Root 
sorting  is  employed  to  track  modes.  These  capabilities  are  shown  in  Figure  5 which  also  simulates  an  actual 
application  case.  Figure  6 shows  a simulation  of  an  actual  display  of  data  retrieval,  flutter  speed  estimates 
and  diagnostic  information.  The  acceptance  criteria  reject  the  estimates  because  the  required  flutter 
speed  IS  too  close  to  the  minimum  flutter  speed  for  the  parameter.  Also,  the  estimated  flutter  speed  could 
be  adversely  affected  by  damping  uncertainty  caused  by  a shallow  g-V  curve.  Figure  7 shows  some  comparisons 
of  FACES  data  with  data  from  other  methods  and  the  trend  of  flutter  speed  with  bomb  ejection  sequence. 

However,  it  could  be  an  illustration  of  a future  capability  in  machine/flutter  engineer  interaction  and  fast 
flutter  evaluations. 

As  previously  mentioned,  the  wing-store  flutter  analysis  system  called  FACES  has  been  checked  out  by 
many  comparisons  with  results  from  other  comfuter  programs.  It  is  in  extensive  use  today  at  several 
industrial  firms  and  US  Government  organizations,  and  is  an  extremely  valuable  tool  for  screening  many 
store  combinations  to  determine  satisfactory  or  critical  store  combinutions.  The  time  and  cost  savings  of 
the  FACES  program  have  been  evaluated  relative  to  a similar  method.  The  same  engineering  skill  ind  wing 
finite  section  approach  was  used  to  provide  a fair  comparison.  Figure  8 shows  that  a 2 to  1 reduction  in 
costs  and  core  occupancy  in  favor  of  FACES  results  when  carrying  two  stores  per  side.  The  advantage  increases 
to  3 to  1 for  four  stores  per  side.  Other  special  purpose  programs  could,  of  course,  reduce  the  advantage. 

Future  improvements  could  include  a mini-finite  element  structural  analysis  method  for  wing,  wing  root, 
pylon  and  rack.  Any  fast  finite  element  structural  model  should  avoid  extreme  detail  in  favor  of  speed. 

Other  future  possibilities  are  incorporation  of  low  supersonic  Mach  number  aerodynamics;  improved  pro- 
cedures in  data  storage-retrieval,  diagnosis,  and  interpolation  and  extrapolation;  incorporation  of  pitch 
inertia-mass  carpet  plot  capability  {displaying  points  or  contours  of  constant  velocity  in  the  inertia- 
mass  plane);  a multi-surface  fit  for  estimating  flutter  char  ‘eristics  for  new  stores  in  lieu  of  a curve 
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fit;  and  display  of  plots  of  store  deflections  relative  to  the  wing  on  the  CRT.  More  applications  of  data 
storage,  retrieval,  new  store  flutter  evaluations  and  the  CRT  capability  are  certainly  needed  and  some 
will  be  performed  by  AFFOL. 
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FIGURE  CRT  GRAPHICS  VIBRATION  DATA  FOR  FIGHTER  WITH  STORES  ON  TWO  PYLONS  PER  SIDE 
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FIGURE  5.  CRT  GRAPHICS  FLUTTER  DATA  FOR  FIGHTER  WITH  STORES  ON  TWO  PYLONS  PER  SIDE 
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FIGURE  6.  DISPLAY  OF  CRT  DATA  STORAGE-RETRIEVAL  AND  CALCULATION  INTERPOLATION  DECISION 
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FIGURE  8.  Ti«E  AND  COST  COMPARISON  GENERAL  FLUHER  PROGRAM  AND  FACES 
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THE  PERTUREATIOH  METHOD 


Another  approach  which  could  expedite  the  flutter  analyses  of  wings  which  carry  many  store  combinations 
is  the  perturbation  method.  The  potential  of  this  method  for  dynamic  analyses  has  been  discussed  by  van  de 
Vooren,  S>.rbin,  flax,  and  Siegel  and  Andrew  (References  14  to  19).  Cases  exist  where  changes  to  planform, 
inertia,  or  stiffness  parameters  are  sufficiently  small  so  that  corresponding  changes  in  the  flutter  and 
vibration  eigenvalues  and  eigenvectors  are  also  small.  Then  the  new  values  can  be  approximated  through 
perturbation  methods  using  a considerable  amount  of  the  previous  data  and  approximating  the  increments  with- 
out eigenvalue  extractions.  This  can  be  accomplished  with  considerably  less  work  and  cost. 

At  the  request  of  AFFOL,  Cross  and  Albano  of  the  Northrop  Corporation,  developed  a tnird  orde*"  perturba- 
tion method  for  predicting  the  vibration  and  flutter  characteristics  of  aircraft  with  wing  stores  (Reference 
20).  It  simulates  the  free  flight  condition  and  accounts  for  rigid  body  freedoms.  Perturbation  flutter 

analyses  are  made  only  if  the  zero  airspeed  vibration  modes  have  converged.  Second  order  terms  for  the 

eigenvectors  are  used  in  the  perturbation  vibration  and  flutter  analyses  to  estimate  effects  of  moderate 
changes  in  parameters.  The  third  ordei  flutter  eigenvalue  term  is  used  to  evaluate  convergence.  The 
perturbation  method  developed  was  applied  to  eight  series  of  wing-store  combinations  having  63  elastic 
degrees  of  freedom  and  3 "igid  body  freedoms  for  each  symmetry.  Up  to  four  underwing  stores  per  side  were 
used  in  some  applications  having  24  store-pylon  degrees  of  freedom. 

Some  of  the  results  obtained  by  Cross  and  Albano  will  now  be  presented  and  discusseo.  For  many  cases 
there  was  very  litMe  difference  between  tesults  from  exact  (classical)  methods  and  those  from  the  perturba- 
tion method.  However,  in  some  cases  involving  larger  perturoations  or  complications  cue  to  proximity  of 
modes,  the  perturbation  method  did  not  give  converged  results. 

Results  for  variations  in  the  chordwise  location  of  a tip  missile  are  shown  in  Figure  9.  The 

perturbation  method  results  are  in  excellent  agreement  with  corresponding  data  from  classical  methods. 

For  the  minus  20%  chordwise  location  of  the  tip  missile  eg  there  is  a switch  in  flutter  mode  compared  to 

the  basic  0%  chr-dwise  case.  Figure  10  shows  the  g-V  curves  for  both  cases  and  the  mode  switch.  Both  modes 

are  well  predicted  using  the  same  basic  case  as  a start.  Use  of  more  modes  (12  in  lieu  of  9)  improves  the 

agreement. 

Figure  11  shows  another  good  agreement  between  methods  for  the  case  of  tip  missile  inertia  variation 
(for  different  missiles)  and  when  an  outboard  store  is  carried.  Similar  good  agreement  is  obtained  for 
variable  tip  store  pitch  inertia  while  carrying  an  inboard  store. 

Coupling  between  inboard  and  outboard  stores  can  produce  cloie  proximity  of  baseline  frequencies  and  a 
lack  of  convergence  in  the  perturbation  techniques.  Smaller  perturbations  could  obtain  better  trends. 

Figure  12  shows  non-con/erged  results  when  baseline  vibration  frequencies  have  3%  separation.  The  160 
slug-ft2  pitch  inertia  symmetric  perturbation  case  compares  well  with  data  from  conventional  procedures. 
However,  there  is  a rapid  rise  in  flutter  speed  on  each  side  of  the  basic  case.  The  lack  of  convergence 
indicates  a steep  gradient.  The  antisymmetric  case  gives  somewhat  better  results  but  agreement  is  handi- 
capped by  changing  conditions  near  a pitch  inertia  of  110-120  slug-ft^.  However,  flutter  velocities  are 
extremely  high  and  well  above  practical  values. 
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In  connection  with  the  above  series  of  calculations,  one  perturbation  proceeded  from  a 2000  oound  borab 
at  the  centerline  to  a full  tank  at  the  centerline  with  the  same  weight  but  with  2.4  times  the  pitch  Inertia. 
This  produced  a lower  frequency  mode  differing  radically  from  those  baseline  modes  with  the  centerline 
borab.  Hot  unsurprisingly,  the  perturbation  method  produced  nonconverged  results. 

Figure  13  shows  relatively  good  agreement  for  the  case  of  tip  missile  with  outboard  store.  The  antl- 
synaetrlc  mode  8 Is  reasonably  well  predicted  ana  needs  to  be  located  to  assure  high  flutter  speeds  since 
It  Is  violent.  The  higher  pitch  Inertia  conditions  probably  produced  perturbations  which  were  too  large. 
Improvements  are  expected  with  smaller  perturbations. 

Figure  14  shows  the  difficulties  encountered  when  two  symmetric  modes,  the  fourth  and  fifth,  are  within 
1.45  In  frequency  so  that  convergence  is  poor.  Furthermore,  these  modes  coalesce,  possibly  switch  and  are 
lightly  damped,  adding  to  the  difficulty.  Mode  3 diverges  more  rapidly  and  Is  better  predicted.  The  anti- 
symmetric modes  are  well  separatee?  (55)  In  frequency.  An  agreement  between  perturbed  and  exact  results  Is 
generally  excellent  except  for  the  1005  pitch  Inertia  change.  Flutter  speeds  are  very  high. 
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FIGLRE  13.  COMPARISON  OF  RESULTS  FOR  TIP  MISSILE  AND  VARIABLE  OUTBOARD  STORE 


The  following  conclusions  are  drawn: 

a.  Wnen  the  changes  in  parameters  are  small  so  that  eigenvalues  and  eigenvector  changes  in  turn 
are  small,  the  perturbation  technique  yields  reliable  results. 

b.  A 905  savings  in  execution  time  is  possible. 

c.  The  perturbation  method  is  very  effective  in  screening  and  locating  critical  cases.  It  can 
also  be  employed  to  determine  best  pylon-store  dynamic  characteristics  (suspension  tu.ning). 

d.  A complete  flutter’  trend  versus  a controlling  parameter  (Figure  9)  can  be  obtained  in  one 
minute  using  a modern  computer  and  could  ba  displayed  on  a cathode  ray  tube  for  qu'ck  engineering  decisions. 

e.  Difficulties  are  caused  in  convergence  when: 

(1)  A flutter  velocity-frequency  curve  crosses  another  mode  and  flutter  mode  dampings  are 


(2)  A new  lower  and  radically  different  mode  is  introduced  by  a large  perturbation. 

(3)  The  baseline  frequencies  are  in  close  proximity. 
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FIGURE  14.  COMPARISON  OF  RESULTS  FOR  TIP  MISSILE  WITH  VARIABLE  INBOARD  STORE 


f.  Smaller  perturbations,  more  modes  and  more  reduced  frequencies  could  alleviate  some  of  the 
above  probless. 

g.  The  perturbation  method  can  be  a useful  and  economical  tool  for  flutter  analyses  involving 
changing  parameters  or  determining  optimum  designs. 
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